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ABSTRACT 


This report describes work done to, analyze and evaluate 
the attitude control performance of the Solar Electric Propulsion 
System (SEPS). A thrust vector control system for powered flight 
control is examined along with a gas jet reaction control system, 
and a reaction wheel system, both of which have been proposed for non- 
powered flight control. Comprehensive computer simulations of each-, 
control system were made and evaluated using a 30-mode spacecraft 
model . 


Results obtained indicate that thrust vector control and 
reaction wheel systems offer acceptable smooth proportional control. 

The gas jet control system is shewn to be risky for a flexible structure 
Such as SEPS, and is therefore, npt recommended as a primary control 
method. 
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SECTION 1 

INTRODUCTION AMO SUMMARY 

This report presents the results of a comprehensive study i 
undertaken to establish the attitude ^control performance provided by 
the Solar Electric Propulsion System (SEPS) for the planned International 
Comet Mission. 

1.1 BACKGROUND 

The "Halley Flyby/Tempel 2 Rendezvous (HFB/T2R) International 
Comet Mission" presents a rare opportunity to fly by the comet Halley 
enroute to rendezvous with a second comet (Tempel 2). This mission 
is possible by usinq a Solar Electric Propulsion System (SEPS) with a 
modest performance level based on current technology. An ecliptic 
plane projection of the overall heliocentric trajectory is illustrated 
in Figure 1-1 . 

The scientific objectives of the mission are the determination 
of the chemical composition and physical structure of cometary nuclei 
and comae, investigation of the dynamic processes occurring ^t The 
surface of nuclei and within the comae, and investigation of the inter- 
action of comets with the solar wind, including study of the formation 
of and structure of cometary tails. 

The mission plan calls for a single flight system to be launched 
by the space shuttle/inertial upper stage in the surnner of 1985. ►he 
flight system (Figure 1-2) is composed of three basic elements: 

1) THE SPACECRAFT, which Includes all rendezvous science 

electronics, booms, and platforms, will provide deployment, 
pointing, power, commanding, and d^ta collection for the 
science payload. It will also provide sequencing control 
for the probe and the SEPS. 

11) THE PROBE, which will be released at Halley, will contain 
Its own set of science instruments. 

ill) THE SEPS, which is a standard space transportation vehicle 
capable of carrying planetary or Earth orbital payloads 

^International Comet Mission: Halley/Tempel 2 Mission Baseline, 
NASA/ESA Doc jnent JPL 626-2, November 1979 (JPL Internal Docunent). 
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Figure 1-1. HFB/T2R Typical Transfer. Trajectory 






on missions for which ion propulsion is appropriate. It 
includes the mercury ion drive thrusters, mercury storage 
tanks, and two large (each 32 m long) solar pane!s 
which provide the energy for the electric propulsion. 

Besides power and propulsion, the SEPS provides flight 
system attitude control and maneuvering and contains 
associated celestial sensors and gyros. 

The fast flyby (57 km/s) encounter with Halley occurs only four 
months after launch, when Halley is approaching perihelion about 1.5 Al) 
solar ranye. The Halley probe is released from the spacecrai t to acquire 
data on a ballistic path through the atmosphere of Halley, while the 
rendezvous spacecraft is deflected by the SEPS for safe flyby to the 
sunward side of Halley. The spacecraft makes remote observations of 
Halley during flyby as well as serving as a relay link for the data 
acquired by tho probe and transmitted in rea' time to the spacecraft. 

The SEPS provides thrust virtual ly continuously throughout the 
heliocentric transfer portion of the mission, with thruster-off periods 
restricted to brief Intervals required for Halley probe deployment, 

Halley encounter data acquisition, and a few special engineering and 
science instrument calibration sequences during cruise. The number of 
thrusters operating and the thruster throttle level are programmed for 
optimum utilization of the available solar array power as the solar 
distance changes. 

Rendezvous with Tempel 2 requires a total flight t<me of nearly 
three years, and is scheduled to occur about two months before comet 
perihelion. The standard mission is scheduled to extend for a year after 
Initial rendezvous with Tempel 2 to permit extensive exploration of the 
cometary environment and nucleus until Tempel 2 Is 3.0 AU from the 
Sun. 

Attitude control and maneuvering throughout the mission is provided 
by the SEPS, both during powered flight to and around. Tempel 2, as well 
as during unpowered orbital operations about Tempel 2. 

Attitude control during pewered flight is achieved by gimballing 
of the ion drive thrusters under what will be called "thrust vector 
control." This control method is used virtually continuously during 
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the first 3 years of the mission while in the heliocentric transfer 
portion of the mission. After Tempel 2 rendezvous the ion engines 
must be turned off and will only be used infrequently (a few hours 
a week) to permit the desired circumnavigation of the nucleus. 

With the thrusters off most of the time following rendezvous, 
a second means of attitude control must be*provided. Two candidates 
that have been proposed are 1) gas jet reaction control, and 
2) momentum wheel control . 

1.2 STUDY OBJECTIVES AND SCOPE 

The SEPS/ICM Flight System represents a new kind of space 
vehi.de. What makes it different from past spacecraft, from the point 
of view of attitude control, is its shear size and flexibility. The 
2 large solar panels (required to provide the large amounts of energy 
needed to operate the ion thrusters) give the vehicle a wing span of 
~75 m (246, ft) and result in very large vehicle inertias (200,000 kg~m ). 
Available control torques to handle these large inertias are low. 
Constraints on system weight force the solar panels to be thin, 
lightweight structures* and, hence, quite flexible. 

The objectives of this study have been to understand the nature 
of the control problems associated with this large flexible vehicle and 
to evaluate the attitude control performance provided by the baseline 
SEPS for the International Comet Mission for the 3 control methods under 
consideration: 

1. Thrust vector control 

2. Gas jet reaction control 

3. Reaction wheel control. 

The baseline attitude control performance requirements for SEPS/ICM are 
given in ref. [1J. Those requirements are to provide a pointing of the 
vehicle within + 1° of a celestial reference direction with a maximum 
inertial rate not to exceed l°/s. 

The approach used has been one of time domain analysis 
(through digital computer simulation) of a comprehensive flexible 
dynamics model of the vehicle, augmented with the corresponding control 
laws and control hardware models. Each of the 3 control methods has 
been studied by simulation of the following four specific maneuvers: 

ll ^Ibid 
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1. Complete vehicle turn about the sun-line. 

2. Reorientd* 'on turn of the SEPS bus while simultaneously 
rotating the solar panels so as to keep them pointed 

to the sun, as the bus turns. 

3. Acquisition of small initial position errors and rates. 

4. Attitude control in the disturbance environment created 
by Science Scan Platform slewing'. 

Section 2 presents the model of SEPS/ICM along with a 
description of the spec i f ic maneuvers considered, and several other 
study assumptions. The systems description, and performance evaluation 
for' thrust vector, gas jet, and reaction wheel control are included in 
sections 3, 4, and 5. respectively. 

1.3 CONCLUSIONS ANO RECOMMENDATIONS . 

The main contribution of this study is that it establishes 
the Attitude Control Performance characteristics of SEPS for the 3 
proposed control methods: TVC. Gas Jets, and Moment Wheel Control. 

The performance characteristics obtained are given as maneuver perfor- 
mance plots in sections 3, 4, and 5. 

In addition to the performance characteristics, 
considerable insight and experience has been gained durinq this study, 
pertaining to the control of flexible spacecraft in general, and of 
the proposed flexible SEPS/ICM, in particular. Advanced simulation 
programs have been created which simulate SEPS flexible dynamics and 
constitute a sophisticated analysis and decision making tool inmediately 
available for further use, should it become desirable. 

A principal conclusion of this study is that both of the 
proposed TVC and Reaction Wheel attitude control (A/C) systems 
constitute acceptable control options for SEPS/ICM in that they both 

1. Are based on currently existing technology. Require 
no advanced enabling technology. 

2 . Meet baseline performance requirements. 
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3. Do so in a very smooth and satisfactory manner, 
without exciting vehicle flexible dynamics signifi- 
cantly and without exhibiting control/structure 
Interaction. 

4. (because of 3) both control methods are not judged 

to be very sensitive to SEPS structural uncertainties, . 
a very important advantage. 

For the 2 Gas Jet Reaction Control Systems considered, it 
is concluaed that 

1. One of the systems appears to meet performance require- 
ments, the other does not. 

2. Performance for both is characterized by high vibration 
levels of the structure with substantial excitation of 
the vehicle flexible appendages. These high vibration 
levels are due to the impulsive nature of the loads 

, produced by the nearly-square gas jet pulses. Such 
loads have a very high frequency content and, thus, 
tend to excite all resonant frequencies in the structure. 
Thfc higher vibration levels are very dangerous and, in 
fact, in 2 out of the 6 simulations, have been shown 
to lead to severe control /structure interaction 
resulting in exceedingly high limit cycle rates, and 
unacceptable propellant consumption. 

3. Gas jet RCS systems, as a class, tend to be highly 
interactive with the structures and their performance 
highly sensitive to model uncertainties. Performance 
will be reliably known only after launch, when it may 
be too late. 

4. Because of the above, RCS systems are judged very risky 
for SCPS, and therefore, are not recoomended as a primary 
control method. 

To further elaborate on the above, it must be understood 
that the attitude control problem of guaranteeing stability and 
performance of the SEPS with the proposed large, flexible solar arrays 
is a sizeable engineering challenge. 
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It has been our experience at JPL that flexibility (even 
slight) can cause serious problems. To name a few. Explorer I, Mariner 
Venus -Mercury; and Voyager, all exhibited flexible appendage effects 
which caused serious in-flight attitude control anomalies (unstable 
behavior, excessive gas consumption, etc.) and ultimately, resulted in 
costly and extensive "fixes" and/or mission limitations. All this 
happened in spite of considerable ground testing and efforts to 
accoimioda te model uncertainties. These S/C were only slightly 
‘flexible.- The SEPS on the other hand, is quite flexible. The solar 
arrays are complex non-linear structures and their large size, mass, 
and flexibility will have a predominant effect on the total vehicle 
dynamics. The arrays ape, in fact so complex and non-linear that even 
the best available structural dynamics modeling techniques yield 
models with uncertainty levels too high for reliable attitude control 
work. Thus, we have a paradox where the array characteristics are 
more difficult to model than for previous spacecraft, but their 
potentially disastrous effects on attitude control make accurate 
modeling more critical . 

Structural Vwiysts provided the solar array dynamics model 
used in this report. It is the best available model that can be 
developed using current state-of-the-art modeling techniques. More 
accurate dynamics models of the solar arrays are unlikely to become 
available in the near future because the development of such accurate 
models would require costly in-orbit dynamic testing 1 J of the solar 
arrays. 

In view of the above, it is clear that the designer of the 
A/C system must aim for a control method which is tolerant of model 
uncertainties. It is precisely in this area that gas jet systems 
show their bad side-effects. They are highly interactive with 
flexible structures, and hence not forgiving of model uncertainties. 

The risk associated with an RCS system in a highly flexible vehicle 
is simply too great. 

^ ivar, A.F., and Garba, J.A., "Attitude Control Considerations for 
Enhanced Instrumentation on SEPS Solar Array Flight Experiment," • 
Engineering Memorandum EM 347-30, 31 May 1979, Jet Propulsion 
Laboratory Internal Document. 
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TABLE 1-1. RCS & RWH Advantages and Disadvantages 
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1.4 


PROPOSAL FOR AN ALTERNATE CONTROL METHOD 


The purpose of this study has been to evaluate the 3 
proposed A/C methods, and not to generate new candidates. In the process 
of this evaluation, however, we have gained a considerable amount of 
new information which suggests that all of the objectionable drawbacks 
associated with a single gas jet or reaction wheel system, may be 
eliminated if we judiciously incorporate both of them into the design.. 
Table 1-1 presents a summary of advantages and disadvantages for 
individual RCS and reaction wheel systems. 

After examining Table. 1-1 in detail, it becomes clear 
that neither a Gas Jet RCS, nor a Reaction Wheel (RUH) system by 
itself will meet all non-powered flight A/C requirements in a satis- 
factory manner, nor will provide a reliable, low risk, control method. 

The following combination of both may provide an acceptable 
solution by taking the best from each method: 

1. Use Reaction Wheels for normal cruise momentum 
management and maneuvering. (This provides smooth 
and reliable control and considerable mass savings 
for turns and slewing.) 

2. Use a much smaller RCS system to provide high control 
reserve needed to handle Initial rate reduction and 
emergency situations, as well as to unload RWH when 
saturated. (This provides the reserve when needed.) 

The following additional option may Improve SEPS A/C 
problems substantially at very minor cost: 

. 3. Provide a center of mass/center of solar pressure trim by 
including a capability of realigning solar array axes 
(2 to 3°). This will 

• eliminate gas needed to offset solar pressure 
torques 

• maintain wheels at lower speed Improving their 
life and reliability 
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• eliminate almost entirely the need for wheel 
momentum dumping 

• allow possibility of correction for SA deployment 
misalignments 

• enhance multimission applicability 

A possible mechanization^ for achieving this vehicle 
turn capability is shown in Figures 1-3 and 1-4. 





Figure 1-3. CM/ CP Trim Concept 


t^Cork, M.J., et al., "Planetary Mission Requirements, Technologies and 
Design Considerations for a Solar Electric Propulsion Stage," A1AA 
Paper 79-0908, presented at the AIAA Conference on Advanced Technolo- 
gies for Futi *e Space Systems, Hampton, Virginia, May 8-11, 1979. 





f U* j * t 



Figure 1-4. Solar Array wing CM/CP trim possible mechanization 



SECTION 2 


SEPS/ICM DYNAMIC MODELING 

In order to evaluate the proposed attitude control schemes 
the SEPS/ICM vehicle must be modeled Into a form suitable for digital 
computer simulation. This section presents the assumptions made In 
this modeling and the abstraction process which takes the vehicle into 
a dynamic model consisting of a series of hinge-connected rigid bodies 
with flexible appendage*. 

The model to be developed here will correspond to vehicle 
characteristics near the time of rendezvous with comet Tempel 2, since 
that is the time when gas jets or reaction wheels, will provide the 
primary attitude control. At that time most of the mercury propellant 
has been spent and the vehicle Is significantly lighter than the 
Initially launched mass. 

This section is organized as follows: 
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SEPS/ICM DESCRIPTION AND ASSUMPTIONS 


- 2.1 

^ — ~ The SEPS/ICM Vehicle will beassumed to consist of the 

following components: 

1) The SEPS Stage, consisting of 2 articulated flexible 
lockheed/MSFC solar array wings plus a rigid SEPS bus 
(comprising the remainder of SEPS, including the inter- 
face adapter, propulsion module, etc.). 

2) The ICM spacecraft, consisting of a rigid ICM bus plus 
an articulated (2 degree-of -freedom) rigid science 
scan platform ( SP ) , as shown in Figure 2-1. 


Scan 



Fig. 2-1. SEPS/ICM Components 
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For the purposes of this analysis, we wil l lump the. SEP5> bus. 
Interface adapter, andICM bus into one single rigid body b Q which will, 
henceforth, be designated as the "vehicle bus" or, simply, "bus". Under 
these assumptions, the SEPS/ICM vehicle will be modeled as shown in 
Figure 2-2. The model consists of the following four substructures. 

1) A rigid vehicle bus (b Q ) (made up of the SEPS bus, 
interface adapter, and the ICM bus). 

2) Two independently articulated flexible Lockheed/MSFC 
solar arrays. (Flexible model is described in 2.2.) 

3) A two-degree-of-freedom articulated Science Scan 
Platform (SP) . 

2.1.1 Coordinate Systems 

» A total of four dextral coordinate systems are necessary to 

describe the SEPS/ ICM as shown in Figure 2-2: one for the bus, one for each 
solar panel, and one fair the scan platform. Throughout this study we 
will use the following definitions for these coordinate systems. 

Bus Coordinate System (X p Y g Z D ) 

The bus coordinate system consists of three mutually perpen- 
dicular reference axes embedded in the vehicle bus: bus pitch (X g ), 
bus yaw (Y g ), and bus roll (Z g ), and is defined as follows: Let the 
"interface plane" be defined as the imaginary plane formed by the mating 
surfaces of the SEPS interface adapter with the base of the ICM space- 
craft, then 

The bus roll axis (+Z g ) is the normal to the interface 
plane passing through the geometric center of the 8 
thruster array, with +Z g in the direction of the 
exhausting propellant. The intersection of the Z g axis 
with the interface plane defines the origin of the bus 
coordinate system. 

The bus pitch axis (+X g ) lies in the interface plane 
and is parallel to the thermal radiator surfaces, with 
-X g in the hemisphere containing the scan platform. 






< 

V) 
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Figure 2 - 2 . SEPS/ICM Model and Coordinates 



The bus yaw axis (+Y g ) lies in the interface plane and 
completes the dextral coordinate system. 

Solar Array Coordinate Systems (X^Y^Z^) and 

The solar array panel located in the +Yg hemisphere 
(Figure 2-2) will be here designated as "SA^.” Similarly, the -Yg solar 
panel will be designated "SA^." Since each panel will be articulated 
independently in this study, each of the two solar panels will have a 
separate coordinate system embedded in it - (X A1 Y fll ) and 
(X A 2 Y Z^)» respectively. Both panels will be articulated about a 
cornnon hinge line parallel to the Y g direction but displaced 0.2 m 
below the interface plane in the *Z g direction. The point of inter- 
section of the solar array common hinge line and the *Zg axis will be 
termed the solar array (SA) hinge point, and its location is 
(0, 0, +0.2)m in bus coordinates. 

The origin of coordinates for (X ftl Y ftl Z A1 ) and . 

(X ft 2 Y a? Z A2 ) is common and will be taken at the "SA hinge point" 
defined above, with 

♦Y ft1 and +Y A2 both coincident and parallel to +Yg 

♦Z^ and * 1^2 Are defined by the normals to the — 
sensitive (solar cell) sides of the undeformed solar 
cell blankets 

♦X A j and +X A2 are defined to complete the two dextral 
systems. 

Scan Platform Coordinate Systems (CL CN) and ( X 5p Y sp*SP ^ 

Two coordinate systems will be useful in the description 
of the scan platform, which is a two degree of freedom platform 
articulated by means of "clock" (CL) and "cone" (CN) actuators, as 
shown in Figure 2-2. The platform coordinates are defined in terms 
of the look direction L of an imaging instrument mounted on it, as 
follows 

(0° CL, +90° CN) corresponds to L looking in the 
+X B direction 

(>90° CL, +90° CN) corresponds to L looking in the 
+Yg direction. 
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The ( X^p Yjp Z<jp) d extral coordinate system has its origin at. the 
intersection of the clock and cone hinge lines and is such that it is 
parallel to (Xg Yg Zg) when the scan platform is at (0 Q CL 0 a CN)., 

2.1.2 Assumed Nominal Configuration 

Since the solar panels and scan platform orientations 
relative to the bus change throughout the mission, the SEPS/ICM vehicle 
exhibits a variable configuration throughout the mission. For the 
purposes of this study, we have selected a nominal configuration which 
is defined as shown in Figure 2-3. 



The selected nominal . configuration has the SEPS thrusting normal to the 
sunline with +Xg, +7^, and + 1 ^ looking at the sun and the scan 
platform at (0° CL, 0° CN). The orientation of the rig axis in inertial 
space is arbitrary. 

2.1.3 Solar Pressure Disturbance Torques at Tempel 2 

Solar pressure disturbance torques can be calculated 

from 

T dist a 1°' 6 0+e) , N-* . 
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where 

€ 

3 surface reflectivity 


A 

3 area , m 2 


t 

3 offset between center of pressure and center of 



mass, m 


r 

3 distance from the sun, AU. 


Using two 32 x 4 in panels yields A * 256 m 2 . Assuming the center of 
pressure to be at the geometric center of the two solar panels* i.e. 
at the SA "hinge point" located at (0, 0, +0.2)m In bus coordinates, 
yields i 3 0.3 m since the vehicle center of mass is at (0, 0, +0.5)m 
In bus coordinates. Using a reflectivity of e 3 0.2 and an average 
sun distance r of 1.43 AU at Temple 2, yields 


3 +2.09 10" N-m about Yg. 


Solar pressure torques about the other two axes are zero. The effect 
of this solar pressure disturbance torque Is the following: 


Thrust Vector Control Mode . Assuming n operating thrusters 
with a thrust of 0.125 N each, all simultaneously gimballed by an angle 
a and located In a plane 2.71 m away from the vehicle CH, results in a 
steady state gimbal angle o required to offset the solar torques which 
can be obtained from 

0.125 n i sin a 3 2.09 10' 4 


I.e. 


o = arc sin 


,2.09 10" , 
rT V 


which gives, with i = 2.71 m. 


0.0044 degrees for n 3 8 operating thrusters 

, 

0.0176 degrees for n = 2 operating thrusters 


for the gimbal angles required to offset solar pressure torques. 

Reaction Wheel Control Mode . If we assune the use of a 
NASA Standard Reaction Wheel (SRW) having a momentian storage capacity 
of 20 N-m-s @ 2200 rpm, the Y-axis SRW will speed up continuously 
from 0 rpm to saturation (2200 rpm) in a time T given by 

T „. -?0 . . = 

2.09 10- 


95693 seconds 



that is, the Y-axis SRW will saturate in 26.58 hour.; due to solar 
pressure, at which time it must be unloaded. Note that the only way to 
reduce the frequency of unloadings is to reduce the center of mass to 
center of pressure offset by bringing the two closer together. 

Reaction Control System Mode . In a 24 hr. period the gas 
jets must overcome a one sided solar pressure torque producing. a total 
angular impulse of 

H = 2.09 10" 4 x 24 x 60 x 60 = 18.0576 N-m-s 

. If we assume that the gas jets - center of mass lever arm 1st 8 2.71 m, 
the mass consumed to offset the solar pressure can be computed from 

m = £ = 9.06 gr/day 

sp 

where we have used I ■ * 75 sec (N g gas). Note that the only way to 
reduce this consumption is to bring the centers of mass and pressure 
closer together. 

2.1.4 Mass Properties 


The total vehicle mass properties assumed in this study 


are as follpws: 

. r ~ r — 



Total mass 

2200.3 kg 



CM location in 




( X b Y B Z b) coord * nates 

(0, 0, +0.5)m 



Self inertia natrix in 

1 179000. 

0. 

-60. 


nominil configuration 

0. 3500. 

0. 

kg-m 2 


-60. 

0. 

178500. 



The mass breakdown for each component is as shown in 


Table 2-1. 


The solar arrays a.*e assumed to be connected to the bus 
through a 1 dof line hinge located at (0, 0, +0.2) m in bus coordinates. 
The scan platform is connected to the bus through a 2 dof hinge located 
at (-1.3691, -0.37, -0.7893)m.in bus coordinates; this implies that the 
center of mass of the scan platform is located coincident with the . 
intersection point of its two hinges. 
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TABLE 2-1. SEPS/ICM Mass Properties-(Nominal-Configuration) — 


Mass, kg 

CM Location^ 
(X,Y,Z) , m 

Self- 

(?) 

inertia Matrix 
kg-m 2 

BUS 6 0 

+0.0703 

1824.304 

+50.133 

+114.634 

1756.474 

+0.0190 

+50.133 

2635.642 

+46.2,49 


+0.6266 

+114.634 

+46.249 

1769.269 

SCAN 

PLATFORM b 3 

-1.3691 

13.300 

-2.100 

+0.200 

90.190 

-0.3700 

-2.100 

7.500 

+1.000 


-0.7893 

+0.200 

+1.000 

16.500 

♦ Y SOLAR 
ARRAY fa. 

0.0000 

22084.000 

-158.020 

0.0 

176.818 

+19.3764 

-158.020 

234.610 

• 0.0 


+ 0.2000 

0.0 

0.0 

21876.000 

-Y SOLAR 
ARRAY b 2 

0.0000 

22084.000 

+158.020 

0.0 

176.818 

-19.3764 

+158.020 

234.610 

0.0 


+ 0.2000 

0.0 

0.0 

21876.000 

TOTAL 

VEHICLE 

0.0000 

179000.000 

0.000 

-60.000 

2200.3 

0.G000 

0.000 

3500.000 

0.000 


+0.5000 

-60.000 

0.000 

178500.000 


(1) referenced to (X g Y g Z g ) system 

(2) reference to a coordinate system parallel to (Xg Yg Zg) but 
centered at the center of mass of each structure 
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2.2 LOCKHEED/MSFC SOLAR ARRAY FLEXIBLE MODEL 

This subsection describes the hybrid modal model used to 
represent the flexible solar arrays In the computer simulations performed 
In this study. The hybrid-mode approach [A] is a cross between 
the free-free approach and the cantilever mode approach In that the rigid 
base of each solar array wing is considered free in translation but fixed 
In rotation. For the translational equations of each array, the free 
rigid base Includes the mass of the remainder of the vehicle (as if It 
were rigid. Including the other solar array wing). 

A hybrid modal model of each array was generated [5] 
using the finite element discrete model shown in Figure 2-4. This 
figure shows the +Y solar panel (SA*) In the cruise configuration, 
attached to the vehicle bus at node 1 through a rigid link RL. The mass 
of the rest of the vehicle was lumped at that node, which In turn was 
connected to the center of mass of the panel by a massless rigid link. 
Node 1 was allowed to have translational motion only, the rotational 
degrees of freedom being set to zero. Both pre-stressed membrane and 
line elements have been used for the discretization. The associated 
geometrical stiffness matrices Kq for the solar panel membrane elements 
were derived by applying pre-stressing along the length of the panel 
(Y direction), while restraining the two longitudinal edges 
(Z - + 2 meters) from any possible lateral movement. Such a measure 
was necessary to prevent any numerical Instability caused by the 
Poisson's ratio effect, resulting in small or negative diagonal elements 
In the geometrical stiffness matrix terms pertaining to the lateral 
(Z direction) motion of the panel. The program EISI/SPAR utilized 
for the analysis was then found to be effective for the solution of the 
problem, although the frequency values obtained are slightly higher 
than expected. 


[4] 


[5] 


Likins, P.W., Dynamics and Control of Flexible Space Vehicles, 
Technical Report 32-1329, Jet Propulsion Laboratory, Pasadena, Calif.,' 
Jan. 15, 1970 

Gupta, K.K., SEPS Solar Array Panel Natural Frequency Analysis, 
Interoffice Memorandum 3545:79:167, Jet Propulsion Laboratory, 
Pasadena, Calif., May 30, 1979 (JPL Internal document). 
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The mass breakdown assumed for solar array wing Is given 
In Table 2-2. This model has blanket dimensions of 4.0 x 31.6 meters 
with a total w.1ng structure length of 35.0 m from node 1 to node 92, 
and a center of mass for the wing located at 

(+0.1242, +16.5474, 0.0) meters 


Table 2-2. Mass Sunmary - Solar Array 


Item 

kg 

1b 

Flexible boom 

18.05 

39.80 

Canister 

18.73 

41.29 

Container 

11.62 

25.62 

Position boom 

3.18 

7.01 

Blanket* 

115.83 ' 

255.36 

Box cover 

8.57 

18.89 

Tip fitting (Mast) 

0.84 

1.85 

TOTAL 

176.82 

389.82 


*Array Guidance Mechanism and Tension Mechanism were 
assumed to be uniformly distributed along blanket. 

with respect to node 1. The rigid link RL is used to connect node 1 
to the SA hinge point hp (0, 0, +0.2) so that the total length from 
hp to the SA\ wing tip Increases to 37.829 m and the center of mass 
of the SAj structure Is positioned at 

(0., +19.3764, 0.) m with respect to hp. 

This places the CM of the SA's on the hinge line. The resulting 
hybrid modal data for SAj is given in Table 2-3 which lists the 
frequency and the rigid elastic coupling coefficients. The corres- 
ponding mode shapes (eigenvectors) are shown in Figure 2-5. The data 
for SA 2 was obtained by rotating the rigid elastic coupling coefficients 
and the mode shapes by a 180° rotation about a line parallel to Xg 
passing through hp. This simply changes the sign of the Y and Z 
translations and rotations. 


-rX-- 


\ V 
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TABLE 2-3. Lockheed/MSFC SA Hybrid Model Data (2200.3 ko vehicle) 

Mode • . Frequency, . Rigid-elastic Coupling Coefficients 

Nn *yP e u, li ,kq*S u u 6 .kcrtm A 
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Figure 2-5. Solar Array Eigenvector? (Mode-shapes) - Exploded View 
(Motion of blanket and boom not to same scale) 
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Figure 2-5 (contd) 
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Figure 2-5 (contd) 
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2.3 SOLAR ARRAY AND SCAN PLATFORM ARTICULATION AND CONTROL 

Each wing of the solar array Is assumed to be Independently 
articulated about Y ftl and by means of Independent actuators. The * 

details on the various actuator characteristics and control schemes 
for the SA considered In this study can be found in subsections 3.2.4 , 4.2, 

4.3, and 5.3.3. ' v 

The scan platform model and control method, on the other hand, 
have been constant throughout this stud and will be described here. 

Each of the two SP axes (CL and CN) will be controlled 
by the simple scheme shown in Figure 2-6. 



Figure 2-6. Scan Platform Block Diagram 

The torque on the scan platform Is given by 

T * - lc(e - e c ) - B e 

where IT a K )^. 

Such a block diagram can be thought to arise from an 
actuator drive mechanism having control gain K, motor torque constant 
and a viscous damper at the output shaft having a damping constant 
of B N-m per rad/s. 


I 

a 
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The pertinent constants are: 


CL CN 

I 16.5 7.5 kg-m 2 * 

Y 1500 1500 N-m/rad 

B 220.25 148.50 N-m-s/rad 

The high value of K is typical of scan platforms flown on earlier 
missions (Viking, Voyager, etc. ) and is required to meet pointing 
accuracy requirements in the presence of frictional forces, although 
such forces have been neglected in this study. The choices for B 
provide 0.7 damping. With these values, the natural frequencies for 
the scan platform are 


'CL S 

9.53 rad/s 

<CL " °' 7 

CN = 

14.14' rad/s 

{ CN *' 0 7 


2.4 DESCRIPTION OF TEST MANEUVERS ■ 

In this study, we will use the following four (typical) 
maneuvers for evaluation of the control schemes: 

i) X TURN : Starting with the vehicle in the nominal configu- 

ration, of Figure 2-3, the complete vehicle is commanded 
to turn about the sun line (pitch). The command is a 
position ramp with a slope of 0.005°/s. The solar arrays 
are conmanded to hold their position relative to the bus. 

ii) BUS YAW TURN . The vehicle is commanded to turn about Yg, 
while the solar. arrays are commanded to remain pointed 
at the sun, as the bus moves. The command to the bus 
is a position ramp with a slope of 0.25°/s. 

1 i i ) MULTIAXIS ACQUISITION . The attitude control system is 
enabled at a tine when the vehicle has the following 
initial position errors and rates:. 



Axis 

Initial Position 


Error 

Initial Rate 

pitch 

+1.0° 

♦0.005VS 

yaw 

+2.0° 

♦0.250°'/ s 

roll 

o 

o 

0 

-0.005°/s 


iv) SCAN SLEWING . The dynamic disturbance .caused by scan 
platform slewing depends on many factors such as the geometry and mass 
properties of the platform and vehicle, the size and direction of the 
slew, the slew rate, and the initial' clock and cone orientation of the 
platform before slewing. In particular, these factors determine the 
angular rates and momentum imparted by the slew and which must be 
controlled and absorbed by the attitude control system - whether it is 
based on reaction wheels, thrust vector gimballing, or on a gas jet 
system. 

Since it is impossible to study all possible sets of scan 
slewing conditions, we have assumed a slew test sequence. This is a 9°x2° 
box or raster type slew, at l“/s. It consists of four 9°-clock slews 
interspersed with four 2°-cone slews and 2 second wait periods in 
between each slew; this 8 slew sequence takes 58 seconds to complete. 

For a graphical description of the sequence, the reader is referred 
to figure 3-14 v), in the next section. 

2.5 COMPUTER SIMULATION PROGRAMS 

Three separate flexible dynamics simulation programs were 
developed, named "TVC," "RCS," and “RWH." The programs are tailored to 
simulate each of the three control systems considered in this study: 
thrust vector control (TVC), gas jets (GAS), and reaction wheels (PWH). 
Each program incorporates the corresponding control system and hardware 
models described in Sections 3, 4, and 5. 

The three simulation programs have several similarities 
such as the use of a powerful subroutine called FBDYFL. This hybrid- 
mode subroutine written by Fleischer and Likins [6] solves the equations 
of rotational motion for a system of hinge-connected rigid bodies with 

^Fleischer, G.E., and Likins, P.W., Attitude Oynamics Simulation Sub- 
routines for Systems of Hirge-Connected Rigid Bodies with Nonrigid 
Appendages, Technical Report 32-1598, JPL, Pasadena, CA, Aug. 15, 1975. 
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flexible appendages. For inputs, the subroutine uses the mass properties 
and vehicle configuration described, in 2.1, as well as the hybrid SA 
model of 2.2. Each SA panel was modeled by its first IS hybrid modes 
and corresponding 'frequencies, yielding a 30-mode flexible model for the 
complete SEPS/ICM vehicle. 

The simulation programs were alt written in the "Continuous 
System Simulation Language" (CSSL III). As the name implies that 
language is used to digitally simulate continuous systems. A complete 
listing of the three programs appear i'r. appendices A (TVC), B (GAS), 
and C ( RWH ) . The four maneuvers simulated by each program were 
described in 2.4, and the results of each are given in sections 3 (TVC), 

4 (GAS), and 5 (RWH). 


SECTION 3 


THRUST VECTOR CONTROL 

The thrust vector control "TVC" system is a three axis 
active controller to be used for vehicle control during powered 
flight. The TVC system utilizes the ion engines to generate vehicle 
control torques. This is accomplished by gimbaling each engine 
on a two degree of freedom gimbal. Because of the low thrust provided 
by the ion ngines the control torques are in general small. This 
coupled with the fact that the vehicle has large rotational inertias 
results in slow response and lengthy maneuver times. With suitable 
control laws though, these small control torques provide very smooth 
vehicle maneuvering and minimal excitation of the vehicle's flexible 
solar panels. 

This section will discuss the TVC attitude control system 
for the SEPS/ICM vehicle. The performance of this system Is also 
detailed and evaluated. The section is organized as follows: 
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3.1 


DESCRIPTION OF THE TVC HARDWARE 


The SEPS/ICM vehicle will use a thrust vector control 
system for attitude control during powered flight. 

The proposed TVC hardware configuration is shown in 
Figure 3-1. Thrust is provided by eight 30 centimeter mercury ion 
engines. In this study we have assumed that all eight engines are on 
and at full thrust, yielding a sum total thrust force F c 3 1 Newton. 
Vehicle control is accomplished by gimbaling each engine on its own 
2 degree of freedom set of gimbals, as shown in Figure 3-2. Each set 
.of gimbals is assumed to have two actuators and two sensors. 

Since each engine Is gimbaled individually. It Is possible 
to control each engine Independently of the others by Issuing a total of 
8 gimbal angle command pairs. This approach leads to complex control 
logic. In an effort to reduce this complexity, it is desirable to 
design a system which Issues the same commands to as many thrusters 
as possible. In this study we have made use of a scheme which reduces 
the number of different gimbal angle command pairs to the minimum, 
which is two pairs. This Is done as follows. In reference to 
Figure 3-1, let the four engines labeled 1, 2, 3, and 4 be known 
collectively as Engine dank #1 or, simply. Bank 1. Similarly, engines 5, 
6, 7, and 8 become Bank 2. We will assume that all eight engines, 
gimbal hardware, controllers and associated sensors and actuators are 
identical. Then, it Is possible to control the vehicle attitude by 
simply issuing two pairs of coomand gimbal angles (c* c x ,e c j ) and (« c2 ,0 c2 ) 
to the engines in Bank 1 and Bank 2, respectively. Since all four 
engines in each Bank receive the same conmand pair and since the 
controllers and hardware are Identical, all four engines in each bank 
will move simultaneously and undistinguishably from each other - as a 
bank Indeed. 


3-2 




3.2 


TVC LEAD-LAG CONTROL SYSTEM 


A very simplified schematic of the TVC control system used 
In this study is shown in Figure 3-3. This schematic. illustrates the 
basic control loop configuration which was provided to JPL by MSFC for 
analysis and evaluation.* 

The proposed controller is of the "lead -lag'' type and it 

achieves three-axis control as follows. The position error e is obtained 

by subtracting the (celestial sensor) position e $ from the comnanded 

position e c . The position error is then filtered. This filter is 

Included to provide some isolation of the loop from the higher frequency 

glmbal dynamics which might be sensed and fed back through the celestial 

sensor. The output of the filter e f Is then put through a lead-lag 

compensating network to derive the control error e £ . The lead-lag 

network ll— * can be thought to consist of a filter (t 2 s+1 ) -1 
t 2 s+l 

yielding a "doubly filtered position error signal" e ff followed by a 
"signal plus derivative" compensator 1+tjS so that e c can be thought of 
as the "position plus rate" of the twice-filtered position error e ff . 

The control error e is then multiplied by a gain K to obtain y , the 

c c 

"thrust vector (TV) command angle." The three TV command angles 
(one per axis) are then fed Into a "Glmbal Angle Commander" black box 
which generates two pairs of glmbal angle commands (° c i»8 cl ) and 
(<»c 2 • &C 2 ) for each of the two banks of engines. Each engine has two 
controller subloops (one per glmbal) and associated actuators and 
sensors to move the gimbals to the commanded angular positions, thus, as 
the engine banks move, control torques are applied to the vehicle, as 
Indicated by the "force and torque" block, which cause the desired 
attitude control. 

We shall explain next how the control torques are obtained 
from the y c 's by describing In detail the 3 black boxes labeled 
"Glmbal Angle Commander," "Glmbal Control Subloops," and "Force and 
Torques." 


*A complete block diagram is shown in Figure 3-10, page 3-1% 












3.2.1 Glmbal Angle Commander 

As we have described earlier, we made the decision to 
glmbal the 8 engines as two banks of 4, called Bank 1 and Bank 2, 
because of its attractive simplicity. As it turns out, such gimbaling 
philosophy results in even further simplification of the analysis, since 
the resultant thrust of the four engines in, say. Bank 1 can be thought 
as arising from an imaginary "Equivalent c ngine (EE) 1" having a thrust 
of 0.5 N, i.e., 4 times that of a sing 1 thruster (f = 0.125 N), and 
located at the geometric center of the Bank 1 thruster pattern. 
Similarly, "Equivalent Engine (EE) 2“ arises from Bank 2 and has 
analogous characteristics. This simpl Ification Is shown in Figure 3-4. 

Also shown in Figure 3-4 are the coordinate frames and 
glmbal angle conventions used for each bank. Note that the frames 
for EEl.and EE2 are related to each other through a 180° rotation about 



Figure 3-4. Bank 6imbaling Conventions 
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Figure 3-5. Method of Obtaining Pure Positive 
Torques about each of the 3 axes 


Figure 3-5 shows how the gimbals can be moved to produce 
pure torques about each axis. Thus, If we wish to obtain a positive 
control torque about Xwe can command a positive TV control angle y cx 
by glmbaling banks 1 and 2 by that angle, as shown. The resultant 
thrust vector forms an angle with the Zg Xg - plane of r cx «nd results 
In a torque about X of 

T x *♦ F L s,n ’ex (: * F L ’cx • ,or ’cx>- 

3- 7 * 

\ 





k 


Similarly If we command a positive r , we obtain 

T y =+ F L sin Y cy (= ♦ F I. Y cy , for small Y cy ) 
and by comnanding a positive y cz > we obtain 

T z s+ F i sin y cz (= + F L y cz . for small y cz ) 

By using Figures 3-4 and 3-5 we can, by inspection, obtain the polarity 
of the glmbal. angle commands (® C1 »B C1 ) ( a C2 *8 C2 ) necessary to produce 
each pure torque commanded by y cx » Y cy » y cz * This relationship Is 
summarized In Table 3-1. 


TABLE 3-1 

8IMBAL COMMANDS TO PRODUCE PURE POSITIVE TORQUES 


Desired Glmbal Commands Needed 


Desired 

Torque 

TV Command 
Angle 

°ci 

“C2 

B C1 

B C2 

about X 

+ v 
T CX 

0 

0 

’ T CX 

%x 

about Y 

+ Y 

Cy 

” Y cy 

^cy 

0 

0 

about Z 

♦ Y 
T C 2 


^CZ 

0 

0 


/ 

K 



f It can be observed that 8-gimbaling (the outer gimbal) 

Is used to provide pitch control, while a-gimbaling (the Inner gimbal) 
Is used for both yaw and roll control. Since, in general, we will 
need simultaneous control In all three axes, it becomes necessary to 
■mix" the y cx , Y cy , and y cz coomands. This can be done, simply, by 
superposition so that 


®ci 


- r 


cy 


♦ Y 


cz 
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/ 


+ 


C2 


= + 


v + y 

'ey ’cz 


s cl 


- Y, 


ex 


This "glmbal angle commander" logic can be equivalently described 
by the schematic shown in Figure 3-6. 




0 



Figure 3-6. Gimbal Angle Ccnanander 


/Sc2 

°k2 


3.2.2 Glmbal Controller Subloops 

The 2-gimbal command pairs just generated are then fed into 
the 16 gimbal controllers (2 per thruster). The controllers energize 
the gimbal drive motors causing the inner and outer gimbals to move to 
the commanded positions. Each gimbal axis will be assumed to be 
driven by a controller of the type shown in Figure 3-7. This scheme 
applies a torque to the gimbal T * Ke which is proportional to the 
gimbal position error. The rate feedback path provides damping. In ^ 

reality the system need not use actual glmbal rate feedback, the same 
damping can be achieved if a suitable viscous damper Is Included. The 
limits of ♦ 30 degrees for position and ♦ TO degrees/sec for rate are < 

included to model expected gimbal hardware limitations. 

* 
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Figure 3-7. Gimbal Control Subloop (1 of 16) 


In the nonrsaturating region the gimbal control loop behaves 
as a simple second order system with natural frequency and damping 

o 1C I ■ >* 

(j ^ 3 ■ • — and £ 3 » K u , 

n * 7 rp n 

The frequency u» n should be chosen high enough so that the gimbals 
track the commands <* ci or ^properly but, at the same time, not so high 
that gimbal dynamic transients cause unwanted excitation of the 
flexible structure. Ue have chosen 

w n *0.2 rad/s (0.032 Hz) . 

which Is approximately one octave below the lowest flexible mode of 

the structure. This results in K * Q.Z 2 I„. 

e 

Each thruster will be assumed to have an Inertia 
I e *0.475 kg-m 2 about its gimbaling axes (the gimbals themselves are 
assumed to have negligible mass). This results in a control gain 
of K * 0.019 N-m/rad. Selecting e * 0.707 for damping Implies a rate 
to position gain of * 7.07 s. This completes the description of the 
the gimbal controllers. 
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3.2.3 


Generation of Forces and Torques 


The total force acting on the bus due to the gimbaled 
thrusters can be computed, simply, as the vector sum of the forces 
produced by each thruster. 

The total torque produced by the thrusters on the bus, on 
the other hand. Is somewhat more difficult to obtain. It must be 
recognized that the gimbaling of the thrusters produces torques on the 
bus which arise from two Inherently different sources: 

PRIMARY CONTROL TORQUES. In the case of pitch and 
yaw these are the torques which cause vehicle 
rotations when the resultant thrust vector does 
not pass through the vehicle center of mass, C.M. 

For roll these torques are produced by unbalanced 
force components In the X Y plane even If the net 
thrust passes through the CM. Note that these 
torques are a function of the Instantaneous 
pointing (position) of the engine thrust axes. 

DYNAMIC REACTION TORQUES. These torques reflect 
the transient torques that are experienced by 
the bus as we rotate the gimbals. They are a 
consequence of the engine Inertial reaction to 
rotation, and depend on the glmbal accelerations 
and velocities. 

In order to compute the resultant forces and torques we must first 
obtain the set of coordinate transformations between bus coordinates 
and each set of thruster-fixed coordinates. 

Thruster Coordinate Transformations 

The glmbal coordinates and configuration shown In Figure 3-8 
will be assumed to be comnon for the 4 engines in Bank 1. 



9S- 


Figure 3-8. Bank 1 Glmbal Coordinates 



As shown in the figure, let { b 1 } . {b 2 }, (b 3 > denote three reference 
frames defined by triads of unit vectors x^, y. , z., 1=1 ,2,3 embedded 
on the bus, the outer gimbal and the thruster (i-e. the inner gimbal), 
respectively . The gimbal angles of all the engines in Bank 1 is (a lt 3i), 
so that frames {b 1 } and {b 3 } are related to each other through a 
rotation about X! followed by a o x rotation about y 2 so that the 
transformation! T relating both frames can be expressed as 

cosa 0 -sina 1 0 0 

(b 3 > = T {b 1 } - 010 0 cosB sine 

sina 0 cosa 0 -sinB. cosB 
b J L 

and 

. COSa 0 Sina 

{ b 1 } » b 3 } * sina cosB cosB -cosa sir.8 

-siha COSB -cosa sinB cosa COSB 




Force Computations 


Let f denote the magnitude of the force produced by one 
single thruster. Then, this force, expressed in vector form, is given by 
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for each of the thrusters in Bank 1. Letting Fi denote the resultant 
force vector for all 4 engines of Bank 1 yields 


F = (b 3 ) 1 


0 

0 

-4f 


in frame (b 3 ). 


The components of this equivalent engine 1 thrust, expressed in bus 
frame { b 1 > are then 


F E) . u«>* 

v 

F Eiy 

- 

* 

0 

0 

= { L» 3 > t 4f 

0 

■» 

-sina! 

+cosa i sinBj, 


— 1 
M 

LU 

U. 

1 


-4f 

0. 


-cosa i cosBi 


For thruster Bank 2 the whole thruster assemblies are rotated 180° 
about Zg, as shown in Figure 3-4, so that the resultant thrust vector 
F e2 from Bank 2 can be obtained simply by taking F^ evaluated at 
(®2 > 82 ) and changing the sign of the x and y components: 


F E2- 


{bi} 1 


' F E2x‘ 

* { b 1 } l 4f Q 

P - 

+sina 2 

F E2y 

-cos* 2 sins 2 

* 

fsl 

C\J 

LU 

U. 

J 


-COSa 2 COSB 2 

» 


so that the net force acting on the bus due to the thrusters is 

r 


FB - F £ ,. F £2 = (b>) T 4f 0 ; 


-sinaj + sino 2 

cosoj sinBj - COS02 sfn0 2 

-(cOSc*! COSBi + COSo 2 COS8 2 ) 


Primary Control Torques tB c 

The primary control torque about the bus center of mass 
(CMg) can be obtained as follows: Let and R^ denote the vectors 

from CMg to equivalent engines 1 and 2, respectively. Then the 
primary control torque. is given by 


TB_ 


R E1 * F E1 * *E2 


E2 
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Dynamic Reaction Torques 


Here we shall compute those torques experienced by the 
bus as we rotate the gimbals. They are caused by the inertial 
resistance to rotation of the engines. 


We shall assume that each thruster is seated o.i the gimbal 
in such a manner that its center of mass is, coincident with the cross- 
point of the gimbal axes, so that gimbal ing of the thrusters does not 
produce any translation of their centers of mass. It will also be 
assumed that the gimbals themselves have negligible mass so that the 
inertia of each engine/gimbals assembly is equal to the self-inertia 
matrix of the engine itself, which will be assumed to be 


0 

I 


0 

0 

L 


in frame ( b 3 ) , 


That is, the inertia for the center of mass of the thruster is I g 
about axes x 3 , y 3 , and z 3 . This assumption results in the interesting 
fact that the self-inertia of the thruster with respect to any other 
frame has the same diagonal ,(I , I 


e » * e » 


I e ) form. 


Let denote the angular velocity of, say, thruster #1 
with respect to inertial space and Tet •'/ denote the inertia dyadic 
of the thruster. Then, the angular momentum of the thruster about its 
CM is given by 

H = ./•« 

The external moment on the thruster T can then be computed from 

T = rf 

Let (b°l denote an inertial frame of reference. Let tb 1 }, tb-’l, and 
(b J l denote the frames defined earlier. Let denote the angular, 
velocity of frame i relative to frame j, then the inertial angular 
velocity of the thruster u = -» 30 can be computed by the chain rule 


= «. ♦ ^10 t 

where the superscripts refer to the frames lb 1 ), i = 0,1, 2, 3. The 
term ^ ,0 represents the inertial angular velocity of the bus. Since . 
the bus moves extremely slow when compared to the gimbal velocities 



and accelerations, we will make the simplifying assumption that the 
velocity w 10 can be neglected, that is, the torques produced on the 
gimbals as a consequence of the bus rotations in inertial space can be 
neglected. Thus, we can assume that lb 1 ) itself is fixed in inertial 
space end, therefore, 

y 3 *31 a + yjl 

By inspection of Figure 3-8, we can write 

(j) 32 s i,y 3 » i 3 y 2 
«a 2 * 3 tS j X2 3 Bi X| 


so that 


u)3l= tSjXj + «iy2 • 


Making use of the fact that 


y 2 =,cost»iyi + sinBjZ! , 


we can write 

w * w M 3 SjXj ♦ ajcose yi * sine z i 

letting u« | 1 denote the components of u> 31 expressed in frame 1 we can 
write 

apd w = y 31 * ( b 1 

The angular momentum can be expressed as 

H 3 ./ • w 3 IbMt i {b 1 1 • (bM 1 u»f 1 * lb 1 ) 1 I 


«»i! s 


U| COStf 

a, costs 


since I is the thruster inertia matrix in all frames. The external torque 
is given by _ _ . 


T 3 ~ (/Aw) = Ib 1 )^ I wj 1 



•*1 

'* • • 
a jcostSi - o l 0 l sine, 

_a l SintS| ♦ ajBi costs 


The reaction torque exerted on the bus by the 4 engines of Bank 1 is 
then, simply 


TB 


gi 


4T = {b 1 F 41 


"®i 

-aj COSBi + ajBi SinBi 
~(Xj sihBi - 3 j 6 j cosBi 


By an analogous procedure, the reaction torque exerted on the bus by 
the four engines of Bank 2 is given by 

a 2 COSB2 - a 2 B2 SinB 2 

|-a 2 sinB 2 * cosB 2 


TB 92 • (b>)* 4I e 


and the total dynamic reaction torques on the bus are 


TB = TB , ♦ TB , . ‘ 

9 9» g2 

The total external torque acting on the bus is 

TB * TB t ♦ TB 9 * T(ji st 

i.e.,the sum of the primary control torques, the dynamic reaction 
torques, and the external disturbances. 

This completes thedescription of the TVC loop. We shall 
next describe the Scan Platform and Solar Array Controllers. 


3.2.4 Solar Array and Scan Platform Controllers 

Each solar array wing has its own separate control loop 
to articulate it relative to the bus. Since each solar array wing 
controller is an Independent system, the wings can be made to 
rotate separately or together. They can also be made to follow the 
bus as it rotates, or not to follow yaw bus rotations 0 2 by commanding 
-e 2 into the solar array wing controllers (this is what has to be done, 
for example, when we wish to reorient the bus in yaw without moving 
the solar panels from the sun). 

The controller assumes that suitable actuators are provided 
between the bus and each wing of the solar array to control their relative 
positions. 
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The actuators Mill be controlled by another lead-lag type 
controller as shown in Figure 3-9. 



* 


Figure 3-9. Solar Array Lead-Lag Controller (1 of 2) 


The position of the solar array relative to the bus (as measured by 
a potentiometer or encoder) is subtracted from the command to obtain 
the position error e which is then put through a filter and a lead-lag 
compensation network to obtain the control error e c . This is similar 
to the TVC loop controller described earlier. The control error Is then 
multiplied by a gain K sa and by means of the motor torque constant 
K^ a we apply a torque to the solar array wing 


where 


T SA ■ *“ »c • 


K sa • K 54 


The scan platform controller used has been described previously. For 
details the reader is referred to subsection 2.3. 


3.2.5 System Block Diagram and Parameters 

Assembling the controllers described in the preceding pages 
together with the ICM scan platform described in 2.3, results in the 
TVC control system block diagram shown in Figure 3-10. The system 
parameters are listed in Table 3-2. 
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Figure 3-10. TVC "Lead-Lag" System Block Diagram (1 axis shown) 









TABLE 3-2. TYC "LEAD-LAG* SYSTEM PARAICTERS 
(um for <11 mm unloss otherwise Indie* ted) 


PARA7CTER 

VALUE 

UNITS 

DESCRIPTION 

I . LEAD- LAS CONTROLLER 


1 

t. 

1.0 

s 

Celestial sensor tine 

s 



constant 

V 

1S.0 

s 

Filter tine constant 


200.0 

s 

Leed-lag 

*2 

10.0 

s 

Lead-lag 

l 

3.S2 

dlnensionless 

X-axis posItion-to-TY 
angle gain 


0.11 

dlnenslontess 

Y-axis posItlon-to-TV 
angle gam 


15.70 

dlnensionless 

Z-axI s posItlon-to-TV 
angle gain 

II. GIMBAL CONTROL SUBLOOPS 



K 

0.019 

N-e/rad 

Pos 1 tlon-to- torque gain 

«rp 

7.07 

s 

Rate-to -position gain 


0.475 

Ko-n 2 

Single thruster Inertia 

' Llltj 

♦10.0 

deg/s 

Gintel rate Unit 

LI "2 

♦30.0 

099 

Global position Unit 

III. SOLAR ARRAY ANO SCAN PLATFORM CONTROL 



0.02557 

H-m/rtd 

Posltlon-to- torque gain 

x t 

15.0 

s 

Filter tine constant 

T 1 

200.0 

s 

Lead-lag 

t 2 

20.0 

s 

Lead-lag 

* ,p «£ p 

1S0Q.0 

N-n/rad 

Scan posltlon-to-torque 
gain 

K ,p 

0.14683 

s 

Scan rate-to-posltlon galn(clock) 

rp 

0.09901 

% 

Scan rate-to-posltlon galr.(coee) 


3.3 TVC PERFORMANCE 

The performance of the TVC lead-lag system just described 
was evaluated through computer simulation of the four test maneuvers 
described In 2.4. Full details of the flexible dynamics models, 
simulation programs and assumptions can be found in Section 2. 


The results of the simulations are shown in Figures 3-11 
through 3-14 and discussed below in some detail. The figures them- 
selves are also captioned with pertinent explanations and conmentary. 
The nomenclature used In these figures is listed in Table 3-3. 


The simulation program was written In such a way that It 
assumes the vehicle to be initially moving in a steady state (constant 
linear or angular velocities). As the simulation starts and the ion 
engines provide thrust, the vehicle experiences a change to a state 
of constant linear acceleration, which brings about a constant 
deformation of the solar panels and possible shifts of the center of 
mass (and accompanying rotations). This dynamic transient is due to 
what we might call the "engine turn-on transient.” It is, generally, 
very small but was found to be large enough to mask the vehicle 
response for some of the planned maneuvers. In an effort to reduce 
this undesirable tum-on-trahslent effect, all TVC simulation programs 
were modified to ramp up the thrust in a very smooth fashion during 
the first 20 seconds of the simulation. The thrust is assumed to 
increase smoothly from 0 to 1 Newton as follows: 


F * 



for 0 < t < 20 s 


, for t < 20s . 


This ramping up Is so smooth compared to the solar panel resonances 
that It produces practically no significant oscillatory transients 
In the solar panels. 


The 1 degree pitch turn maneuver simulation (Figure 3-11) 
demonstrated the TVC system's capability to make smooth turns with 
minimal excitation of the vehicle's flexible dynamics. The vehicle 
achieved the turn rate of .005°/ s in approximately 150 seconds and 
completed the maneuver in 800 seconds. There Is a small amount of cross 
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axis motion In yaw and roll for this maneuver, as is for most maneuvers. 
This is a result of the vehicle dynamics and not a controller generated 
problem. The maximum gimbal angle was 3.72° in 8 while the maximum 
gimbal rate was ,12°/s also in g. The solar array panels experienced a 
small rotation due to the cross coupled motion of the bus in yaw but 
the controller was able to correct this condition. The generalized 
coordinates plots of SA modal deformation indicate that their flexible 
body motion was smooth and non-vibratory. All coordinates settled 
to some non-zero steady-state value as a result of thrusting in the 
-Z direction throughout the length of the maneuver. Note that most of 
the modal deformations occurred within the first 20 seconds of the 
simulation. 


The bus yaw turn simulation (Figure 3-12) shows again that 
the TVC system can make smooth turns. The bus position has one over- 
shoot and settles to the final 30° command position in about 800 
seconds. Pitch and roll exhibit very small cross-coupled motion for 
this maneuver. The bus achieved the commanded turn rate of .25°/s 
In 120 seconds then continued to Increase until It reached a peak 
value of .34°/$. The yaw control torque was amooth and the gimbal 
subloop tracked the desired torque very well. The roll torque had a 
small transient due to thruster turn- n effects during the first 20 
seconds. The maximum gimbal angle during the maneuver was about 6.3° 

In a and the maximum gimbal rate was ,22°/s also In a. During this 
maneuver the SA wings were commanded to remain in their initial position 
(normal to the sun line) while the bus was turning. They performed 
as commanded with negligible error. The modal deformations in the SA 
panels were primarily due to the linear acceleration of the vehicle in 
the -Z direction. 

The performance of the acquisition maneuver is highly 
dependent upon the initial conditions specified, but for the set used 
in this simulation the performance was very good (Figure 3-13) The 
set of initial conditions is listed below: 


Axis 

Initial Position 

Initial Rate 

pitch 

♦1 .0° 

+0.005°/s 

yaw 

+2.0° 

+0.250°/s 

roll 

o 

o 

0 

-0.005°/s 
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The vehicle was stabilized to a zero final condition in approximately 
700 seconds. The maximum position error, 13.3°, occurred in yaw (the 
axis of least Inertia) at 88 seconds. Since pitch and roll axes have 
higher inertias their positions did not drift as much as yaw. The 
maximum pitch and roll angles were 1.2° and 0.4° respectively. The 
engine gimbal subloops tracked the desired torques very well and 
exercised maximum gimbal angles of 21° in a and 7° in B. As the vehicle 
moved predominantly in yaw, the SA wings also experienced large angular 
rotations since the SA actuators have very low control authority 
(low positlon-to-torque gain). The modal deformations reflect most 
the response in yaw while constant deformations due to engine thrusting 
are also present. 

For the scan slewing simulation (Figure 3-14), the TVC 
system was comnanded to maintain the vehicle In its initial zero state 
while the science platform was scanning. The box slew sequence Intro- 
duces high frequency disturbances Into the vehicle, with yaw being the 
axis most affected. The engine gimbals reached a maximum angle of 
.11° in a to compensate for the .014° movement in yaw. The most 
noticeable characteristic of the controller during this maneuver Is 
Its low pass filtering of the Input disturbances. The roll torque 
(plot n) shows that the controller tracks the low frequency components 
of the input well but filters out the high frequency transients. For 
this maneuver the SA actuators were comnanded to hold the wings in 
place. The solqr panel vibration observed In the generalized modal 
coordinates is by in large due to the Impulsive type disturbances 
transmitted to the vehicle during science platform slewing. 

In summary, the TVC control system provided very smooth 
proportional control of the vehicle for the four maneuvers studied, with 
negligible excitation of the vehicle flexibility. The required gimbal 
ingles were acceptable. 
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ALPHl ,2,3* 

AnguUr position of the bus (b Q ), rad 

wwn 

System angular moment uo magnitude, N-«-s 

C0ERR1.2.3 

Control error e c , rad 

COACL, COACH 

Position cooaands Into Scan Platfone, rad 

C0ASA1 .2 

Position coownds Into SA wings, rad 

ERR0R1 ,2.3 

Sensed position error e, rad 

ETAU to ETA18 

Panel deformation generalized coordinates 
(pane) 1, first 8 nodes) 

SIAAI.2 A GIAB1.2 

Global angles a,. a,. a., s 7 , respectively, 
rad 14,4 

GIMA1Q.2D & GM10.20 

Global angle rates a.. a,. i 7 , 

respectively, rad/s 

GM1. 2.3.4 

Hinge rotation angles for Solar Array 
Paitels 1,2, and Scan Clock, and Cone, 
respectively, rad 

RATE1 ,2,3 

Bus angular rates, rad/s 

T0ESR1 ,2,3 

Desired control torques. N-n 

TH1 .2,3,4 

Hinge torques for Solar Array Panels 1,2, 
and Scan Clock, and Cone, respectively, rad 

TJWCOA 

* c position input (turn comnd), rad 

TSCI .2.3 

External torques about vehicle CM, N-e 

# 

1,2,3, denote the pertinent exit: 

X.T.Z, respectively. 
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PHC*1 


r.cx •'.oati tvc-ilzco/ sc ruuRC 3-n 

1*X TURK 




1* tern with ipproilMttly 28X orars h oot (highly d n pnndnnt on cowad final position). 
Snooth transient with no signs of aacltatlon of fltilblo dynwics. 
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■. z * kcoj.: 


• v:- k L2cc/ic 


FIGURE 3-11 

1* X T'jRK 



pitch turn took «0out 800 s to complete, c) Position reechos comwM position In 
•bout 260 s, thon overshoots It by -28S. 0) Rato Uk#s about 150 s to «cceler«te to 
comnd turn rot* (0.005*/s), thon overshoots by -32S. 
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1*X TJFS 



u) through j) 

SA defoliation In panrralijed coordinates stowing predominant (and expected) node 4 In- 
plane defoliation of SA, and tone snail coupled out-of-plane notion. Niote steady state 
panel defoliation due to constant linear acceleration of veftlcle produced by ion 
Ur* tars. Extremely low residual vibration levels. 
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roll fLC* no^c: s icOC'G FI6URE 3-13 

accumTiOK 



m) through i) 

list* *t**dy state tolar panel defoliation due to linear vehicle acceleretloo due to 
thrusting. Very lou residual vibration levels. 
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SECTION 4 


GAS JET REACTION CONTROL 

A Gas Jet Reaction Control System (RCS) achieves attitude 
control by selective firing of tiny gas jets suitably located on the 
vehicle. 

RCS systems can provide high control authority at the cost 
of mass expenditures (system weight). As a class, however, they are, 
generally, undesirable for the control of flexible spacecraft because 
their nature is to Impart Impulsive loads to the structure. These 
loads are, essentially, square pulses having very high frequency content 
and, thus, all structural resonant frequencies tend to be excited. 

Because of this, RCS systems frequently exhibit undesirable interactions 
of the control system and the structure, resulting In degraded performance. 

In this section we will discuss two RCS control systems for 
the SEPS/ICM vehicle and will analyze their performance. The section 
Is organized as follows: 
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4.1 . DESCRIPTION OF THE RCS HARDWARE 

To obtain control torques on the vehicle two clusters of 
gas jets will be used as shown In Figure 4-1. 



Figure 4-1. Location of Attitude Control Jets and 
Firing Polarities 


With the above arrangement of jets, it is possible to obtain Independent 
control torques about the three axes by selective firing of the gas 
jets as Indicated In the figure. Note that for pitch (or roll) two gas 
jets are fired providing a pure torque couple on’ the vehicle, while for 
yaw two thrusters are fired which provide torque as well as a transla- 
tional Impulse to the vehicle. This configuration of thrusters has the 
advantage of clustering the gas jets and, thus, minimizing plumbing 
problems. It also has some inheren' disadvantages, for example. 


4-7 
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the yaw jets provide the maximum torque about Y, the axis of least 
inertia, while at the same time producing a translational impulse 
along X which will induce out-of-plane deformation of the solar panels 
and, thus, tend to aggravate the structure/control system interaction 
problem. This configuration may also have plume impingement problems. 
It is thus recognized that this approach may not be the best and, 
although it will be assumed in this study, further analysis and trade- 
offs should be made in the future, to determine if further improvement 
can be attained by trading plumbing complexity against improved perfor- 
mance. 

The gas jets considered in this study were assumed to 
provide a thrust of 0.0675N (15 mlb) with a minimum-on-time (HOT) of 
20 ms, providing a minimum impulse bit of 1.35 10' 3 N-s. The jet 
clusters will be assumed to be located at (0, + 1.27,. +3. 21) m in bus 
coordinates as shown in Fig. 4-1. Since the vehicle center of mass is 
located at (0, 0, + 0.5) m the lever arms for each cluster, in m, are 



Cluster 1 

Cluster 2 

*x • 

1.27 

1 .27 

V 

2.71 

2.71 


1.27 

1.27 


4.1.1 Criteria for Selection of Attitude Control Propellant 

The scientific objectives of the SEPS/ICM mission of 
determining the chemical nature and physical structure of the nuclei, 
atmosphere, and ionisphere of comets places very severe limitations on 
the chemical composition of the propellants which may be used for 

t 

attitude control. For example, we cannot use propellants which (in 
either pure form or by forming other compounds) would compromise the 
scientific objectives of the mission by confusing the background 
concentration ratios. ^ 

Table 4.1 is a list of materials expected to be present in 
the vicinity of a comet. These species are tabulated by mass, and their 

T 7J McMinlmy, W. , Reaction Control Gases for the HFB/T2R Comet Mission, 
I0M 3136-79-167, July 6, 1979, JPL Internal Document. 
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TABLE 4-1. A Search List for Cometary Gases 
Species (1og 10 of Abundance, Normalized to log 10 (H 2 0) a 10) 


Mass Species (aoundance) Mass Species (abundance) 


1 ’ 

0(10) 

It 

, 

2 

0,(1). 0(»> 

11 

0(4) 

J 


40 

a», c-co(7). c*(t>. ax (4) 

A 

0.(4) 

41 

CO, 00(7) 

% 


42 

CO. COd) , KM.CB(l) 

4 



HOCO(7). H,hC0(7) 

7 

■ 

44 

CS(?>. C0,(»). CO,CBO(7) 

• 


4* 

K0.CB-0<6). Ul, -00-01, (4) 

« 


44 

HC0-00(7). MjCS(l), 

co)-oo-oO}(»), »s(»>. 

10 



C ,0,00(7). C0,0C0,(7) 

u 


4 7 


u 

ia C(t) 

4t 

c*(4). »«<*> 

u 

eo<»». u e(t> 

• 49 

C t H(») 

u 

■<7). d,(t) 

*0 

0,0(4) 

u 

*(7). ^BO), 0,(7) 

u 

CN‘C*C»<4» 

u 

o(io). co t (t). »,(7) 

12 

CK4) 

17 

*,(»). oouo) 

11 

C0,-CH*C0(4) 

u 

0,0(10). u 0<7| 

54 


l* 


J! 

144(4) . CO, CO, 01(4) 

20 

B»(4) 

34 

r«(l). C,H,(4) 

U 

) 1 

17 


12 


It 

01(4) 

21 

04(7) 

11 

Co»4) 

24 

C,(l) 

*0 

00(1). CO(BB2>}(4>. 

21 

“c^t). ecioi 


CO,OOCO») 



4i 

00. COOOI4) 

• It 

»(•). (^0,(6) 

42 


27 

oco(0) 




41 


20 

00(10). 0,(1). C J 0 l <4) 

44 

S,(i>. SO, (7) 

21 

OM0(7). OOO(t) 




4) 

Cu(4) 

JO 

0^0(4), 00(7). CjO^O) 

44 

*50,(11 

11 

00,00,(7). 000(4) 

72 

c,o u (i) 

12 

*(•), 00,00(7). *,•*,(?>. 



0,(t) ' ‘ 

71 

O.KUI1I. 0C,0(4> 

11 


44 

0*(4> 

It 

1^(1) 

w 

00.0(4) 

11 


U1 

00,0(41 

» 

0,(7) 

112 

«o(7) 

17 





^ See reference In pace 1-1 
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expected relative abundances to water are shown. In order to not com- 
promise the aforementioned scientific objectives, it is necessary that 
the propellant gases do not introduce into the environment species 
having a relative abundance higher than 7, following rendezvous with 
Tempel 2. 


Examination of the table shows that all gases containing 
Qarbon, Nitrogen, or Oxygen are not acceptable for attitude control 
after Tempel 2 rendezvous. The noble gases (Helium - He, Neon - Ne, 

Argon - A r. Krypton - Kr, and Xenon - Xe) all satisfy the concentration 
constraint. However, there are certain restrictions in their use. 

In particular, if the proposed Gamma Ray Spectrometer is included in the 

OC 

ICM science payload, Kr would not be used because Kr produces 0.514 MeV 
y photons which could alias the y ray observations. The remaining gases 
could be used, although Xe Is less desirable because of its higher 
expected abundance (7), and He is also less desirable because of its 
small molecular size which makes it difficult to prevent from leaking. 

In view of the above, Ne and Ar would be optimum from the 
standpoint of being about 6 orders of magnitude down from water in 
abundance, and are mc r e easily maintained in storage. . Neon appears 
to be a good choice because it provides a specific impulse of I sp =75 s 
which is twice the I jp of Argon. It should be noted, however, that the 
optimum choice depends also on other factors besides I $p (such as 
gas densities, tankage and plumbing weight and, ultimately in total 
system weight). From a systems viewpoint, therefore, Ne may not be such 
an obvious best choice. In fact, it may turn out that the best choice 
will vary from mission to mission depending on whether the total RCS 
system weight is dominated by tankage weight nr by propellant mass 
requirements for the mission. 

For the purposes of this study we will assume a Ne system for 
gas consumption calculations. It should be noted that this assumption 
in no way affects the applicability of other results shown herein should 
another gas be decided upon. Change of propellant will only change 
propellant consumption. 
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RCS LEAD-LAG CONT'OL SYSTEM 


The Lead-Lag RCS control system block diagram (based on loop 
configuration and data supplied by MSEC) is shown in Figure 4.2. . It is 
similar to the lead lag controller described in, Section 3 for TVC 
except that the TVC hardware has been replaced by RCS hardware. Also, 
the control gains have been replaced by deadband logic to fire the 
thrusters. The filtering of the error has been eliminated since it 
was intended to filter the gimbal high harmonics which are no longer 
present. The scan platform and solar array controllers are exactly the 
same as those described in the preceding section for tVC. A list of 
the pertinent system parameters is given in Table 4-2. Most of them 
were described in Section 3, a few new ones, such as the deadband 
characteristics, will be described herein. 

The block diagram shuws the position error being put through 
a lead-lag network to generate the control error e c . This control error 
is then fed intc the deadband logic which determines how the gas jets 
are to be fired: If the control error is less than 0.5° no action is 

taken. If the error exceeds. 0.5° the gas jets are fired to bring the 
vehicle back into the deadband. The scheme proposed by MSFC uses the 
following deadband logic in an attempt to make the system a pseudo- 
proportional bang-bang system. If the control error exceeds 0.5° the gas 
jets are turned on for an amount of time tg N which is proportional to 
the control error as shown in Figure 4-3. For example, at e c » 0.5° 
tQ N a 20 ms, at e c = 0.6° tg N 3 36 ms. and at e c = 1° t QN = 100 ms. 

In this study we have assumed that the RCS logic Is imple- 
mented in ari onboard computer as a routine having a fraimtime of 50 
milliseconds, so that every 50 ms a decision will be made by the routine 
as to how long the thrusters will be on during the next 50 milliseconds. 
Thus, if at some Instant of time the error is 1.5°, for example, the 
logic would say "fire for 180 ms;" in reality, the controller would leave 
the jets on for 50 ms and then perform another logic decision at the next 
50 ms frametime. If the error still exceeded the deadband we would • 
continue firing, if, on the other hand, it had returned into the deadband 
then no further firings would take place. .Therefore, the 50 ms frametime 
decisions override decisions made previously. 


“SCAM 


\ 
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Figure 4-2. RCS "Lead-tag" System- Block Diagram (1 axis shown) 







TASLE 4-2. RCS *L£A0-LA6’ SYSTEM PARAMETERS 
(mm for <11 <ms units* otherwise Indicated) 


PARAMETER 

VALUE 

UNITS 

OESCRIPTION . 

(. IEA3-LAC CONTROLLER 



T, , 

1.0 

% 

Celestial senior tine 




constant 

V 

1S.0 

S' 

Filter tine constant 

T i 

200.0 

s 

Lead-lag 

*2 

10.0 

% 

Lead-lag 

It. OCAOBANO AW 

GAS JET CHARACTERISTICS 


OB 

♦0.5 

d®9 

Deadband size 

PW 

160.0 

M/deg 

Pulse-width Modulation 

Slope: ton • 20 ♦ I60|e_| 
for |e c l»0.5* c 

FTl* 

50.0 

M 

Frequency of execution of 
Jet firing decision routine 

MOT 

20.0 

M 

Gas Jet nlnlnun-on-tlM 

f 

0.0675 

• 

N 

Single Jet force (IS olbf) 

L r 

(0,*1 .27.2.71) 

m 

CM to Jet cluster 1 vector 

4 

(0.-1.27.2.71) 

m 

CM to Jet cluster 2 vector 

III. SOLAR ARRAY AW SCAN PLATFORM 

CONTROL 


E« 

0.02557 

N-«/rad 

Posit Ion -to -torque 9*1n 

T f 

15.0 

s 

Filter tine constant 

T 1 

200.0 

s 

Lead-lag 

T 2 

20.0 

s 

Lead- lag 

. K sp K* p 

1500.0 

N-n/rad 

Scan posltlon-to- torque 
gam 


0.14683 

s 

Clock rate-to-oosltlon gain 

n? 

0.09901 

$ 

Cone rate-to-posltion gain 
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Figure 4-3. Deadband Logic with 20 ms MOT 


The t 0N signals coming out of the deadband logic* are then 
fed Into the Thruster Selection Logic which selects the appropriate / 

gas jet thrusters to be fired as shown In Figure 4-1. It also generates 
the necessary electrical signals to open and close the valves of those 
gas jets. 

The force and torque acting on the bus can then be Computed 
as follows. Let and fj 2 (t),1 * 1,2, ...6 denote the vector 

force pulses produced by each of the 6 jets of clusters 1 and 2 
respectively, all expressed In some common coordinates, say bus 
coordinates. Let Rj and R 2 denote the vector from the bus center of 
mass CMg to clusters 1 and 2 respectively, in bus coordinates. The 
total force and torque produced by the jets on the bus are then given 
by 

FB(t) - l f »(t) + l f 2 (t) 
i-1 1 1*1 1 

TB^t) - R X - | V(t) + R 2 *f V(t) 

. 1*1 1 1*1 1 
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In order for the simulation program to integrate accurately 
this pulse-width-modulated train of torques* one has to use an extremely 
small, integration step. This poses serious numerical and cost problems 
if we consider that simulation times of 500 - 1000 seconds are typical 
due to the slow response of this large vehicle. Because of this* the 
pulse width modulation logic was not simulated directly in that form but 
converted to an amplitude-modulation scheme, such that the amplitude 
of the pulse was made constant over the 50 ms and proportional to the 
error.. Thus, instead of firing 15 mlb for the first, say, 27 ms of a 
50 ms interval, we implemented the proportional logic by firing for the 
full 50 ms with a force of 

5 ® 1 x 1 5 ml b . 

This amplitude modulation approximation makes sense because during each 
50 ms interval it imparts the same momentum to the vehicle as the pulse 
width modulation approach it simulates. 

This completes the description of the RCS lead-lag control 

system. 


4.3 RCS RATE + POSITION CONTROL SYSTEM 

As we will see in the following subsection, some of the RCS 
maneuvers simulated with the lead-lag controller scheme produced 
unsatisfactory results. In an effort to determine if this poor perfor- 
mance could be improved, a second control scheme has been considered. 
This is the so called rate plus position system in which the control 
signal e c is obtained as the negative of the weighted sum of the 
position error and the actual vehicle rate 

e c << 9 S * V + K rp 

The “rate to position" gain K rp provides damping. This 
scheme is similar to the lead- lag system in that the lead-lag produces 
a control error based on the position error plus the rate of the error. 
Our present controller, however, uses the vehicle rate itself .*ather 
than the rate of the position error. The Solar Array controllers were 
also changed to rate * position systems. 
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The resulting block diagram Is shown In Figure 4-4. It Is 
similar to the block diagram shown In Figure 4-2 but the lead-lags have 
been replaced with "rate plus position" signals. The system parameters 
are given in Table 4-3 which reflects a change from the lead-lag 
controllers to rate plus position control systems. 


























TABLE 4-3. RCS ’RATE ♦ POSITION* SYSTEM PARAMETERS 
(um for oil ooot unloss otherwise Indicated) 


PARAMETER 

VALUE 

UNITS 

DESCRIPTION 


I. RATE ♦ POSITION CONTROLLER 




s 

49.96142 

s 

(U to- to -position gain 


t 

1.0 

s 

Celestial sensor tlae 
constant 

I 

11. OEAOBANO AMO BAS JET CHARACTERISTICS 


A 

OS 

*0.5 

deg 

Deadband size 


PUN 

160.0 

as/dog 

Pulse-width modulation 
slope: Con * 20 * 160|e | 
for |e c l»0.5* c 

■/ • ■ 

FTI* 

50.0 

» 

Frequency of execution 
of Jet firing decision 
routine 


NOT 

20.0 

■s 

Gas jet elnlM-on-tlae 


f 

0.0675 

N 

Single jet force (15 a1b f ) 

■ 

L 1 

(0,*1.27,2.71) 

• 

CM to jet cluster 1 vector 

‘ , 

4 

(0,-1.27,2.71) 

■ 

CM to jet cluster 2 vector 

/ 

: 

/ 

/ 

III. SOLAR ARRAY 

AMO SCAN PLATFORM CONTROLLERS 



15.0 

Una/ rod 

SA posl tl on- to- torque gain 



5.53676 

s 

SA rate-to -position gain 



1500.0 

N ho/ rod 

Scan position- to-torque gain 


«s 

.0.14663 

s 

Rate-to-posl tlon gain 



0.09901 

S 

Rate-to-posltlon gain 

4 / 


■V. 


.s' / 
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4.4 RCS PERFORMANCE 

The performance of the lead-lag control systems 
described above was evaluated through computer simulation of the four 
test maneuvers described in 2.4. The bus yaw turn and scan slewing 
maneuvers were simulated also with the rate + position controller 
described in 4.3. Full details of the flexible dynamics models, 
simulation Programs, and assumptions can be found in Section 2. 

The results of the simulations are shown in Figures 4-5 
through 4-10, and are discussed below in some detail. The figures 
themselves are also captioned with pertinent explanations and commentary. 
The nomenclature used in these figures is listed in Table 4-4. 

The 1° pitch turn maneuver simulation (Figure 4-5) 
shows acceptable performance for this type of maneuver. Some struc- 
tural vibration is evident in the rate plots d), h), and, especially, 
in the solar panel deformation coordinate plots u) through x). Note 
structural vibrations to be coincident with gas jet firings. 

The 30° bus yaw turn (Figure 4-6), on the other hand, 
exhibits totally unacceptable performance characterized by uncontrolled, 
undamped limit cycling resulting in a steady state limit cycle gas 
consumption of 540 gr/hour (-1.19 lb/hr) following turn stop. This 
serious situation is not an uncommon one. Several spacecraft hav? 
exhibited this type of behavior in flight. The Voyagers 1 and 2 
are a recent example of exactly this type of situation [8]. This 
phenomenon is thought to occur as a consequence of structure/control 
interactions, where the solar arrays are excited into oscillation and 
their motipn causes (phase-shifted) gas jet firings which sustain the 
oscillations indefinitely. In an effort to determine if this poor 
performance could be improved upon, this same bus yaw maneuver was 
simulated a second time, this time after changing the lead-lag (LI) 
controller to the "rate-plus-position" (R+P) controller described in 4.3. 
The R+P system provided a very clean and reasonable bus yaw turn 
(Figure 4-7) with RCS-typical vibration levels but no evidence of 

^.Tolivar, A.F., "Voyager Oscillations - Some Data and Comments" Inter- 
office Memorandum 343-319, Nov. 16, 1977, JPL Internal Document. 
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control /structure Interaction. Gas consumption was much smaller than 
before and the residual limit cycle rates after the turn are very 
small . 

The multi-axis acquisition maneuver was simulated with the 
baseline lead-lag controller and the resulting plots are shown in 
Figure 4-8. Acquisition (for the particular set of initial conditions 
considered here) takes about 700 seconds and results in unacceptable, 
undamped, high limit cycle rates in yaw, entirely analogous to the 
situation observed for the bus yaw turn. 

Simulation plots for the scan platform slewing test are 
given in Figures 4-9 for the LL controller, and in Figure 4-10,. for 
the R+P controller. In both cases, the vehicle was assuned to be 
initially just "sitting 11 on the edge of the deadband on all 3 axes. 

Both sets of simulation plots appear quite similar. In the sense that 
most of the dynamic disturbance experienced by the vehicle is caused 
by the slewing of the platform, not by the gas jet firings. However, 
a very significant difference is the amount of (Ne) attitude control 
gas consumed for the slew sequence considered; the LL system used 
0.63 gr, while the R+P system used 2.25 gr (almost four times as 
much). This is because the R+P system considered here Is more sensi- 
tive to the rate spikes (produced by the scan slewing and jet firings) 
and results in considerably more multi pul sing when at the edge of the 
deadband. 


In summary, the RCS system has been shown to produce 
vehicle attitude control characterized by high vibration levels and 
excitation of the vehicle flexible appendages. These high vibration 
levels are due to the impulsive nature of the loads produced by the 
nearly-square gas jet pulses. Such loads have a very high frequency 
content and, thus, tend to excite all resonant frequencies In the 
structure. The higher vibration levels are very dangerous and, in 
fact, in 2 out of the 6 simulations, have been shown co lead to severe 
control /structure interaction resulting In exceedingly high limit 
cycle rates, and unacceptable propellant consumption. 
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TABU 4-4. KEY TO VARIABLE NAMES FOR FIGURES 4-S to 4-10 


ALPNl.2,3* 

Angular position of the bus (b Q ), rad 

ANGN 

System angular momentum magnitude, N-rn-s 

COERR) ,2.3 

Control error e c . rad 

COMCL, CONCH 

Position comands Into Scan Platform, rad 

C0NSA1.2 

Position coeMnds Into SA wings, rad 

ERROR! .2.3 . 

Sensed position error e, rad 

/ 

ETA11 to ETA18 

Panel deformation generalized coordinates 
(panel 1, first 8 modes) 

Oil .2,3.4 

Hinge rotation for Solar Array 
Panels 1,2, and Scan Clock, and Cone, 
respectively, rad 

RATE! .2,3 

8us angular rates, rad/s 

TH1 .2,3,4 

Hinge torques for Solar Array Panels 1,2, 
and Scan Clock, and Cone, respectively, rad 

TPONX.Y.Z* 

Gas jet on polarities 

TRKCON 

» c position Input (turn command), rad 

TSC1.2.3 

External torques about vehicle CM, N-m . 

X.Y.ZGASGR A TGASGR 

Neon gas conswptlon In grams (X.Y.Z, and 
total) 

1,2.3, or (equivalently) X.Y.Z. 

denote the pertinent axis. 
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controller of Figure 4-9. 
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SECTION 5 


REACTION WHEEL ATTITUDE CONTROL 

A Reaction Wheel Altitude Control System achieves vehicle 
control through the use of motor driven flywheels. The motors are 
energized to develop the torques which control the vehicle. These 
torques are felt by both the vehicle and the wheels causing each to 
accelerate (In opposite directions) in Inverse ratio to their respective 
moments of Inertia. 

Reaction wheels with suitable control laws can provide 
extremely accurate and smooth proportional control. This Is a very 
Important asset In the control of flexible vehicles as It reduces 
significantly the problems of control/structure interactions that one 
encounters with Impulsive type A/C systems* such as the RCS system 
described In Section 4. 

In the present section we will develop an attitude control 
system for SEPS/ICM using NASA Standard Reaction Wheels. He will also 
evaluate Its performance. The section Is organized as follows: 
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5.1 DESCRIPTION OF THE NASA STANDARD REACTION WHEEL (SRW) 

The SRW (Figure 5-1) consists of a flywheel, a two phase 
induction drive motor, redundant tachometers, a temperature sensor, a 
pressure transducer and a sealed housing. Table 5-1 lists Its perfor- 
mance characteristics. 



Figure 5-1. NASA SRW Cross Section 


Figure 5-2 shows the torques available at various wheel 
speeds and motor voltages for dual-phase excitation of the motor (this 
means that, for a given speed, the torque Is proportional to the square 
of the voltage). It can be observed from Figure 5-2 that In the range 
of ♦ 1800 rpm - and for a given applied voltage - the torque Is approxi- 
mately constant. 

The synchronous speed of the motor (<u s ) at 400 Hz excitation 
Is 2400 rpm. The wheel, however, loses its torquing capability near 
this speed and should not be operated above 2200 rpm.* In this study. 

It Is possible to drive the SRW at different frequencies and thus alter 
Its momentum storage capability, torque characteristics, and maxlmtan 
speed. 
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Table 5-1 . PERFORMANCE CHARACTERISTICS 


Typical Performance 

Description Specification Requlrtnant at Ambient 


20 N-*-s minimum 


Angular Ho— n tun at 24V 

Wheel Speed Range at 
400 Hi Excitation 

Reaction Torque 

Normal Mode at 24V 
High Mode at 33V 

Tachometer Output 


Momentum Vector Mltal I gn — n t 

Torque Noise 

Wheel Breakaway Torque 

Po— r Consumption 

Max lawn Torqulng at 24 V 
Constant Speed (at .9 
Synchronous Speed) 


Zero to «2200 rpm 

.15 Aei minimum 
.30 N-a minimum 

240 Pulses per revolution 
with rotation direction 
Indication 

15 arc nlnute maximum 
±1.5° atitaun above .1 rad/s 
.003 N-ai max la- 

75 watts maximum 
10 watts maximum 


21.8 N-la-t 
Zero to >2200 rpm 

.175 N-m 
.339 N-m 

240 Pulses per revolution 
with rotation direction 
Indication 

2.48 arc minute 

41 . 0 ° above .1 rad/s 

.0023 N-m 

62.4 watts 
6.6 watts 



Figure 5-2. Torque-speed characteristics 



we have aimed at using the SRW at less than 2200 rpm ( 92% of w $ ) and 
have taken advantage of a constant torque approximation over this range 
of speeds, as shown in Figure 5-2. This approximation is given by 

T r * "i. *„* si S" > 5 - 1 > 

where 

= 0.045 in-oz/v 2 (0.3177 N-m/V 2 ) is the "motor 
torque constant," 

a the voltage applied to the motor, in volts, 
sign(V m ) 3 the "sign of V m ". 

Neglecting frictional torques, the SRW can then be modelled 
as shown in Figure 5-3, 



Figure 5-3. SRW model 


2 

0.3177 N-m/V is the motor torque constant, 

2 

0.0868117 kg-m is the flywheel and motor rotor inertia, 
the inertial spin rate of the SRW, in rad/s. 

5.2 THREE-AXES CONTROL IMPLEMENTATION USING SRW'S 

To obtain control torques on the vehicle three SRW’s will 
be used. They are mounted on the bus (b g ) with their spin axes parallel 
to the bus XYZ axes as shown in Figure 5-4. Note that their actual 
location within the bus is imhaterial, they can be placed anywhere they 
fit. 


where 


«m 

J 


U> 
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Figure 5-4. Bus with 3 SRH's 


Let u 3 (u>x « 2 « 3 ) T denote the body fixed angular rates of the 

bus b 

• ® T 

1 > 3 (^ x i> 2 * 3 ) denote the spin angle of SRW's relative to 
b o 

i 3 diagonal (lx I 2 I 3 ) with 1^ the spin axis Inertia of 
rotor 1 

T„ 3 (T T. T denote the vector of applied torques to 
r ri rz ra 

rotors 

TB„ 3 (TB„ TB_ TB„J* denote the vector of torques on the 
r n f 2 r3- 

bus b Q due to SRM's 
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The torques on the rotors and on the bus are given by 



T r = I(* + o>) 

(5-2) 


TB r = - [I(* + o)) + o> * I(* + o))] 

(5-3) 

In our case, for three Identical wheels Ii * I 2 = I 3 * J, these equations 
yield 

T^ 3 J(*i + ) , TB^ = - T rj - J(-«3+2 + 102 * 3 ) (5-4a) 


+ u 2 ) . TB^ 3 - T r2 - J(o)3*j - u>i*3 ) 

(5-4b) 

T |d e ^(*3 

+ W 3 ) • TB n * - T pj - J( -o>2*l + 0)1*2 ) 

(5-4c) 

The terms 

T gi a J(-0)3*2 + 0 ) 2 * 3 ) 

(5- 5a) 


Tg 2 a J(+W3*l, ■* “ 1 * 3 ) 

(5-5b) 


T g3 a j(- 0 ) 2 *l + “ 1 * 2 ) 

(5-5c) 


represent the gyroscopic cross-coupling reaction torques on each axis 
caused by the vehicle and wheel spin rates of the other two axes. These 
gyrpscopic terms are very small under normal conditions, becoming notice- 
able only when the vehicle rates are high and, simultaneously, the wheels 
are spinning at high velocity. 

It should be noted that the tachometers in the SRW’s provide 
measurements of *, the spin rate relative to the SRW housing (i.e., b 0 ). 
To obtain the inertial spin rate o> r all that is needed is to compute 
uy. ° * + hi where * io obtained from the tachometer and o> is obtained 
from the gyros mounted on b Q . 

The final SRW model (including gyroscopic torques) is shown 
in Figure 5-5. 
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■ Figure 5-5. SRW model including gyroscopic torques 

5.3 DESIGN OF THE ATTITUDE CONTROL SYSTEM 

In order to design the Reaction Wheel A/C System for SEPS 
we will make some simplifying assumptions which will allow us to design 

the systems as If the vehicle were rigid. This Is a temporary assumption 
which will only be used for the design phase of the A/C loop and the 
sizing of its parameters. The A/C system developed under these assump- 
tions will then be Integrated with a fully flexible dynamics SEPS model. 
Their joint performance will be evaluated through computer simulation 
of several maneuvers In subsection 5.4. 

Under certain conditions, to be described below, this is 
a viable and desirable design approach because of its simplicity. It 
can be justified by the following considerations. If the control loop Is 
designed so as to provide a control frequency well below. the lowest 
flexible mode frequency, then the vehicle will move slowly enough that 
structural deformations will be minimal and, therefore, vehicle flexibility 
will not play a significant role; the vehicle will move as an essentially 
rigid body. In more mundane and graphical terms, this philosophy can be 
sumnarlzed by the following principle: "If moved slowly enough, even a 

bowl of Jello will behave as a nearly perfect rigid body." 

We should, thus, aim for a control frequency fi c sufficiently 
below the lowest modal frequency of the solar panels - so as to not 
excite the structure significantly - but, at the same time, high enough, 
so that the speed of response of the vehicle is adequate. The lowest 




modal frequency of the panels was given In Section 2 as f^=0.065226 Hz 
S* (0.409827 rad/s, with a period "M5.33 seconds). We will select for n c 
a frequency 10 to 20 times lower than f L ', i.e. , somewhere between 
0.0065 and 0.0032 Hz, let us say 

fl c ■ 0.004 Hz ■ 0.025 rad/s (T = 250 seconds) , 

which is 16 times (4 octaves) below the lowest flexible mode. Such value 
Is considered to be sufficiently below f^ to satisfy the above principle 
and, yet, high enough to provide adequate response times to meet perfor- 
mance requirements. The adequacy of this choice (on both counts) can 
be appreciated by examining the A/C performance plots given in subsection 
5.4. 

The attitude control scheme to be considered here Is a 
"rate-plus-position error" type controller as shown In Figure 5-6. The 
controller commands the SRW's to apply a torque T c on the vehicle which 
Is proportional to the weighted sum of the vehicle rate and Its 

• '■ t 

position error 

T c - - K((e s - e c ) ♦ K rp ^) 



Figure 5-6. .A/C loop using "rate ♦ position error" controller 
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where 


8 C = position conmand 
e 3 vehicle position 

0 3 vehicle position as sensed by a position - sensor - (here 

5 1 

assumed to be -~r) 
xs+1' 

e t 3 0 C - 0 $ . Note that -ei is the position error. 

*2 = ® c ” ® s ■ ^ r p“* Note that -e 2 is the "position error 
plus (weighted) rate error." 

K fp 3 "rate-to-position“ gain (determines damping) 

K 3 control gain T c 3 K e 2 (determines control frequency) 

3 external disturbance torque on vehicle 

I 3 vehicle inertia 

In order to make a preliminary evaluation of the control scheme, let us 
make two simplifying assumptions for the moment: 

a) A perfect celestial sensor (t= 0 so that 0 $ 3 e) 

b) A perfect reaction wheel SRW(s) 3 1, so that the SRW 

can provide a torque equal to its input T c » whatever 
it might be. 

Under these simplifying assumptions, the closed loop transfer function of 
Figure 5-6 (with T d ^ st 3 0 and u 3 w) is given by 




K/I 

s(s + K rp K/TTTK7T 


(5-6) 


that is, a simple second order system with characteristic polynomial 


where 


S 2 ♦ 2£w_ s + w 2 a 0 
n n 

“n <• fl c> 


i K rp 


(5-7a) 
(5— 7b) 
(5— 7c) 


so that a suitable controller can be designed by choice of K and K^. 


In practice, however, SRW(s) is not a "straight wire" but 
contributes its own non-linear dynamics to the system. Such dynamics 
are studied next. 
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5.3.1 Design of SRW Torque Controller 

We will begin by designing the SRW Torque Controller, l.e., 
the controller that takes the desired torque T c as Its input and generates 
a torque TB f on the vehicle which Is reasonably close to T . A method to 
accomplish this Is shown In Figure 5-7. It uses an Integrator followed 
by an accurate speed control loop for the wheel. 


ling 



Figure 5-7. SRW torque control loop 

• i 

The torque command Is first converted to a SRW speed command u by 

Integrating -T c /J with a limited Integrator. The difference between 

the commanded and the actual (Ir.ertlal) speeds Is fed through a compensator 

I^O+k/s) to produce a voltage to be applied to the motor. The square 

root of this voltage Is put through a + 30 V limiter (protection to the 

motor) and then applied to the motor. The tachometer output plus the 

bus gyro rate are added to obtain the SRW Inertial spin rate which 

completes the loop. The square root block Included In this loop has the 

purpose of cancelling the nonlinear square characteristics of the motor. 

Note that the scheme attempts to compensate for the disturbance-like 

effects of the gyroscopic torques T by commanding, not the “desired" 

• y 

torque T £ , but T c ■ T c ♦ Tg Into the loop. This causes T f to track 
- T c , and TB r to track T £ proper. 

There are two "limited Integrators" included In the SW< 
loop. A limited Integrator Is a device which acts as normal integrator 
until Its state (output) reaches a preset limit (LIM) at which point 

S 
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it stops accumulating. Its output stays saturated at LIM until the input 
changes polarity, causing the integrator to come out of saturation 
instantly. 

The first limited integrator in the loop is the one which 
generates the RW speed command and is limited to + LIMi = + 230.38 rad/s 
(2200 rpm) to avoid generating speed commands greater than the capability 
of the SRW, and to avoid the lag problem that would occur if such greater 
commands were generated. 

The second limited integrator is in the feedforward path of the 
compensator. This feedforward path is included so that two problems are 
. taken care of 


a) Maintain SRW speed in spite of its frictional forces. 

b) Maintain a non-zero accelerating voltage input to the 

SRW to compensate for external disturbance torques on 
the vehicle, without the need for a steady state error. 

The integrator is limited also to avoid lag problems. Its limit + LIM 2 
will be set so that the output of the saturated integrator will never 
command more than 30 V at ? c (Figure 5-7) i.e., 

LIM 2 x K r <30 volts. 

This implies 

|LIM 2 | <27.447 rad/s. 

In this study it will be assumed that the vehicle rate 
sensor (rate estimator or gyro) of Figure 5-6 provides a perfect measure- 
ment of the bus rate u, and that the tachometer provides an exact 
measurement of the SRW speed J». Then, for the non-saturating region 
of operation (|V C | <^30 volts, |w c | <_ 230.38 and |a» e ^ | <_ 27.447 rad/s), 
the SRW torque controller can be simplified to that shown in Figure 5-8. 


by 


The transfer function of the SRW Torque Controller is given 

$RM (s). Ms*) . (5-8) 

V s ' s 2 ♦ K(s+k) J 


I 
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Figure 5-8. S mplified SRW torque controller for 
nm-saturating region 


Note that we can place the poles of 
choice of gains and k 

K * 

$(s) a s 2 + (s+k) 

which implies k K 

- 2 _ *r**n 

“n “3 

•r- 2 E “r, t 


the SRW loop anywhere we wish by 

=* s 2 + 2lu> n s + u n 2 , 

(5-9a) 

r . 1 V 

* * ?r • • 

( 5-9b) 

’ 2X * 

(5-9c) 


Note that the zero is automatically fixed at -k. Also note that 
lC » 2^. By trial and error It was determined that a good choice for 
the SRW loop parameters is 


w„ » 0.06 rad/s and C » 1 
n 


That Is, a critically damped SRW loop with w n about 2.4 times faster than 
the control frequency n c * 0.025 rad/s. This choice results in 


k • « 0.03 

2s 


•V ’ 2 * w n ^ 


32.7990066 volts/rad/s 


where we have used J a 0.0868117 kg m 2 and » 0.3177 10’ 3 N-m/V 2 . 
This completes the design of the SRW Torque Controller. 



5.3.2 


Design of Rate + Position Error Controller 


Incorporating the SRW Torque Controller transfer function 
given by equation 5-8 into the block diagram of Figure 5-6, we obtain the 
linearized A/C loop as shown in Figure 5-9. 

The open loop transfer function for the inner loop is 


G^s) Hj(s) 


K K rp K(s+k) 

1 s(s2 ♦ K(sH)) 


(5-10) 


T dist 



Figure '5-9. Linearized A/C loop 


The corresponding root locus is shown in Figure 5-10 (calibrated in 
terms of the dc-gain 

K K_. - , K 

(5-11) 


K . K *ro it k „ t *rp 

* I If k ‘ 


Selecting K^ c = 0.03 yields 


K K rp = 0.03 I , 


(5-12) 


and inner closed loop poles at 

p t =- 0.021307 and -0.04947 + 0.05165j (Ci * C.6917, «i = 0.071519), 
wiwh an inner loop transfer function 





Figure 5-10. Root locus for -Inner loop 
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I K(s+k) 

(s + 0. 02 1 1 307 ) ( s ^ + 2c 1W1 s + m" 7 ) 


(5-13) 


Note that the inner (rate) loop only determines the product 
of the gains K K rp (which is appropriate, since K K rp can be viewed 
as a "rate to torque" gain). 


The outer (oosition) loop is shown in Figure 5-11. 




F i <iu re 5-11. Oli t er 1 oop 


The root locus of this system is shown in Figure 5-12 calibrated in 
terms of the dc-qain 


: . K K 

*dc ' J 

Selecting our operating point at K dc 
(our control frequency >.' c ) at 


— i-r (5-14) 

PW 

= 0.02 yields a dominant pole pair 


-0.0180986 + 0.01 71 644 j (f « 0.725, a • 0.02494 rad/s) 

c c 

so we do indeed obtain a control frequency of about 0.025 rad/s as was 
our intent. 

The other pole pair is at 

-0.042 + 0 . 04 1 5 j (f. = 0.711, u> = 0.059 rad/s) 
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The tot^ue on each solar array wing Is given by 


where 


T a - K (a 
sa s a sa 


“sa ^ * 



» K sa K^ a . 


Such a s>stem can be thought to arise for either of the following two 
mechanizations. 


MECHANIZATION 1) An actuator drive mechanism having control 
.gain K Sa , motor torque constant K^ a , and a 
viscous damper at the output shaft having 
damping constant B = K Sa K^ a K^ a N-m per 
rad/s. 


MECHANIZATION 2) A rate-plus position type controller similar 
to the one described earlier having position- 
to-torque gain ? sa a K sa K^ a and rate-to- 
' position gain K* a . , In this case thedamping 
is provided by the controlloop and no 
damping hardware is required; Instead, it 
requires a measurement (or estimate) of 
the SA rate. 

In either case, the controller results in a second order 
system with ui and given by 

5d So 


'sa 1 


sa 

‘si 


sa 


1 *sa 
I *rp " 


sa 


Using I ■ 234.61 kg nr/wing and selecting 

Sd 


K $a « 15 N-m/rad, Kjj ‘ « 5.53676 s , 


results In 

« 0.2528 rad/s (0.0402 Hz, T a 24.85 s period 1 

Sa 




sa 


0.7 
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This frequency was selected because of the following 
considerations 

1} It is about 10 times faster than the vehicle control 
frequency n c , so that the SA's and their control loops will have 
minimal effect on the vehicle controller. 

2) It is about 1.6 times slower than the lowest modal 
frequency f^ (which corresponds to the first-out-of-plane bending of the 
solar array). However, and more importantly, it is 8.3 times slower 
(-3 octaves) than the frequency of the lowest torsion mode (0.3337 Hz) 
and, therefore, rotational excitation of the SA at. will result in 
minimal excitation of the SA structure. Because of this, the SA control- 
lers will "see" a SA which appears quite rigid. 

5.3.4 System Block Diagram and Parameters 

, Assembling the controllers designed in the preceding pages 
together with the ICM Scan Platform described in Section 2, results 
in the attitude control system block diagram shown in Figure 5-14. The 
system parameters are listed in Table 5-2. 
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TA8LE 5-2. REACTION WHEEL CONTROL SYSTEM PARAMETERS 
(mm for all axes unless otherwise Indicated) 


PARAMETER 

value 

UNITS 

DESCRIPTION 

I. RATE ♦ POSITION CONTROLLER 



t 

1.0 

s 

Celestial sensor tlM 
constant 

K 

107.48293 

N-a/rad 

X-axis position- to- torque 
gain 


2.101620 

N-a/rad 

Y-axis posltlon-to-torque 
gain 


107.48293 

N-a/rad 

Z-axls posltlon-to- torque 
gain 

•w 

49.96142 

s 

Rate- to -posit Ion gain 


II.. SRM TORQUE CONTROLLER 



i 

8.68117 10* 2 

Mm . 2 

SRH Inertia 

LW, 

♦230.38 

rad/s 

SRM speed cound Unit 

k 

0.03 

s' 1 

reed forward Integrator gain 

LIMj 

♦27.447 

rad/s 

Feedforward Integrator Unit 

•v 

32.79006 

V-s/rad 

SRM rate-to-voltage gain 

V LIM 

♦30.0 

V 

Motor voltage Halter 

Si 

0.3177 10" 3 

N-a/V 2 

Motor torque constant 

III. SOLAR ARRAY 

ANO SCAN PLATFORM CONTROLLERS 


15.0 

N-a/rad 

SA posltlon-to-torque gain 


5.53676 

s 

SA rate-to -position gain 


1500.0 

N-a/rad 

Scan posltlon-to-torque gain 

R* 

n> 

0.M683 

s 

Scan rate- to -position galn(clock) 

0.09901 

s 

Scan rate- to -position galn(cone) 
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5.4 SRW ATTITUDE CONTROL PERFORMANCE 

The performance of the SRW control system just described 
Mas evaluated through computer simulation of the, four test maneuvers 
described In 2.4. Full details of the 30-mode flexible dynamics 
.model, simulation programs, and assumptions can be found In Section 2. 

The results of the SRW simulations are shown in figures 
5-15 through 5-1.8, and are discussed below. The figures themselves 
are also captioned with pertinent explanations and commentary. The 
nomenclature used In these figures Is given In Table 5-3. 

Figures 5-15 and 5-16 show the performance for 1° pitch 
turn and 30° bus yaw turn, respectively. Both of these figures 
exhibit similar characteristics in that they both serve to illustrate 
the very smooth proportional control and vehicle response achievable 
with reaction wheels. Excitation of the flexible appendages Is very 
small. 

Figure 5-17 shows the results of the simulation for the 
multi-axes' acquisition maneuver. The initial conditions assumed were: 


Axis 

Initial Position 

Initial Rate 

Corresponding 

Momentum 

pitch 

♦1.0° 

+0.005°/s 

15.62 N-m-s 

yaw 

♦2.0° 

+0.250°/s 

15.27 N-m-s 

roll 

0.0° 

-0.005°/s 

15.58 N-m-s 


The performance plots show a marginally acceptable acquisition 
performance, because the initial conditions used are near the limit of 
the SRW momentum storage capability. This causes the wheels to speed 
up to (or very near to) their saturation point, where they loose their 
torqulng ability. The 20 N-m-s momentum storage capacity of the NASA 
SRW limits acquisitions to cohdltions where the initial momentun to be 
absorbed Is <15 N-m-s. 

For the scan slewing simulation (Figure 5rl8), the SRW 
system was commanded to maintain the vehicle attitude In Its Initial 
zero state while the scan platform was slewing. The box slew sequence 


5-22 


Introduces high frequency disturbances Into the vehicle, with yaw 
being the axis most affected. The control system manages to keep 
yaw within 165 vrad of the Initial position. Mote that the control 
system attempts to preserve the Initial attitude by slowly varying, 
low-frequency smooth nudging of the vehicle towards null, not by trying 
to react to and control each Individual transient caused by platform 
articulation. 

In summary, the SRW system has been shown to provide very 
smooth proportional attitude control of the vehicle, with very small 
excitation of the vehicle flexibility. Its excellent dynamic performance 
is similar to the smooth TVC performance seen In Section 3. It must 
be remembered, however, that all reaction wheel systems are "momentum 
management" systems which have the limitations Imposed by the momentum 
storage capacity (sl 2 e) of the wheels; once the wheels reach their 
maximum speed (saturation) they become "loaded" and must be unloaded 
by a secondary torquing method such as TVC or a gas jet system. Thus, 
even though an SRW system must also have, say, a secondary RCS system, 
the main (and very significant) advantage of SRW system is the mass 
savings to be realized. This is because all zero net momentum maneuvers 
(not exceeding the wheel capacity) may be performed by the wheels with- 
out any propellant consumption at all. This is true for turns, as well 
as for scan platform slewing since the net average momentum on a long 
sequence of typical slews tends to zero, as the number of slews 
Increases. 
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TABLE 5-3. ttT TO VARIABLE NAMES FOR FIGURES S-IS to 5-18 


ALPH1.2.3* 

ANGM 

C0ERR1.2.3 
concL » CCMCN 
EPROR1 ,2,3 
ETA)! to ETA18 

GM1 ,2.3.1 


RATE) .2.3 


TB1.2.3 


TOESAX.T.Z* 


TM) .2.3.4 


TRNCON 
TOM.t.Z 
VWt.T.2 
URMCX.T ,Z 
HRWX.T.Z 


Angular position of the bus (b ), rad 

Systea angular aoaentix* eagn ttude, H-o-s 

Control arror t c> rad 

Position comands Into SA wings, rad 

Sonsad position arror a, rad 

Panel defonaatlon generalized coordinates 
(panel 1, first 9 nodes) 

Hinge rotation angles for Solar Array 
Panels 1,2. and Scan Clock, and Cone, 
respectively, rad 

Bus angular rates, rad/s 

Externa) torques about bus CM (frea SRH's 
plus external disturbances), N-a 

Desired control torques T , N-a 
c 

Hinge torques for Solar Array Panels 1,2, 
and Scan Clock, and Cone, respectively, rad 

* c position Input (turn coimand), rad 

Torques on SRH's, N-a 

SUM no tor applied voltage, V 

SRW speed coaasnd ** c , rad/s 

SRH actual speed « r , ra d/s 


7 


1,2.3, or (equivalently) X.T.Z 


denote the pertinent axis. 


























»r - 


V* 


l •S'. 


Torques and deflections at hinges (solar array and scan platforw) values are negligible. 
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Nagllglbl* Mng* daflactlms 

























f Cl - 


r _L>: MDGlI hlTri RWrt FI6URE 5-16 

30*8U3 : T 'JRh 



c) 1m takas -120 s to accelerate to turn rata (0.25*/s) at uMeh tfee turn Is coaaended 
to stop. It takas another 120 s for vehicle to slow down. d) Shows control error 
conandlng the vehicle to accalerata and decelerate. 
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5 





































600 seconds required tor this perticuUr set of ecquisttion Inltltl conditions. 


5 





X RW stops vehicle and starts to bring It back but RW saturates at T'lOO, loosing its 
torquing capability. Vehicle coasts towards null with wheel saturated until T«530 seconds 
at which tiae control error changes sign and conaands RW to decelerate, causing it to 
coae out of saturation and regain Control of vehicle. See also i). j). o) and r). 
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C: WRWC 


US MRta' 


g i 



o) Saturated pitch 
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'lCa IOjl: «;*r sue FIGURE 5-17 



Me «sc tx- 

T 


u) SA aaalnua hinge angle transient of >0.07*. *) Scan platfom locked In place for 
U»1s none wa r. 
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u) through *) 

Higher viferetion level* in the structure. Prepayment aode i* »1 (out-o* -pltne) 










APPENDIX A 


rt- 


I cr9dfl/tvc-hast7r"1 

PROGRAM THREE AXES TVC CIO FLEX 2 SOLAR PANELS UNOEPLOYEO MAG BOOM 
TWO DEGREE OF 'FREEDOM SCAN PLATFORM IrfOOY B3> 

STOWED HIGH GAIN ANTENNA 
CONDITIONS AT TEMPLE .? RENDEZVOUS 
COMMENT 

REVISED 1C. OCT 79 

METRIC UNITS THROUGHOUT 

BUS DATA FROM H. PRICE'S MOM PROGRAM 

PAfiTL DATA TROM K. GUPTA'S FINITE EL MODEL (SPAR OCT 79) 
THRUST VECTOR CONTROL LOOP WITH LEDLAG OR KRP (NOV 79) 

F J* c, SPECIFIC MANEUVER, PARAMETERS SEE LINE 193 
COMMENT 

INTEGER NC*NF*Ht3,2 >*F (?*3) *PI Cb> •J*K,PLTNRfi *L ,M* N*NN* MANNBR 
ARRAY 


C '•FL V ,( 3 > 

♦ 

CMSLVJ3.3) « 

H!NGLV(3*3) * 

TEMPSTC8) 

• • • 

M'-:< 7) 

• 

MS (3, 7J , 

PB( T» 3) * 

PS( 3*3*3) 


r. ; 9 * 3 ) 

• 

TH < 9 ) , 

TB ( 3 ) * 

IS ( 3 *3 ) 


F n 1 3 ) 

* 

FS ( 3 • 3 ) « 

G M ( 9 ) * 

G MD ( 9 ) 


GMjOf 9 ) 

t 

TIG! 2*6*15) • 

REC(?*6, 15) , 

RF (2*1*3) 

♦ • • 

JF 12*15) 

♦ 

ZE(2*15) * 

T F (2*1*3) * 

FF(2, 1 ,3) 

t • • 

rtl2,l5) 

* 

ETO«2*l3) • 

UO( 3 ) * 

UFHZI2.15) 

♦ • • 

CM t 2*3) 

% 

CMSCLV1 3) * 

SENG1(3) • 

SENG2C3) 

t • • 

E NG1LVI3) 

• 

E%*G?LV(3> * 

RENG1C3) • 

RCNG2(3) 



nrUBLC PRECISION WP0T(7)*ET DO (2*15) 
DATA (HI1*J)«J=1 ,2)/j*l/ 


* LOGICAL NOPLOT t DROP 


(M(2,J),J=1,2>/0,1/ 

• H (3tO) • J = 1 « 2 )/0 . 2/ 
( r (l*J)«jTl*3)/ 1*1*15/ 
°1/G,3,l,l*l/ 


DATA 
DATA 
DATA 
DATA 
CO'MFNT 
COMMFNT 
DATA M«... 

/lHZA.'GA, 2635. 692, 1769.269* 
DATA ( M S(1«J)*J=I*7)... 


* DATA (F (? *J) ,Jsl*3)/2*l* 15/ 


DEFINE THE MASS PROPERTIES OF THE A PODIES 


-50.133* -119. 639* -96.299, 1756.979/ 


/227R9.M, 239.61, 21876. 00* ♦158.02, C.GO, 3.30, 176.818/ 


DATA (MS(2*J),J=1*7)... 

/22389.00, 239.61* 21876.00* 

OATA <M$(3*J),J=1 ,7) . . . 

/13.3 * 7.5 , • 16.5 * »2.1 

COMMENT 

COMMENT POS VEC SPECIF LOCATION OF 
OATA C*!rlV /.0703* .019C, .6266 / 

(CMSLVI l,J)*d=l*2)/ 0.00* ♦19.37695 
<CMSLV<2,J)*o-l ,3)/ 0.0C* -19.37695 
(CYSlVf3,J>*J=l,3>/-1.3&91, -0.370* -0.7893/ 

(HINGLV 11 « J>« J=1 ,3) / a. 00 ,0.00 ,V0.2/ 
<HIN61V(2*J)«J = 1,3)/ C • 9 0 *0.00 *♦0.2/ 
C‘«INGIVT3*J).J = 1,3)/-1 .3091, -0.37Q, -C.7R93/ 


DATA 
DATA 
DATA 
DATA 
f'ATA 
OATA 
FG W MENT 
COMMENT ION 
DATA CNCILV 


-158.32* G.00, 3.00, 176.818/ 


* -0.2 • -1.0 * 90.19 / 

CM'S & HINGES-IN S/C COORDINATES 


*♦ 0 . 2 / 

.♦C.2/ 


ENGINES LOCATIONS - IN 3/C COORDINATES 

/(l . C* ♦(' .5 * 3. 2 1 / l OATA ENG2LV / 0 . 0 *-0 . 5 * 3 . 2 1/ 


COMMENT 

COMMENT mINOE ORIENTATIONS 
DATA <G(1*J)*J = 1«5>/ 0 .CCtGGO*M .0U000C , G. 000000/ 
DATA '6(2*J)*J=1*7>/ C.39COCO*'1.03COCC* 0.000030/ 
DATA »G (3* J) *U=1 *3 )/ 0.000CC3* l< .3 00000 *♦! .COCOOO/ 
DATA <G(9«J),Jrl,2>/ 0 .C C C 0 00 * ♦ 1 .0 0 0 0 C 0 * 0.000000/ 


qS£J! 


/ 


COMMf NT 

COKHfNT RIGID CLASTIC COUPLINC C DC F ► I C I f NT 5 FOB IHl 2 PANCLS ... 

H'A» Wort L °T K. C-U° T A OF OCT 7*5 ... 

CONDITIONS AT T»*PlC. 2 atNOfZtfCuS 
DATA ( <KFC(1«J.K7 ,J=1«67 ,K = 1.1 ?)/... 


11.18700. 

.Ob 995. 

-.C4629. 

-.55259, 

.23969, 

68.52603,... 

-.61909, 

-.0°91°» 

• 0 tJ 7 H2 9 

1.02466, 

-.34126, 

-86.99000,... 

-3 • 13* bC • 

. C 89 2 J . 

-.19753, 

-2.21670, 

.03269, 

-4.5 395 C , •• . 

-.1328?, 

-.0 1 15P • 

9.57620, 

104.930UU, 

-.6598C. 

.76180,... 

.ac 2 -*p. 

•96203 . 

.11763, 

1 ,25 e C 0 » 

•243/4, 

25.2P5C?*... 

1 .*63S? • 

-.0*18* . 

-.C25C3, 

-.24615, 

-1 .37720, 

18 . 396 07,... 

.16*77, 

-. 001 * 8 , 

.00886, 

• 0 ‘j C 6 9 * 

12.1700b, 

1 . 12110 ,... 

.or f,' r , 

' . 2 24 5 7 , 

.01526. 

.11299, 

.15319, 

.54573.... 

• r i 7 6 7 ? , 

.LT97®, 

.01564,, 

. 1 3 b 2 7 , 

A J i 8 *A i | 

1.1520C,... 

• n I s • 


*4 Oi l 5, 

.06512, 

-.47612, 

11.24900,... 

7 393 . 

• C C 8 3 7 « 

.02850, 

.10454, 

.26855. 

-1.5186C ... . 

-.:>3 346 , 


,03038, 

.264 2 *, 

2.25483, 

-2.59760/ 

DATA ( <F FC C 1 ,J,K7 . 

1,6),*' = 13,1 

5 7/... 



.Pl*f °, 

- .0 5626 . 

. JC957, 

.09616, 

§HT 

T #5t**G « •• • 

-.11517. 

-.0 1331 . 

-.03 152. 

-.05192 , 

-.33515, 

-1.13 713,... 

.298,- 7 . 

• r 01 P 2 , 

. I t ' 3 7. 

.00157, 

.23342, 

t. Sb*,-/? / 

''ll t ( (AKCI?*J.K) .J: 

1, 0 ),*=}, 127/.. • 



1 1 . 1 57C ? , 

*.C 695 a . 

.24629, 

-.95239, 

-.33969. 

-68.5260. ... . 

-.61995, 

.29919, 

-.18782. 

1.324,6, 

.74128, 

88,99000. ... 

- 1 . 13*P- . 

-.:■ P9??. 

.19730, 

-2.2167?. 

-. M 2 h 9, 

4.63953.... 


.01158, 

-9.676 25. 

134.931 . 2 , 

.(.59 8 0 . 

-.7 <263, ... 

.«:?• p. 

-.362 C* . 

-. 11 7P 3 , 

1 .25 >1 2 . 

Q ; 4 • 

-?9.266tr ... . 

l • <* - ? f r «. 

.0*»lb4. 

.02533, 

-.24815, 

1.37720, 

-10.596C0 .... 

. 1 *'!• * 7 . 

.0 21 , 


.09169. - 

12.17310, 

-1.12113,... 

. •‘•C* 3. 

- . j ‘ 4 6 7 , 

-.0192%, 

.11299. 

-.13214, 

-.54573.... 

.L 7 F. 7 ? * 

9 79, 


.13027. 

-4.3334. f 

-1.15239,... 

.:*9i v . . 

-•5226C, 

* • 1C ^ 1 5 9 

• 6 1 ? t 

.47812.. 

- 11 .2490C* •• . 

735*. 

-.cre37. 

-• 0 2 660, 

. 1 34 34 . 

-.26895, 

1 . 6 ie» 0 ,... 

-.23 .*5, 

-.ci* 6 r . 

-,C3: 38, 

.1 6 * :?4 . 

*A • *♦ f> 0 f 

2.59780/ 

DATA ll&EC 

f 2 » ») • 8 7 f vl - 

1,67,6=13,157/... 



. 8 1 46 <5 * 

.00625, 

-.CC957, 

.0661 8 « 

-.87656, 

—7 • 364 4 0*... 

-.’1917. 

•21331, 

.03352. 

-.05150, 

.33515, 

1.1571 3,,.. 

3. 

-* 0 r . 162, 

-.jr*>57. 

.T0137, 

-.23342, 

-5.36093/ 

rc g ' r Nf 






•: c cm pan r l r» ► «i utc nc j» s iin hfptzt ant, damping 


TATA IwFHZI 1«J»,D=1 • 

157/... 





, 0.126146 

, C. 1925 0 7 

• !. 22 fcfel? 

• • • • 


: .?» i r 74 

. G. 329296 

« 0.333715 

• C ,35 1 K 29 

t • • • 


:.i7>>3?7 

, .7.393/48 

, 0.410395 

• 0.443883 

• • • • 


^ « 4 fi 3 1 1 H 

, C. 4 78 3 I 6 

• C.6«HR?0 

/ 



DATA < * F HZ ( 2 , J 7 « J = 1 . 

157/... 




'.'.Oh 6226 

, 0.126145 

, C. 142507 

, 3 .226612 

• • • • 


9.26107* 

. C. 329296 

- C. 33 371 3 

• C. 351529 

t • • • 


P.I7.i:;/ 

, 0.393746 

, 0.41 C39 5 

« 0.443863 

«••• 


r . .46*11 H 

, 0.478315 

, 6.521621 

/ 




DATA < Z F (I«.l)«J=l ,157/1 5*9.004/ 
DATA <?f <2 ,J 7, J = 1 ,157/15*0 .CON/ 
DO L? *=j .IN 

VFf 1 ,<0=WFHZC1 ,K> .2..PIF 
. f C2 . K > =WF HZ »2 • « 7 *2 . *P I K. 
12.. rCNTIfiUt 
COPMCNT 



A-3 


CuMPFNT CALCULATE THE S/C CENTER Of MASS 

no cm i k=i«s 

CMSCLV CRJsMrt* 7 >«CNBLV 4 ft »**S 4 1 *M »CMSLV 41 •K|*HS 4 ? » 7 >«CMSLV... 
12 •* » *PS 43«7 >*CMSLV 4 *•«» 

CMSCL V 4 K l=CNSCL¥ 4 * I/ 4 PB 4 M *MS 4 1 • 7 I **!S 42 * 7 > ♦MS 4 3 # Tl > 

C * 1 1 • • COM INUt 
WPITf (b «F 3 ) CMS CL V 

fO» 4 (Nt rO«PuTAT ION Of PH AND PS YCCTORS AMO CH 8 TO f MG «CCTO«S 
00 LI * = 1*5 

*f NCI <«i =ENG1L YCK»-C*HL¥4*» * M NC24* > s| NG2LV4K » -CM8L¥4«f 
c .| *(,! 4*»=fNr.lLY4*»-C*'SCL¥4«» i S t NG? i A » sfN6?L¥ 4« l-CMSCLY 4 « I 
P"4I*» )SH|MU V4!,K»-C*HLV4«> 

P». 4 ;>.«»sHIMCLV 42 f H»-fP«L¥ 4 M» 

P*U A«H):H|NCL ¥4 3 «* »-CHHLV 4 «> 

p S 4 1 • 1 «K ) SH | NGL ¥ 4 1 I • CHSLV41««> 

PS4?*?.a|smi\cL¥42#««>-CMSLV4^*H> 

* SI*-, , ,«»rM|Nr.LV 4 ?#HI-CPSL ¥4 *•«» 

LI.. -ONIIMUC 
►• 1 C*- UP 


1*1 11 tL 

rc*r.f mi 


J 


'M F I N'f 1 «f PU‘ flMO AXES AS FOLLOMS: 


AMSi = 7 4 ION C MG ! Ml PLUHL IXHAUST AXIS I 
AXIS? = Y 4 SOL A* ARRAY AXIS I 
AXIS 1 r, » 


sum TMf PANEUYER HV ITS NUMBER 


1 - T UK N tklim 'SPACECRAFT AriOUT AXIS1 

? - TURN ENTIRE SPACICKAFT AROUT AXIS? 

* - TURN TNT IKE STATECRAFT ABOUT AX I S 3 

A - TURN . hUS ONLY AH&UT AXIS? 

»*’ ILC KTIPINC ARRAYS ON THE SUN 

- ACQUISITION OF CEL REFERENCES 

f - SCAN PLATFORM SLEta OR BOXSCAN 


= 1.0 


CONS r ANT 
r 0 M Hf NT 


PA.NEUV 


A -4 


CONSTANT 
CCNRf NT 


f. J\ST CM 

ro*n v! 


CONSTANT 

NT 


< O^R! N T 


r ir f UN! 
CON *UN! 
ro *'***' nt 


ri/.MAM 

r c M * , r nt 


CONSTANT 

rov**tNT 



S' I INITIAL C.NOITIOWS 


U1S1 


A K I S ? 


AK 1 S3 


•’! 1ST UKHA A-Cf TOKQurS 

TOISTiiJ.Q • !IJlSI?:,ZO'rt-A t 13lST5;0.0 

initial • at t s 

A A T I 1 I - T »P • * A T I ? 1 =0 • 0 • RATfTUC.O 

INITIAL ANGULAR POSITION CRRORF 
ALP^-il^r.P » ALR*?I=0.0 • ALPMM-O.C 

TUmN COR^ANPS 


T»NTI*= PCC.CO 

TKV TkT = H. 7; t t.*l -f • TUN; *» T =A.3 g 3.'?1 - T#T‘0(3M T = 8.7 208 4f 


SOLAR ARRAY 
splAR ARRAY CORRANHS 

SAPSUN = .,.? • SALUCK = 1.0 • SASL! 8=0.0 

SCAN PIATFOH* 

SCAN PLATTORN CORRANPS 

SCA*4CN=r.O « SLfWON = 0.0 * 0OXON =G.O 


• • • 
• • • 
• • • 
• • « 
• • • 
• • • 


• • • 
• • • 
• • • 
• • • 
• • • 
• • • 
• • • 


• • • 
• 0m 


• • • 


• • • 
• • • 


co*wr wt 

C OftPCNT 
co M *r nt 



f 

l 


A-5 


- t 


.r 


•f 


TRN10N=C.O S TRN20N=0.0 S TRN30N=0 . Q 

SAPSUNsf.9 * SC ANON = 0*0 

MANNBR = MANE UV ♦ .5 

CO TO (Ml ♦M2«N3*H4*Mb.N6». MANN HR 


Ml .. 

TRNlONsl .0 
CO TO MCNO 


• 

M2 • . 

T RN2 0N= 1.0 
GO TO Mf NO 



M3.. 

T RN, SONS 1.0 
• CO TO MENO 

$ 

S AL OCK r 1 . 0 


T RN?ON= 1 .0 
r.o TO. mcno 

l 

SAPSL'N=1.0 

N* • . 

CO TO MENO 



MS • . 

SCANON=1.0 

% 

SLE WON = U .0 $ BO XON = 1.0 

ME NO.. 

CONT INU» 




COMMENT 
COMMENT 
COMM» NT 
COMMENT 
COMMENT 


•••••• TVC CONTROLLER PARAMETERS 

SPECIFY 'FIlUR'Skf. LL AD-LAGS* ANO ' KRP PARAMETERS 


CONSTANT 


SLOtaLL 1.0 

f 

F AS T LL - 0. 

0# 

RPPON S 0.0 

• • • • 

N** = 3.*C 

t- 

F« *=39.f,G 


K K X s 39 • Ub 

• • • • 

S < Y s 0 . 1 1 

f 

FRY* 0.77 

• 

RKYsO.77 

»••• 

S*2=1S.7 

t 

F R 2 s2 1 A.Nh 

• 

KK2 =214.96 

t • • • 

T AU* “IN. 



i 

KKRPXr*?.96 

• ••• 

STAulL-209. 

f 

FT AUlL=fcO. 

t 

KRRPVS49.96 

t • •• 

S? Af?L=IC. 

* 

F T AU2L = 1 0 . 

• 

KKRPZS49.96 



ILAC-CNsSLOWLL.F ASTLL 
* < = S< * • SLC.LL ♦ FKY.FAf.TLL 

rY = S-Y*SLOWLL* FKY*FASTLL 

•»/ s , S K 2 »SLOwLL* H/.FASUL 

TAU1LL=STAUU • c LO«LL*FT AU1L«FA5TLL 

T AU' ll = ST AUTL »SLOwLL*F T AO?L*F ASTLL^RPPON*'#^. 
CO*«EVT ClMjAlLLO ENMNE INERTIAS (EACH MV I MB ALLE 0 ENGINE 
h ANN OF 4 ACTUAL fNCINES) 

I&r\Cl=4#0.t?b * IGTNG?s4*«i.4 7b 

**C E NGt = 4*B.‘» % MGENG2=4*8.9 

force of r fngincs = i nimton 

f C : 1 . 

»• TVC CON TROLL! R PARAMETERS END •••••* 


IS A 


COMMENT 

CONSTANT 

r o v Mf NT • « 

COMM! NT •• 

COMMENT «♦ 

CONSTANT 

comment •• 

Cd*»**»vT SPfriFV SA FILTER « 
CONSTANT 


• CELESTIAL 
T AUSsl 

>* SA ANO SP 


SfNSOX parameters 


characteristics 

SftF LEAO-LAGS • ANOKRP PARAMETERS 


• FKSA=2000. 


r TAUT A=<>1. 
F TAU2A=20. 


KKSAs. 32397 


KKRPSAs.N7.6 


S*SA=3.41 
RA=C.01 r 

CONSTANT STAU1A=2C0. * 

STAU?A=?0. • 

TAUf SA=15. 

K r A s S*SA*SLOULL» FKSA*FASTLL 
T AU1SA=STAU1A*SL0VLL*FT AU1 A*F ASTLL 
TAU2SA = STAU2A*.SL0VLL4FT AU?A«F ASTLL*RPPON«99. 
KSAKAsKSA.K A/2. 

COMMFf.’T SOLAR ARRAY SFLECTION 
IF (SALPCK.GT.0.9>SASLFW=C.0 
IF (SAPSUN.LT .0.9JGO TO SA1 


• ••• 





SALOCK = P • 

SASLEW=0.0 
SA1.. CONTINUE 

CONSTANT A1SLBG= G.00 • A1SLEN= 50.00 » A 1 SLRT sfl • 00438 

CONST AWT A?SLBG=. 0.0C • A2SL Eft = 50.00 • A 2SLR T =3 • 0 0 4 38 

COMHf NT 

COMMENT SP CONTROL LOOP J ARAMET E RS AND IC»S 
CONSTANT KCL =1500. • *CN =1500. 

CONSTANT RCL = 1*8.50 • BCN s 223.25 

f 0 " ** C N T SCAN SELECTION 
IMSLruGN.GT.0.9»uOXON = 0.0 
IF <SCA» ON.GT .0.9»GO TO SCI 
SLlWON=C.O 

p^xon 

5-Cl ..CENT INUr 


constant 

CLRATCsO.' 17A5T 

■ t 

C N.R A T E = 0 . 

017*53 


constant 

CLSLPOs20.CG 

* 

CL SLE Ns 

5C.00 


constant 

CNfLH C = 12 L- . £ 

t 

CNSLENs 

150. CO 


constant 

CLB XB 1 =20 .0 O' 

t 

CLE'Xf Q = 

30.00 • 

CLOX WO =9 • 0 

CONST AST 

CLBXB2=5t .00 





constant 

CN«* XBC s31 . GO 

♦ 

CNBXFQs 

15. OC • 

CNBXWQs?. 00 

constant 

B 0 X E N 0 = B ft • f 0 





CONSTANT 

TlMTL=e. CO COO 





CONSTANT 

PIF = 3. 1*1 59285 





COMMENT 
f ;o**“r%T 

SFT INITIAL SS 

• 

filtfn, and 

IL VALUES 



rLSAl T=-ALPH?I . sap SUN 
FLSA2T =-ALP«?I •SAPSUN 
LLSA1I =FLSA1 1 
LLSAPI =FLSA? I 
THJM Ir ftLPHl I 
THTA2I = AL e H? I 
THTAil =ALF'H* 1 
F 1LXI=-TMT»1 I 
F TLYl=-T*tT AJ> I 
f IL71 = -THTA «! 

LL X I =F 1 L X f 
LL * I =T 1LT1 
LL7I=F IL71 
CO’^rKT 

C Q**'*F NT SET UP PLOTTING INFORMATION 
NOPLOT = .FALSE. 

CONSTANT PL TNMR-G 

INTRV 1=10.0 
NEuINTslOO.O 
1NTRV2=50.0 
Sta I TC H=Nf 8 TNT 

CONSENT 
NC =7 4 Nf =? 

NXQOYNsC . 

NXOOER = C: . 

NTOTF. R = C • 

NX.GINT -NXGINT*!. 

CO?*MCNT 

CALL NPCYFL INC. H, KB.HS.PB.PS *G. PI, NF,F,E 1G.REC.RF.WF ,2 FI 

Dt BUG 


* . . . 
* . . . 
* . . . 


ni 




C-f*Q 


COMMFNT 

dynamic 

V*R I A8L l TrTINIU 
Cl RPOY=NXODr R-NXPf RL 
NKOf RL=NXOnC« 

\'VCOYV = NXOOYNM . 

IF (NCPLOTI GO TO NOSAVC 

IF CT.CT .SV1 TCH» GO TO INTI 

YT SPIT=PITNPR*INT«V1 

rr TO st I 

1ST?.. CONT INUf 

CMrCKi =rLTN:4M* INTR VI 

I» (rKF C i«.f*.r.>.*.rw!VTi nfv int=nf • int*intrv? 

Yl C PLT s Nr V 1ST 
SI 1.. CONIT INor 
T J ‘*F T ♦ . i. i >j l 

I» I T IMf P.LI .tf SPt T > CP TO NOS A Vt 
U«!Tl II il l PLTM'R.l 
h!KlTr { ». «F ? ) p<'«F 1 •P 2 «P 3 
• 0 I T r <♦.»'» G'M 1 »*( M« 2 I*GM< 3 I*G*M<.I 

■'O SI ? J=1 ,? 
l : • 

;o sir r si *5 

jf 111 <r T ( J,l ♦ *) ,** = I 

r : .. TOP "A T UIO014.H) 

l =l*‘- 

r ir.. CoNTtNUt 

bkIT r <G«FJ) (( :« <S ,NM ,s\r J « 5» **s| *? » 

F 3 • . FOR-AT 

PI TNPKsPl T*'«S*! 

1‘ « T .(• T .Sri I TCM » PT UlNTsNf blNT«lNTSv:- 

AVF.« CONTI NUT 
I* TT.M.f INI INIGO TO r IN 

OUTPUT ?r.r 


T HNCO 1 « 

T RNC 0? * 

T RNC 03 

f RROR 1 , 

rPR0R2, 

t P R OR 3 

AlPMl • 

ALPM? , 

AlPHi 

RATfl • 

R AT f 2 . 

RATL3 

T SC 1 . 

TSC? « 

TSC 3 

T r > 1 • 

TP? • 

T 1*3 

f. M l « 

IT A 1 1 , 

FTA21 

GM2 , 

I.TA12 » 

F T A 2 2 

GN3 • 

CTA13 • 

ETA23 

G«A • 

TTA1A , 

FTA2A 

T HI , 

FTA1S • 

T. T A 2 N 

T*2 . 

FT A 1 S * 

FTA26 

T H 3 « 

ITA17 * 

ETA27 

ThA • . 

r T A 1 8 • 

FTA2R 

A PGM * 

Of RPCY« 

NXGINT 

COMSA 1, 

COMCL . 

NXODYN 

COMS A? « 

COMCN • 

NXGDFR 

GI“A1 , 

GIMP1 , 

NXQTER 

GIMA2 , 

GIMP2 « 

GNAO 

GO A 10 . 

GMMD • 

G**30 

GMA2D « 

GMH20 • 

GM2D 

COF.RR 1* 

C0CRR2* 

GM10 

COFKK 3 




A -8 


PREPAR ?0 

TDFSR1 
TSC1 
ERROR 1 
AtPHl 
R ATF1 

cm 

GH2 

GH3 

GH4 

TH1 

TH 2 

TH3 

T HA 

GR AID 

ANGM 

COfRRl 

PASGf ... 

TRNC01 
ERROR 1 
ALPH1 
R aim 
TP 1 
I SCI 
G-l 
G m2 
S M 3 
GH% 

T M 1 
TM2 
TH3 
TH A 
A NGN 
C CMS A 1 
CONCL 
GIPA1 
GI**A2 
G *• A 1 0 
G*A?'J 
COERR 1 
C PFRR .* 

DERIVATIVE . CIO 

CINTERVAL CI=P.05 
HORROR ... 

0*10 =1.F-R *GN?P 

G*»50 = 1.C-P «G**AO 

R A T T 1 = 1-.I-S #RATC? 
wr RROR ... 

ph i o =i.t -r *r,«»2n 

GN7D = l.F-8 ♦ GNA 0 
RATC1 =1.F-P *R ATF ? 
NOSORT 

NXQPER =N XQDF R* 1 • 
GMOI1)=GH10 $ 

GHQ< 7 . ) = GH?Q 
GN0<3»=GH3D 
GRD<A)=GHAO 


T0ESR2* 

T0ESR3 

T SC? • 

T SC 3 

ERR0R2, 

ERROR 3 

ALPM2 • 

ALPM3 

RATE? « 

RATES 

ETA11 « 

C0NSA1 

ETA12 • 

C0MSA2 

FTA13 • 

CONCL 

ETA1A • 

COHCN 

ETAlb ♦ 

GIHA1 

FTAlfe , 

6 1 HA 2 

ETA17 * 

GIHfil 

ETA18 . 

G1HB2 

GMA20 • 

G HP 1 0 

TRNCOH, 

GHP20 

COCRR?, 

COEPR 3 

TRNC02, 

TRNC03 

ERROR 2, 

r RROR 3 

ALPH2 , 

ALPH3 

RATE? , 

RATES 

TH? , 

TIJi 

TSC? , 

TSC 3 

r r a 1 1 , 

EIA21 

FTR12 • 

F TA22 

ETA13 , 

C T Ac 3 

F T A 1 * , 

E TA2A 

CTAlb , 

E T A?t 

F T A 1 (• , 

F T A?6 

LTA17 , 

F T A? 7 

CTA18 , 

f T A28 

orRppr, 

NXC1NT 

C OHS A 2 , 

NX 00 YN 

COHCN , 

NXQOI R 

OlHpl • 

NXOTER 

1. 1 h r 2 • 

GNAO 

GHP10 • 

GH30 

G^H ?D « 

GH2D 

C OF HR 2, 

G“1D 


=i .r -8 

♦ G“l 

=l.t -8 

= 1 .E -8 

♦ GN3 

= 1 .E -8 

= 1 .r -f. 

♦RATES 

= 1 .L -t» 

=1 .C-R 

• GH1 

= 1 .E -8 

= 1 .r -8 

♦ GH3 

=i.r -8 

II 

*■* 

• 

n 

l 

(7 s 

♦ RATES 

= 1 . 1-6 


GH < 1 ) = C.H1 

* G u t? ) = GH2 

* GN(3»=GH3 
S GM t A ) = GH A 


• GH? 

♦ G**A 


»GM2 
♦ GNA 


=l.E-8 

=1.08 


=l.E-8 
=1 .E-8 


* • • • 
• • • • 


• • • • 
*••• 


rioi i 

« 1 ) sfTDll 

$ 

till 

• 1) “£ T A 1 1 

£ T0( 1 

.21 -F T 01 2 

A 

£T< I 

.21 sCTA12 

F T 0 ( 1 

.3) SFTD13 

A 

mi 

.3) =FTA13 

CTOU 

• 41 =FT014 

A 

C T C 1 

.41 =£TA14 

F T 0 < 1 

.5) =F10VS 

• 

ETIl 

•SI stTAlS 

F TO ( 1 

.0) sr TD16 

A 

Fill 

•6) sfTAlfc 

F TO < 1 

» 7 > SFT017 


Fill 

.7) =t TA17 

F T D ( 1 

.81 =FTD18 

A 

fill 

♦81 SFTA16 

F T 0 ( 1 

.91 -FTD19 


t T i 1 

.9) sfTAlS 

F TO < 1 

• lOIsFTOUO 

t 

r t<i 

•10l=f TA1I0 

r to ( i 

• 1 1 > = p TD11 1 

! 

FTU 

• III- C TA111 

rT J(t 

.12»=l TU112 


r t 1 1 

• 1 21 =F T A 1 1 2 

F T J ( 1 

* 1 3 1 = F T01 \ 3 

4 

LTIl 

. 1 31 =! Till! 

FT 0(1 

. 1 • JsF TL1 14 

4 

FT(J 

.14I=FTA114 

f T0(1 

.t‘«l = r T01 IS 

t 

r T tl 

.ISlsf TAilS 

f T 0 ( 2 

* I 1 = • T 12 \ 

» 

r t t: 

.11 :FTA21 

FT 0 ( 2 

• 2 ) si .T 02 2 

( 

F T I ? 

.21 SFTA2? 

FT 0(2 

.») si TO? 3 

l 

F T ( ? 

.3) sFTA?! 

F T D < 2 

. <• t s* J P?4 

A 

r 7(2 

.41 s r T A24 

f T 0 < 2 

• S ) st TOPS 

1 

CT (2 

,S| = f T A 2 S 

f'Ol? 

.»•) sF T02F 

4 

F T (2 

«* I =T T A2b 

F T 0 ( ? 

. /» =! T02 7 

1 

F T (2 

4 71 sf TA2 / 

F T 0 ( 2 

. f ? =1 T02P 

1 

r T(2 

.hi s *_ 1 A 2 b 

•>’ T 0 < 2 

• n =t ir 2s 

A 

FT (2 

• 91 sf T A? 9 

r Tu(2 

,ir. >S* T021C 

1 

FTC? 

•101SETA2 JC 

rTj(p 

• 1 1 »sr TC2 11 

4 

F T (2 

.lllsf TA211 

F T 0 ( 2 

• l?ls» T0212 

1 

f. Tt? 

.1 ?!=f TA21? 

F T l) ( 1 

.131 = ' T 0 2 1 3 

1 

r T<? 

.1 »lsr TA2 1 » 

r T 0 ( 2 

• 1* >s‘ TC214 

4 

t Tl? 

.141SETA214 

r T D C 2 

.l‘')srT0215 

4 . 

r i(? 

•1S1=F T A?*S 


WOlllsRATFl * UOtnsRATfr A yO(3ls»ATC3 » ANG*=MN 

C*CRO N<fRO SStllO 0 
-ac«o »r"fri\r u,vp 
yOs<f;j?>-u)/0CSI 
W = ! Ml f. <yO«C(3l 1 
ft isyo<4i«un 
"ACRO If..? 

CC’TNT SCAN PLATFG»»* CCNTRPLL' « 

Cl 'LON-PUlM »Cl SLOt^SPCCC.* ICLSLFK-CLSLHG1.T l«SL£UON 
'“•SlON;p*jLSt ICNStbf 4SCC00.4 (CSSLtN’-CNSIBGl 4TI-SLEWON 
U (T .GT .“"XT NOlbOYGNsO .C 
CLPXON= IPULSF <rLbxni«CLnxF3»CL«Xyn*7)... 

-PULSf <CLBXI ?,CL8XFU«CL3XyO*T)>«eOXON 
f*Mi X 0*. s PUU r.C<CNPX{ J G.CNBXFO.CNf*XW0»Tl*PCX0N 
CO w CL = I \T‘ G ( (CLSLf '.‘♦CLHXON) *CLR*Tr «P,» 

C r *CN= l \TF (,( <CNSLON*CNPXONl.CNRAlt,0.1 
fOH-ENT SOLAR ARRAYS CONTROLLER 
4 1SL0r.r«>ULSF t A 1 SLF f- . b D POt . • C A 1 SL F N - A 1 SL HG1 * T ) *S A5LE V 
CC M S A ! = < -* L PH2 • S AP SUN* I NTF U < Al SL ON • A 1 SLR T» 0 • 1 I • < l .-SALOCK » 
A?SLO\=PULSr ( A 2 S L F - 1 ,SfiOOO.. (A2SLFN-A2SLBG1 .1 l*SASLf.W 
C (<PSA2= ( -AL^ H2 • SAP SUN* 1 NT F C I A? S L ON* A2 SLR T » 0 • I I * ( 1 . -S AL OCK ) 
COM“FNT SOLAR ARRAY HlNGF TCROUFS 
CRRSAl-<C0PSAl-6*’l ) 

FRRSA2= (COPS A2-GH2 ) 

F ILSA1 =REALPL«TAUFSA*tRRSAJ f FLSAlIJ 
F ILSA2 = PEALPL«TAUFSA.ERRSA2«FLSA2I I . 


/ 


A-IU 



SSLEAO LLSA1,FILSA1.LLSA1I«TAU1SA«TAU?SA 
SSLE AO LLSA2,F1LSA2,LLSA2I.TAU1SA,TAU2SA 
F.RPPAl sERRSA 1-KKRPSA*GM1D 
ERPPA2rERRSA2-KKRPSA«GH20 

TH1=KSAKA*LLSA1*LL AGON ♦ KKSA*fRPPAl *R PPON 
TH?=KSAKA*LL£A2»LLAG0N ♦ KKSA*f RPPA2*RPP0N 
COMMENT scan PLATFORM HINGE torqucs 
TH3s-KCL*(GM 3-COHCL)-RCL*GM3D 
TH4 r-KCN* <C“4-C0MCN)-8CN*GM40 
COHPfNT 

COMMENT GENERATE TURN RAPP INPUT SIGNALS 
TRNTR‘J=PULSE C20.0« , 'Ci0CC.C«TRNT!P»Tl 
TRNTOl = I N T F G I TRN1GN»T RN1RT »TRNTRU»G. ) 

TRNC02 = INTF G< TRN73N-TRN7R T*TRNTRU*C. » 

TRNCO ^ =INTL&«TRN50N*TRN3RI *TRNTRU«0.» 

TRN'COM = TRNCOl ♦ TRNCG2 ♦ TRNC93 
Cu«P»NT CELESTIAL SENSORS 

1hrTAl5RrALPL»TAUS*ALPHl*TMTAl I) 

THCT A;sRCAlPL<TAUS*ALPH2«THTAr I I 
THrTA5 = Rr ALPLUAUS,ALrM3*THTA31 > 
f RRPRJ rTRNCOl -THflAJ 

rpR0«?=T«NC02 -TMrTA? 

ERROR >:TRNC.C3 -THE7A3 
COMMENT . 

C OP**l NT •**••• T VC CONTROLLER DYNAMICS SICTIOM 
CuPPYNT COMPUTE GAMMA CO-MANOS 
FIL*=RC ALPL <TAUf..f.‘>KOR UF ILXI) 

F1LY=R» AlPLITAuF, ERROR?, FILTI) 

F II 7 =•» CA LPL < TAUF, ERROR 3,F I L7I ) 

SSLfAO LLX.f RRORl ,LLX I , T AU ILL , T A U?LL 
SSLfAD LLY,»RP0R?,LLYI,TAU1LL*TAU2LL 
SSLEAO LL7, ERROR.', LLZ I , T AUILL , T AU?LL 
E R PPX zF R K OR 1 -R x R P X «R A T f 1 
ERPPY=F F «> OS» i -K*RPY*RATf 2 
CRPP7 = ERRO» 3-KRRPZ*RATf 3 
COERR lsLLX*LLACON ♦ ERPPX*RPPON 
CPrRR?=LLY*LLAr,ON ♦ rRPPY*RP»ON 
COE RK?-LL7*LLAOCN ♦ ESPP7*RPPON 
CtPXC = LL X*K V •LLASO*. ♦ ERkP X *KK * «RPP0N 
0 APVC sLL Y*K V *L L AGON' V ERPP Y «KK V *RPPON 
r A**2C = LL7*K7«LL AGON ♦ ER PP 2 *RR 7 • RPPON 
COHTNT OOPPUTF THE DESIRED TORQUES 
THE SRI = GAMXC •FC*SENC l <31 
T0ESR2 = GAMYC*FC*SFNG1<3> 

TCESR3 = GAP?C«FC*SENG1<2> 

COMMENT COMPUTE &IMBAL ANGLE COMMANDS FROM GAMMA COMMANOS 
G1MA1C=-GAMYC*GAM7C 
GIMA?C= GAMY C*GAM7C 
GIMR1C=-GAMXC 
G1 MBOC- GA M XC 

COMMCNT CC M PUTr GIM8AL OYNAMICS 

GMA10n=GIMAlC*.04-U.?83*GHA10-.C4*GIMAl 
GMA2DQ=GIMA2C*.04-C.283*GNA20-.C4*CIMA2 
GMB10D = GIMHIC* . Q4-C . 28 3*GMR 13- . Q4 *G I MB 1 
GMn2DO=G IMR2C* .04-0 .283*GMp?0-, 04*&IMB2 
GM A 1D = L I Ml NT EGMA 10n, 0, ,-0.1 745,0.1745) 

GMA2D = L I MINT < G m A2D D , 0 • ,-G . I 74 5 , 0 • 1 745 ) 




GNB1D=LININT (GNBlOOtO. *-0.1 7*5 * 0.1745) 

GNR2D=L I BI NT (GN82DD«0. *-G.1745»0.17A5> 

G I NA 1 = L 1 NINT IGNA1 0.0. « -0.5235. 0.52 35 > 

GIHA2=LININT «GN A2D • 0 . • -0 .5235. 0 . 52 35 > 

GINPl=l ININT CGNB10*0.«-0.5235t 0.5235) 
r.IHB2=LI«INT (GNB20.0.. -0.S2 35.0.5235) 

C ONNf NT COPPUTE GlNBALLLO ENGINES REACTION torques ... 

TA IL-MAGS-OOG (TWO) EFFECT OF GINBALLING ENGINES 
TWB1X = -IGF NG1«6MR100 

TWP1Y = -IGCNG1 *<CA|00«C0S<GIHH1 )-GH81D*6NA10«SINC6INBl) ) 
TVD1Z = -lr,lNGl*lG*AlDD«SINlGINBl>*GNBlD*GHA10*C0SC6INBi>> 
TyO?X = ♦IGENG2*6HB?00 

Tw??V r ♦IGF NG2* <G P A29P *COS <G I NB 2 > -GNfe^D *G NA20*S 1 Nt 6 1 NB2 ) ) 
TWO?/ = -IGrNG?*IG**A2DD*SlN<GINB?>*GNP?0*GNA20«C0SCGINB2ll 
TwJX = TWO 1 X ♦ TW02X 
TWDV = T WO 1 Y ♦ TwDOY . 

Twrz = t wc i / ♦ Twnrz 

COPPf NT WRING THRUST UP IN 2C SEC 

fcm*lf=fc« ci -cos ip if *T/20.a ))/* .0 
IF «T.GT. IS. 999999) F C>iALF =FC / 2 . 0 
ry**MO.T CONPuTE FORCF A NO TORQUE ON PUS 
FTlxs-FCHUF *SIN«GINA1 > 
rwit=*FC M ALF •SINIGI-Pl )*CPS<GI"Al ) 
r» i/s-fchalf •cos(G]mbi)*coS(G1na) > 

rv, rx = *FCHXLF •SIN(GIM*2> 

FcrY=-rcH*tF«siN«GlP*2 )«C0C«G1*A?) 

F: ?/t-FCHRLF •CGSIGIN‘-2 >*'COS <G1«A2) 

f f»sF e i y*fb; y 
Ft /SFPI?*FP2? 

T* NGlXS4r«l/«ftfN6M?»-FblY*RELGl<3) 
T‘*‘GlYs-F n 17*NfNGlClli*FB|X*RENGMS> 

TF'.Gl/s^FMl Y*»FNGl < 1) -» H 1 X • RE NG l < 2 ) 

UNO?* :.m?7.RrNG2<2)-FH2Y.R£NG2t3) 


n '.'02 Y = -Ff ?/«RFNG2<l>*FW?X*RENr>?<3) 

T* *.'G2 » = *Fp?YwRrNG2«ll«FP2Y*RCNf2l?> 
T^XsKM.lXMENGPXOWPX 
THV=T» \C.lY*Tf NG?V*TW0Y 
TI«Z = Tf \G l/*TE.>iG?2*TWG/ 

rO^'IfNT Cl^POTE TORGUF ON S/C CFNTFR OF NASS 
T r NG'lX = *FWl /*SFNC1 <2 ) -F 5 1 Y • ST N3 1 <3) 
TFVGlYs-FPlZ^SENClU ) ♦FB 1 X *SENG 1 < 3) 

TFNGIZ = *FH1Y«SFNG1M )-FblX«SENGl <2) 

TE NG?X=*rP27*SFNG2C2>-FH2Y*SENG2C3l 
T r VG2 Y = -F P?Z*SFNG?ll > ♦FB2X *SE NG2 < 3 ) 
TENG27s*FH2V«SF.NG?tl > -F R?X *SENG? ? 2) 
TSCX=TCNG1 V ♦ T F N G ? X ♦TwDX 
TSCY = TF*JG1 Y*TFNG2Y4TW0Y 
TSCZ=TFNG1Z*TF NG2Z*TWOZ 


CQ*PE*.'T TVC CONTROLLER SECTION END 

COP HT NT 
COHPFNT 

COMMFNT •••*•***•#••♦•••••*•*••••••• 

COPPE^T *•••** FLEX OYNANICS ••••*• 


COPHrNT FOHCrS AND torques acting on bus ano hinges 
TP lrTbX^TPISTl $ THCDsTBl 

TH2=TBY.TP1ST2 * TB«2)=TB2 


<V' 


V/ v 


A TBC3>=TB3 


1 H3 -ThZ ♦ TO I S 13 
t *C1 

M' I £ > = r ° Y 

Tr»Cl)-THl 
1nC2 > = TH;> 

TH|3)=TH3 
1m(4) = TH4 

CO“*fVI TORGUCS ON S/C 
T ' Cl = T SC X ♦ TOISH 
TbC.’ = TSCY ♦ TOIST? 
ilC3 - ISC' * FOIST* 
c«"»*»r % r 

Co“*t*'T f C L V E f OR SYSTE P ACCELERATIONS 

call “ * a T C C*. C,TH,TP»TS,Ffi*FS»TF « F F • ON ,0*0, C*00 « E T ,E TO • . 
yo.unuT ,r ior*H“,CR) 

cn*nrM 

CONSENT SVCFE" R A C f S AND POSITIONS 

RAtfi si NTFG tyl'OI c 1 )»R All 1 I » 

Hair? -|MtG.tWPOM?)«RATr; I) 

HAltS s'.Mt G IN DOT t * > <*P AT I 3 1 > 
atrhl :|hli[.(RA!U,ALPHin 
atPH? srvrro cratf?»alf *21.) 

AIP«2 - : NTf r. IK A If »,Al PH.\i > 

P'l>,PlU,P?‘',Pjr>=MC« l‘fCn,p; s Pl,P2,P. ,RATfl,RAU?,RATE3) 


p.; 

= INT‘ GtPi'0,1 .i;CC ) 






Pi 

= r*<t‘ iiPio.r.^.i 







= 1 VT‘ G cp?d,o .cor ) 






P3 

s ! \ 1 r CP 30 • 0 ; C D 0 ) 






r **i o 

vi roc* ooicai.c. 

) 





f V ^ 

-1» T'r. cy^CTCM.o. 

) 





(,* »n 

= 1 N T •' MKf.CUi > ,f . 

) 






: 1 \Tt G CWOOT < 7 1 «C . 

) 






=int: sco-io,' . » 






C-? 

= I NT ( r. (0-20, . . > 






*.«• * 

Sl'.T* g(C*J!)t' •) 






g-a 

= I M F ' GCG-^n.i .) 






*, If'U 

= fNT‘G ctlCDU*l) 

,0. ) 

1 

FTA11 =INTFGIETD11 

,0. 


rni? 

: I N T : G If ICBC 1.2) 

• C.) 

A 

CTA1? =INT r GCETD12 

,0 • 


«■ ni ^ 

: 1 NT' G CF TOO < 1.3) 

,0.1 

t 

FTA13 =1NTFGCLT013 

,0. 


t: t o i a 

= !\TFG CtTOOC 1.4) 

,0. ) 

A 

FTA14 =INTEGCriD14 

,0. 


i ; f d i *• 

=1NTFG Cf TDOC 1,S> 

, 0 . > 

1 

ETAlb =1NTFGCETD15 

,c. 


FT, life 

- I NT r G C F TOO C 1,1 ) 

, C . ) 

$ 

ETAlfa =INTEGCETD1R 

,0. 


r t r* 1 7 

= l NTFC, CF TDOC 1.7) 

, C. ) 

s 

ETA17 =INTFGCETD17 

,0. 


FTOlH 

=lNTFGCt T 00< l.ft) 

♦ 0. ) 

l 

F T A 1 8 = INTF. GCCTD16 

, 0 . 


rxnip 

= I NTF G CF TOO C 1 ,° > 

»C. ) 

s 

ETA19 = I NT E G C E TOl 9 

* o • 


noi io 

- I \T r G C E T DO C 1 .10 ) 

,0. ) 

A 

FTAllOrlNTFGCETOll C 

,0. 


f TO 1 11 

= 1 VTFGCCTDOC 1 ,11) 

f C. ) 

1 

ETA111-INFFGCET0111 

, 0 • 


ETDl 1? 

= INTFGCF TOO C 1 ,12) 

,0. ) 

t 

FTA112=1NTEGCET0112 

,0. 


r toi 1 1 

= 1 KT» G CF TOPC 1.13) 

,C . ) 

I 

FTA113 = :iMTI GCET0113 

,0. 


cmi* 

s 1 NT! G CF TDOC 1 ,14 ) 

,0. ) 

A 

FTA114 = 1NTF. GCETD114 

,0. 


E 10115 

= 1 NT t G CF TDOC 1 , 15 ) 

,0.) 

A 

FTAllb=INTLGCETDllb 

, 0 • 


F 1021 

= 1NTFG CF TDOC ?, 1 ) 

,0. ) 

1. 

ETA21 =INTEGCCT021 

,0. 


r 1022 

= IMFG CC T DO C 2 , 2 ) 

,0.) 

A 

ETA22 = I NT TG CE TD2 2 

, 0. 


F TO? 3 

-lNFrGCLFOOC.' > ,3) 

,0.) 

A 

ETA23 = INTFG CET023’ 

,0. 


t T 024 

sINTEGCFT DOC 2,4) 

,0. ) 

A 

ETA24 =I.VTt G CE7024 

• 0. 
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E lOPb -INTCG(ET00<?»5> *0.1 & ETA25 =I«TE6<ET025 *0.) 

ETP’fe =lMTEG(ETDD(?*f>) *0.) * ETA?b =INTEGIET02b «0.) 

ET027 =!MTfC <Et00<?.7) *0.1 t ETA27 =1NTEG<ET027 «0. > 

ETD28 - 1 NT P G (ET00(?«8> ,0.) S E T 428 =1NTE6<ETD28 ,0.) 

ET029 =1NTEG<ET0DI2«9» . * ETA29 = I NTtG <ET029 «0.) 

ET021S=!MTCG(ETOO(2«10)«0.» * £TA210=INTE6<ET0210«0» > 

FT0211=INT£GtETDDI2»ll)tO»> S E 1 A2 1 1 =1 WTE6 (ET0211 « 0. > 
ETD21?=lNTfG IET00C2. 121.0. > 1 E T A212- I NTFG (E 1021 2*0. > 

ET0215 = 1NTEG(ETOD(2«:3>*0.) \ £ T A? 1 3 = I NT r G < E T02 1 3 ♦ 0 . ) 

f TD21A=lNTf GtETOO<?»lA )«0.) % E T A 2 1 A = l NT EG t E TO 2 1 A , 0 . > 

. £ TD21^= INTTGIE T00(2tlS >«0. > I f T A2 1 3 = 1 NT £ G <E 702 1 5« 0. » 
CO*Nr*T •••••• FLEX OTN AN ICS ENO 

CONNENT 

rsn 

r*;o 

END 

I E NH INAL 

UN.. CONTINUE 
Df p UG 

NXOTrR:NXOTER«l • 

MN 1 Tf I f ;F?»T .NXQINT .NXUD»N*NXODJ R.NXOTER 

f \ r 
r vt 



APPENDIX 3 


8-1 


| CR4BFL/GAS-MASTCR] 

PROGRAM THRfT AXES GAS C1D FLEX 2 SOLAR PANELS UNOEPLOYED NAG BOON 
TWO DEGREE OF FREEDOM SCAN PLATFORM IBOOY 83) 

STOWED HIGH GAIN ANTENNA 
rONOITlONS AT TEMPLE 2 RENDEZVOUS 
COMMENT 

REVISED OR NOV T9 

METRIC UNITS THROUGHOUT 

HUS DATA FROM H. PRICE'S MOM PROGRAM 

PANEL OATA FROM K. GUPTA'S FINITE EL MOOEL ( SPAR OCT 79) 
NEON GAS JETS CONTROL LOOP WITH LEOLAG OR KRP (NOV. 79) 

FOR SPECIFIC MANCUVER PARAMETERS SEE LINE M3 


CC'Mf NT 
IMEGCR 


NC*NF*H(3*2)*F(2»3)«PIIf>)*J*K*PLTN8R*L*M»N*NN* MANNBR 


ARRAY CMbLV(3) 

* CMSLV ( 3* J ) 

*HINGLV( 3*3) • 

T EMPST (8) 

9 • • • 

MH( 7) 

« MSI 3*7) 

• 

PB ( 3 * 3 ) • 

PS (3*3*3) 

9 • • • 

C (4 *3 ) 

* TH ( 4 ) 

• 

Tb ( 3 > • 

TS( 3*3) 

9« • i 

FRI 3) 

• FS ( 3* 3 ) 

% 

GM ( 4 > « 

G MO ( 4 ) 

• • • • 

GMOO(4) 

• EIG(?*6*15>* 

Rf C( 2*6*15) • 

RF (2*1*3) 

• • • • 

UF (2*15) 

« ZF ( 2* 1 5 ) 

• 

TF( 2* 1 * 3) • 

FF (2*1*3) 

9 • • « 

CT(?*15» 

• FTD(2*15) 

t 

WO( 3) • 

wrHZ(2*lb) 

9 • « • 

CK(?*3) 

• CMSC LV ( 3 ) 

t 

SJTCM3) * 

SJTC2( 3) 

9 • • • 

JTC1LVI3) 

• JTC2LV(3) 

9 

RJTCK3) • 

RJTC2(3) 


DOUBLE PRECISION WDOT ( 7) *ETOC(? • 

lb > 

* LOGICAL NOPLOT 

X DROP 



DATA 


(H(2,J) ,J=1*?>/0*1/ 

(H ( * » J ) « J = 1 * 2 > / 0 « 2 / 

(F(l,j),jsl,3)/l,l*15/ 

PI/J,0,1,1*1/ 


DATA 
OATA 
DATA 
OATA 
COMMENT 
OATA M F... 

/ 1 8?4 • ^ 04 * 2635.642 • 1 769.269* 
data (*S(l*J)*Jsl*7)... 

/22084.00, 234.61* 21876.00* 
OATA (MS(2*J)*J = 1 • ?)••• 

/22 Oft A • 00 * 234.61* 21876.00* 
data i MS ( 3* J > * J=1 » 7 ) • • • 

/ 1 3 • 3 • 7.5 « 16.5 * *2.1 

COMMENT 


X DATA (F(2,J),Jsl,3>/?*l« 15/ 


-50.133* -114.634* -46.249* 1756.474/ 
♦158.02* 0.00* 3.00* 176.818/ 

-158.02. 0.00* 0.00* 176.818/ 


• - 0.2 


- 1.0 


• 90.19 / 


COMMENT PCS VEC SPECIF LOCATION OF CM'S t HINGES- IN S/C COORDINATES 
n A T A C»PLV / .070 3 • .0190* .6266 / 

DATA ( CMSL V ( 1*J)«J-1«3 )/ 0.00* ♦19.37645 *♦0.2/ 

OATA ( CMSLVC2.J).J=l*3>/ 0.00* -19.^7645 ,*0.2/ 

OATA <CMSLV( ?,J)«J=1*3)/-1.3691» -0.370* -0.7893/ 

DATA (HINGLVU ,J>*J=1*3>/ 0.09 *0.03 •♦0. 2/ 

DATA («INGLV (2*J)*J=1 *3)/ 0.00 *0.00 9*0*2/ 

DATA ( W INGLV(3*J)«J = 1*3)/-1.3691* -0.37G* -0.7893/ 

COMMENT 

COMMENT JET CLUSTERS LOCATIONS - IN S/C COORDINATES 
DATA JTC1LV /0.0*M .27*3.21/ S DATA JTC2LV /C.0*-1 .27*3.21/ 


COMMENT 

COMMENT HINGE ORIENTATIONS 
DATA (DU *J) »J = 1 *3)/ Q.3C'GGJC«*1 .OCODCD* 0.000300/ 

DATA (G (2* J> • J = 1 *3 >/ O.OOQOOO.M. 000000* 0.000000/ !i * v *-v»F 

DATA (G(3*J)»J = i*3>/ 0.000000'* 0 . 0 00 00 G * ♦ 1 . 0 03 0 00 / * 

OATA (G(<i«J> *J=1*3 )/ 0.C0CG00*'1. 020000* 0.000000/ 

COMMENT 


/ 
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COMMENT RIGID ELASTIC 

COUPLING 

COEFFICIENTS 

FOR THE 2 

PANELS ..^ 

SPAR MODEL 6 Y 

K. GUPTA 

OF OCT 79 



C ONO I T I ONS AT 

TFHPLE 2 . 

RENDEZVOUS 



DATA M»FC< 1 » J.K) 

1 , 6 ) ,K= 1,1 

2 )/... 



-11.16*30. 

.0*995. 

- • C 4* 29, 

-.=5239, 

.339*9, 

68.526QC, 

-.519C9, 

-.0991°, 

.08782, 

1 .02460 , 

-.04128, 

-88.99000. 

-3.1348E. 

• Of 922 , 

-.19750, 

-2.21670, 

.03269, 

-4.53950, 

-.13282, 

-.0 1 158. 

9.5762P, 

1 04.930C0 . 

-.65930, 

.76180, 

,80238. 

.C 6 ? 0 \, 

• .11785, 

1.25900, 

. 340 C 4, 

25.20500, 

1 .4*5.4 * • 

-.7 4164 . 

- • 7 * ^ 3 » 

-.24315, 

-i. 37723, 

10.3960C ,. 

. 158 7 /. 

-.00149, 

.0 3S65. 

.05069, 

12.17030, 

1 . 12110 , 

.P6603, 

.00457, 

.or>:e. 

.11299, 

.13019, 

.54573 , 

.•'76 7". 

.0 09 7°, 

• C 2 ** f 4 , 

.13027, 

4.33850, 

1.15200, 

.54915, 

.022*7. 

• Oo 1 5 , 

« u t» 5 1 2 , 

-.47812, 

11.24900, 

-.77393, 

• 0 '*837. 

.C2‘? r C, 

.10434, 

•26855, 

-1,51880, 

-.2304*:. 

.01450. 

.03038, 

. u 8 4 * 4 , 

2.25463, 

-2.39780/ 

0 / IA l fhT 

Cl 1 ,J.K> ,J = 

1 , l > ,.K = 1 3. , 

10 )/.,. 



.81469, 

-.0 :s25. 

.u0°57. 

.Out IP, 

»87o56," 

7 • 3 6 4 4 j. 

-.11917, 

-.01331, 

-.0*362, 

- . C 6 1 5 0 , 

-.53515. 

-1.13710, 

.79f?* , 

.0 Cl 82, 

. 0 n .. 5. 7 . 

. 0 u 1 3 7 , 

.23342, 

5 .3 59 90 / 

DATA < <firr ( ?, J,K) , . 1 : 

l ,« 1 

?.)/•«• 



-H.1R7.;. 

-.•16 Q 95, 

• 0 4(i 29 , 

-.65239, 

-.33969, 

-68.52600 , 


• 9 r ' 9 1 9 , 

'-. 08782 , 

1.02460, 

.34128, 

86.99030, 

-3 . t 3«* rt " * 

-.06922, 

.15/60, 

-2.21670, 

-.032*9, 

4.53950, 

-.13782. 

t 15?*. 

- 6 . 6 7». 2 0 , 

1 .7 4,93003 . 

.65980, 

-.76180. 


-.762'’ ' . 

-.1 ! 783, 

1 . C 3 9 2 C . 

-. 34 314, 

-25.28533, 

1.46 JR *> 

.041*4, 

• 0 7 6 0 3, 

-.24815, 

1.37720, 

-10.39600 , 

. :bH7 7, 

.PC 146, 

-,r C*85, 

.050*9, 

-17.17000, 

- 1 . 1211 u. 

.***,:;*. 

-.P •'457, 

-• u 1 ''28 , 

. 11299, 

-.13019, 

-.5457?, 

.•:7f. 7?; 

-.0 L°7?, 

— .L,*?*4, 

.13027, 

-4'. 03850 , 

-1.1522?, 

. k AVI 5 , 

-.'. r?to. 

-.07' 16, 

.0 6517 , 

.4731 ?, 

-11.24930, 

-,C 7593. 

-.30857, 

-•C2650, 

. 1 0434, 

-• 26855 , 

1.51880, 

-.rscAS, 

- . J l 457 , 

— • u 3 v 38, 

.0*424, 

-2.26.4t J, 

2.39783/ 

DATA ((FI 

C ( 2 ♦ J t * ) * 4 * - 

1 ,6) . = 1 3 . 

r. )/. .. 



.fllHS 

•• J. v bJ ^ • 

-.10967, 

. j 0 6 1 8 , 

-.37656, 

-7 .3644 J * 

-.1191 7, 

• 0 IJM • 

.03352, 

- • Q51 50 , 

.33515, 

1.13710, 

.29 F 2 ’ . 


- . C C :• 3 7 . 

.00137, 

-.23342, 

-5.35990/ 

rC-^EN’T 






CO** M rNT PA 

?; r is *KC.7ur\cics < l n 

1 1 E K T 7 > AND DAMPING 


data ( yrn? ( l ,j » , j=l , 

16)/... 




•‘.06 t .?r«. 

, :. 12*145 

, 1*.1°253 7 , 0. 22**17 

t • • • 


v . 2*1 0 74 

% J. 329: 96 

, 0.33571 

*. , 0.361 52° 

% • • • 


. 1 ,3 7f A?7 

, O. 39*74° 

, C . 4 1 0 3 ° 

5 ., 0.443683 



0.45311 H 

, 0.4 78 3 1 6 

, 0 . 5 • J 1 8 2 

r / 



DATA <WFhZ(?,J),JM, 

16)/... 




' • 0 o5 21* 

, 0 • 1 7 fc 1 4 5 

, C • l c * i< 

7 , • • 2 2 r. b 1 ? 

t ••• 


: . 2*1 r /«* 

, 0.3292*-£= 

« 0.35371 

* , 0.35152° 

t ••• 


?.57h??7 

• 0.59*746 

, 0.41039 

5 , 0.443Gb* 



0.4521 1 fc 

, 0 . 4 7 H 3 ! ‘ 

. 0.5? 167 

: / 




DATA C. ,r Cl»J)*J = l ,13>/150.( •'!<*/ 
DATA tZrt2»J),J=l,l5)/l5«0.004/ 

an l? k * i • l ! 

jM 1 ,* J = W* H / C 1 ,A> *2. «P I ^ 

W c 1 7 « H ) ?VF H z { t K ) * 2 . «P I r 

L2 • • CONTINUE 
ro^TM 
COMMENT 


— ■ I _ 





SLV<1*K) ♦ MS* 2 »7>*CMSLV... 
MS«2t7»*MS<3*7l» 

NO CUB TO ENG VECTORS 

< K ) =JTC2L V < K)-CHBLV(K) 
2CK)=JTC2LV(K >-CMSCLV<K> 


follows: 

•r 

EXHAUST AXIS ) 
) 


l 


CONSTANT 

COMMENT 


SET SIMULATION TIME 


FINTIH = N.9 


SET INITIAL CONDITIONS 


AXIS1 AXIS? 


CONSTANT 

COMMENT 


CONSTANT 

COMMENT 


CONSTANT 

COMMENT 


DISTURBANCE TORUUTS 
TDIST 1=E *0 • T01ST2=.209E- 

INITTAL RATES 

RAT r lI-l»0 t RATC2I=0.Q 

INITIAL ANGULAR POSITION IRHORS 
ALPH1I = 0.0 • AL PH? I =0 • C 

TURN COMMANDS 


CONSTANT TRNTIMs POP. 00 

CC VST A NT TR N1 PT =e.7 ?Gb4 T-b *TRN2» F=R.36i 3PE 

COMMENT * 


UMBER 


SOLAR ARRAY 


CONSTANT 

COMMENT 


CONSTANT 

CO-MFNT 


CONSTANT 

COGENT 


CONSTANT 

COMMFNT 


CONSTANT 
CONSTANT 
C OMH r NT. 


CONSTANT 

COMMENT 


CONSTANT 

COMMENT 


SET SIMULA! ION TIME 


SET INITIAL 

F1NTIM = S.9 

CONDI TIONS 


AXIS1 

A*IS? 

A X I S3 

OISTURBANCE TORUUES 


TO I ST 1 ? ( .0 * 

Tt'lST2=.20SE-4* 

TOIST 3=0. 0 

INITIAL RATES 

• 


RAT r lI=l' .0 « 

RAT C21-0.0 • 

RATE 31 -0. 0 

INITIAL ANGULAR 
ALPH1 1 = C. Q * 

POSITION ERRORS 
AL PH? 1 =0 • C 

ALPM3I=0.0 

TURN COMMANPS . 



T RN T I M= 2 00.00 




TRN1RT =e.7?f-bA f-S»TRN2R T = R . 36 3 32 F-3 * TR N3R T =8 . 7 266* E-5 

• , 

SOLAR ARRAT 

SOLAR ARRAY COMMANDS 

SAPSUN = i'.0 t ' SALOCK=1.0 * SASLFW=0.0 

SCAN PLATE O-R M 

SCAN PLATFORM COMMANDS 
SCANON=r.O « SLFWON=0.0 • BOXON =0.0 


comment 

COMMENT 

COMMENT 
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COMMENT CALCULATE THE S/C CENTER OF MASS 
DO CM 1 K=1*S 

CHsrLVCK)=MP<7>*CHBLVCKMMSCl*7>*CMSLVCl*K> ♦MS <2 «7>*CMSLV.. 
<2,K)«HSC3»7>*CMSLVl3 f K> 

CMSCLVCK »=CMSCLV <K > / ( MB C 7 » *MS < 1 * 7 1 *HS C 2 * 7 > ♦* SC 3 * 71 » 

CHI.. CONTINUE 
WRITE C6.F3) CMSCLV 

COMMENT COMPUTATION OF PB ANO PS VECTORS AND CMB TO ENG VECTORS 
00 LI K-l *3 

RJTC1 CK) =JTCtLVCK)-CMBLVCK> t RJTC2CK >=JTC2LVCK) -CMBLVCK) 
CJTC1 ( K > ~JTC 1LV CK ) -CMSCLV CK ) & S JT C2 C K ) =JT C2L V C K >-CHSCL V C K ) 
PP<l*K )=HINGLVC ltK»-CMBLV <K> 

PHC2.K >=HINGLVC2fK J-CMBLV <K > 

Pec2*K>=HINGLVC3,K)-CMBLVC.O 
PSC1*1 •K»=HINGLVCl «K>- CMSLVC1*K> 

PSC?.2.K>=H1VGLV C2,K>-C"SLVC?«K> 

PSC3.3fK>=HINGLV I3,K)-CMSLVC3«K) 

CONTINUE 


LI.. 

PICKUP 

INITIAL 

COMMENT 


CONST 8NT 
COM*FNT 


DEFINE THE BUS FIXED AXES AS FOLLOWS: 


A X IS 3 = 2 ( ION ENGINE PLUME EXHAUST AXIS > 
AXIS? = Y C SOLAR ARRAY AXIS ) 

AXIS1 = X 


SELECT THE MANEUVER BY ITS NUMBER 


« 

I - TURN 

ENTIRE 

SPACECRAFT 

AROUT 

AXIS1 

• 

2. - TURN 

ENTIRE 

SPACECRAFT 

ABOUT 

AXIS2 

ft 

3 - TURN 

ENTIRE 

SPACECRAFT 

AROUT 

AXIS3 

si- ft 

A - TURN 

BUS 

ONLY ABOUT 

AXIS2 



5 

ft 


WHILE KEEPING ARRAYS ON THE SUN 
ACQUISITION OF CEL REFERENCES 
SCAN PLATFORM SLEW OR EOXSCAN 


HANEUV = 1.0 


C-cj 


w - v 


IF (SALOCK.GT »C .9)5ASLFW = f>.0 
IF tSAPSUN.L T .0 .9»G0 TO SAt 
SAL OCK = 0 ■ 

SASLFU=P.O 
SA1 ..CONT INUF , 

CONSTANT A 1 S L 8 G = 0.00 * A 1 S L EN = bO.OO • A 1 SLR T =0 . 0 0 A 38 

CONSTANT A2SLBG= 01.00 « A2SLEN- 50.00 « A?SLRT =0. 00438 

COMMENT 

COMMFNT SP CONTROL LOOP PARAMETERS AND I C • S 
CONSTANT KCL =1N,;0. « KCN ' =1500. 

CONSTANT HCL = 148 . SO * BCN = 22C.25 

COMMENT SCAN SELECTION 
IF CSLLWON.GT .0 .NJROXON = 0.0 
IF ESC ANGN.f T .0 .9>C0 TO SCI 
SLFW0N = ::.ii 
prxON = o.o 
sr i ..co-w? tnuf 

CONSTANT CLRATr=0.C17A53 « CNR ATI =0 .0 1 7 45 3 

CONSTANT CLSLl*G= 0.00 « CLSLfN = , 30. OC 

CONSTANT CNSLPC.= 100.00 , CNSLEN= 130.00 . 

CONSTANT CLR*R1= (1 . (> , CLEXFQ= 30 .00 • CLHxyD=9.0 

CONSTANT CLP XH 2 = 15. CO 

constant cnhxpg=h.oc « cnbxfo= 15.00 CNuxypr2.oo 

CONSTANT hC xtNO = f>l;.OC 
CONSTANT T IMTL = 0.7 Ot'OCO 
CONST ANT PJf = 3.1415«»2b5 
COMMFNT 

CO-MINT SET INITIAL SS* FILTER* AND LL VAIUIS 
FLSAll =-ALPH?I *SAPSUN 
FLSA2I=-ALPH2I *SAPSUN 
LLSA1 I =FLSA 1 I 
LLSA21 =FLSA?.I 
THT AT I =ALPf*T’I 
THT A2T =ALPM2I 
T«T AM = AL°MT l 
FILXlf-TMTAM 
F1 LVI=tFHTA 21 
FI17I=-TMTA « I 
LLXI =f I l XT 
LLY I =F ILTI 
LL7I=F1L.!1 
CO**MF NT 

COMMENT SF T UP PLOTTING INFORMATION 
NOPLOT = .rALSE. 

CONSTANT PLTNRR=0 «... 

INTRV1=10.0 «... 

NruiNT=ior.n ,... 

T*TRV2=5t).C 
Sy I TCMsNF WlNT 

COMMTNT 

NC = 3 1 NF = 2 

NXODtN = *i. 

NXOOER=r. 

Nxf}T r H =. • 

NXCl NT =NX01NTM . 

CALL NPPYFL ( NC .H.mR.mS *PH,PS*C. f P I ,NF *F *f 1 G * RFC *RF * VF *2F> 

Of PUG 





N.* 


B-5 


TRN10N=0.0 % TRN20N=0.0 S TRN30N=C.O 

SAPSUN=0.0 i SC A NON = 0.0 

MANNER = MA.NrUV ♦ .5 

CO TO «M1 *N2*M3*N4*M5*N6> * MANNBR 


f H 1 ft ft 


T RN ION-1 .0 
GO TO MENO 

• 





M2,. 


TRN?ON=I.O 
GO TO, MEND 





* 

M3.. 


TRN3ON=1.0 
GO iC Mf NO 

l 

S AL OCK = 1 .0 


- 


M4 . . 


TRN20N=1 .0 
CO TO MENO 

% 

SAPSUN-1 .0 

• 





CO TC ME NO 






M*.. 


SCAN0N=1 .0 

» 

SLE WON = 0.0 S BOXON 

= 1.0 


•F NO • 

• 

CONT iNur 1 






CO**ME NT 
COMMENT 

comment 

ft • • 


• ft 


• •• 



comment 


ft 



C0*“F NT 

SPECIFY F ILTCR.SAF 

LEAD-LAGS* ANO 

KRP 

PARAMETERS 


CONSTANT 

SLOWLL = 1 

.0 

• F ASTLL- O.C 


RPPON = 0.0 

• ft • • 



SKXsO.S 


« FKXS0.5 


KK X -0 . 5 

*.. . 



SKVsO.5 


« FKY=O.S 


KKYrO.S' 

• ••ft 



SKZsO.5 


• FKZsO.b 


KKZ=0.5 

*... 



TAUFS15. 


* 


KKR PXs 2 0.0 

* . . . 



STAU1L=200 

• 

* FT AUIL -60 . 


KKRPYS23.0 

*... 



STAU2L=10. 


« FTAU2LSI0. 


kkrpz=2 0.0 



LLAGON=SLOULL*FASTLL 

<X = SKX*SLOWLL* FKX«FASTLL 

Ft = SKV*SLOWLL* FKY*FA3TLL 

«7 = SN?*SLOWLL* FKZ*F ASTLL 

TAUlLL=STAUlL*SLOyLL*FT AU1L «F ASTLL 

TAU2LL=STAU2L*SL0WLL*FT AU2L*FASTLL*RPP0N*99. 

COFMFNT SPfCIFY OTAOPANDS 0P2C ANfl DblOO WM ICH ARE ... 

THE HTN ON TIME <20 MS* OP ANO Tlir FULL-ON IICQ MS* OB 
CONSTANT DB2HSA. 3633203 • OB 1 00? 8 . 726f.bE -3 

COMMENT GAS JET ADJUSTMENT FOR 2" MS MOT 
FONFJTsO •0*45*3 
CONSTANT ISPsTS. 


GASFACs? .,*1 C00./<9.B1*ISP> 

TTONKlsO. t TT0NYl=0. % TTONZ=0. 

COMMFNT ****** GAS CONTROLLER ° AK AMI TENS END ••*••* 

COMMENT ♦**••••*•«••♦****•••••••**•••**•*•*••••••*•** 

COMMENT **•••• CELESTIAL SrNSOR PARAMETERS ••*•*•** 

... CONSTANT TAUSM.O 

COMMENT •*•••• SA AND SP CHARACTERISTICS ••••••*•• 

COMMENT SPECIFY SA FILTER * S<F LEAP-LAGS * AND KRP PARAMETERS 
CONSTANT SKSAs3.Nl • FKS4=*3.16 « KKSAslB.OOO * 

FAsO.OlP 

s CONSTANT S T AU 1 A r 2C 0 • t FTAU!A=feC. • KKRPSAsb.33 * 

STAU2AS20. • FTAU2A=20. r 

T AUF SA = 1 5 • 

F SA s SF SA *SLOWLL *' F KSA*FASTLL 
T AU IS A=STAU1A*SL0«LL*FTAU1A*F ASTLL 


TAU2SAssTAU?A*SL0WLL*FTAU2A*F ASTLL*RFP0N*99. 

FSAKA=KSA*F A/2. 

COMMENT SOLAR ARRAY SELECTION * 


» 





1 
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FRROR 1 *. 

F RR0R2* 

ERR0R3 « 

• • • 





f 


COERRl. 

COERR?* 

C0ERR3* 

• • • 







ALPH1 , 

ALPH? * 

ALPH3 • 








RATE1 « 

RATF2 • ' 

RATE3 « 

• • • 







T iC 1 • 

TSC2 * 

TSC3 • 

• • • 







FRX • 

FRY • 

FRZ « 

• • • 


, 





TP1 « 

TB? • 

TB3 * 

• • • 







G*1 , 

ETA11 • 

ETA21 « 








G f*2 • 

f T A 12 • 

ETA22 * 

• • • 







GP3 • 

F T A 1 3 • 

ETA23 « 

• • • 







r. mr • 

ftair * 

et».;r * 




' 




T HI • 

ETA15 • 

£TA2b « 

• • • 







TM2 ■ , 

£T A 1 6 * 

ETAIE. • 








TH3 * 

ETA17 • 

E T A? 7 * 

• » • 







THR « 

etair * 

E T A23 « 

• • • 







TRHCO** • 

OFRPOr, 

nxqint* 

• • • 







COHSA1 « 

COHCL • 

NXQOYN* 

• • • 







COPSA'-, 

COHCN « 

NXOOER * 

• • • 







tgascr. 

ANGM * 

NXQTER « 

• • • 





, • v 


XE.ASGR* 

YGASGR* 

Z (j ASGR * 

• • • 







TPONX , 

TPONY * 

TPONZ * 

• • • 







TTONX * 

TTONY * 

TTONZ 







PR F PIN ir 

T SC 1 • 

TSC2 • 

TSC3 • 

• • • • 
• • • 







ERROR 1* 

ERROR2* 

ERROR?* 

• • • 


* 





ALPH1 * 

ALPH2 * 

ALPH3 « 

• • • 





_ 


. RATF1 • 

RATE? • 

R*TE3 * 

• • • 







GM1 • 

FT All * 

COMSA1 . 

• • • 







G*2 • 

ETA12 • 

COMSA2* 

• • • 







GM3 • 

C T A 1 3 * 

CO»CL • 




, 




GMR , 

FTAIR , 

COMCN * 

• • • 







TM1 , 

E T A 1 5 • 

TGASGR* 

• • • 






' 

TH? • 

E T A 1 Ei , 

xgasgk « 

• • • 







TH3 • 

ETA17 • 

YGASGR • 

• • • 







THR , 

FTA18 « 

ZGASGR • 

• • • 





— 


TPONX « 

TPONY • 

TPONZ « 

• • • 







TTONX • 

T TONY * 

TTONZ * 

• • • 







angm , 

TRNCOH* 


• • • - 







COERR 1 • 

C0ERR2* 

COERR 3 







RANGE 

TRNCOl , 

TRNCO? * 

TRNC03* 

• • • 







ERROR 1 « 

FRR0R2. 

ERR0R3* 








COERR 1 • 

C0ERR2 • 

C0ERR5* 

• • • 







ALPH1 • 

ALPH2 * 

ALPH3 , 

• • • 







RATE! * 

RATF2, * 

RATF3 * 

• • • 


' 





T SCI • 

TSC2 • 

TSC3 • 

• • • 



• 




FeX • 

FRY • 

F HZ • 

• • • 







TF1 • 

TR? • 

TH3 • 


* 






G«1 • 

ETAU . 

ETA21 • 

• • • 


• 





G“2 • 

ETA12 * 

ETA22 « 

• • • 







GM3 • 

ETA13 • 

ETA23 • 

• • • 





/ 


G“R • 

ETA1R * 

ETA24 * 

• • • 







T H 1 , 

ETA15 • 

E T A?b * 

• • • 







TH2 « 

FTA16 * 

ETA26 • 








TH3 « 

ETA17 * 

ETA27 * 

• • • 







THR » 

FTA18 • 

ETA28 « 

• • • 




* 



TRMCO*** 

OERPOY* 

NXQ1NT • 



‘ 








/ A m — 







\ 




x i 

' ' ' v 

1 





LNO 

COMMENT 

DYNAMIC 

VARIABLE T=TINIU 
DERPOY=NXOPE R-NXDE Rl 
NXOERL = NXQDE R 
NXOOYN=NXOrYN*l . 

IF (NOPLOT) CO TO NOSAVE 
IF (T.CT. SWITCH! CO TO INT2 
YESPlT = PLTN<>R. I.N T R V 1 
CO TO SU 
1NT2.. CONTINUE 
CHrCKP=PLTNHR« INTRVl 

IF (CHLCKP.GL .P.TUINTI NEVlNTsVE WtNT*INTRV2 
YESPLT : MWINT 
Sll.. CONTI NUT 
TIMFPST*.00CI 

ir ITIM1 P.L T .YESPL T I CO TO NOS AVI 
WRITE. »h.FP> PLTNBR.T 
WRITE (b.F2> P7.P1.P2.P3 
WRITE (b,F2l GM(l!.GN<2).GM<3!«bM(4> 

00 SL? J = 1.2 .. 

L = P 

00 SL2 Ksl.3 . . 

WRIT r (b » F2 ) ITT I J,L«-M) ,Hrl,5> 

E2.. FORMAT (lX.bGM.81 
L = L ♦ 

•N L 2 . . CONTINUE 

WRIT? (b.F3> («CMN*NN>.NN=1.3).N=1 *2) 

FJ.. FORMAT Uv.3GM.fl) 

PLTNHR=PLTN0R*1 

I r (T.r.T. SWITCH! NrWINT = NfWlNT*INTRV2 
NOS A VF • . CONTINUE 
IF <T.GT.FlNTIM!GO TO FIN 
COMMENT CALCULATION CF JET EOUIVALFNT ON TIMf 
PWCALC GAMXC.l'62G«PMX.CI.TP0N* .ORluG 
PWCALC GANYC.OPZP.PWY.CI.TPONY .08 ICO 
PWCALC GAM/C.DB2C.PW7.CI.TPON2.DB100 
COMMENT 
COMMENT 

COMMENT MACRO TO CALCULATE PULSC WI07H 
MACRO MACRO PmCAIC P 

MACRO RELABEL GL1.CL2.GL3 
P(SI=0. 

IF(AFS(P( l) l-PIPI! GLl.6L2.GL2 
GL1 .. P(3>=0. 

GOTO (,L3 

GL2 .. CONTINUE 

Pl.M^.C?-*! AHSIPI 1 >>-P (711*0. GB/«P(bl-PI21> 
IE I PI il.GT.PI*! 1 P(3»=P(*1 

^ IF|P|l).C.T.r.» P(S1=1 . 

IFIPil l.LT.O.l P(5>=-I. 

P(*Sl=P(bl*P(31/P(4> 

CL 3 • • CONTINUE 
MACHO L NO 

OUTPUT 2C0 .... 


TRNCOl. TRNC02. 


TRNC03. 


» • • 


8—1 0 


MACRO REDEF I NE U.UD 
U0=(Q«2)-U)/0(5 ) 

W=INTFG(WD«Q(3) ) r 

Q<1 >=W«-a<4>*WD 
HACRO ENO 

COMMENT SCAN PLATFORM CONTROLLER 
CLSLON=PULSE (CLSLBG*5300G.*(CLSLEN-CLSLBG> *T)*SLEWON 
CNSLON=PULSE ICNSLBG.S 0 00 0 . « (CNSLEN-CNSLBG) *T )*SLEUON 
IF(T.GT.FOXENO>BOXON=0.0 

CLBXONs<PULSCCCLBXBl ♦CLBXFQ.CLBXWO.T) ... 

-PULSE <CLBXB2»CLBXFQ*CLBXW0*T >)*BOXON 
CNBXON= PULSE<CNBXBG«CNBXFQ*CNbXUO«T)*FiOXON 
COMCL = INTEG< <CLSLON*CLBXON) *CLRATE»0.) 

COMCN=I NTFGC <CNSLON*CNBXON) «CNR ATE«0.) 

COMMENT SOLAR ARRAYS CONTROLLER 
AlSLON = PULSE (AlSLHG.bOOOC,* < A 1 SL F N-A 1 SL BG > *T>*SASLLW 
COMSA1 - C-ALPH?*SAPSUN* I NT E G C A 1 SLCN* A 1 SLRT * 0. ) ) *C l.-SALOCK > 

A2SLON=PULSF I A2SL«C .5 0 000 . • < A2SLEN-A2SLRG) *T MSA SLEW 
C0MSA2 = <-ALPH2*SAPi;UN*TNTEGCA2SL0N»A2SLRT» 0. ) > • < 1 . -SAL OCX ) 

COMMENT SOLAR ARRAY HINGE TCRQUES 
ERRSA1=(C0MSA: -GM1 ) 

ERRSA?=<COMSA2-GM2 > 

F1LSA1=F£ ACPL<T AUFCA.ERRSA1 .FLSA1I > 

FlLSA2 s RE ALPL'C TAUF SA.ERRSA? .FLSA2I ) 

SSLE AD LLSAl.rRRSAl ,LLSA1 I . T AU 1 S A ♦ T AU2S A 
SSLEAO LLSA2.FRRSA2.LLSA2I ,T AU 1 SA * TAU2SA 
£RPPA1=FRRSA1-KKRPSA*GM1 0 

ERPPA2srRRSA2-KKRPSA*GM2n , 

,THlrKSAKA*LLSAl«LLAGON ♦ KK SA « E RPP A I *RPPON 

TH2sKSAKA*LLSA?*LLAGON ♦ KKSA*ERP°A2*RPP0N 

COMMENT SCAN PLATFORM HINGE TORQUES 
THJs-KCL* <GM3-C0MCL>-BCL*GM30 
THAt-KCN* (G m R— CO m CN)-PCN*G m AO 
comment 

COMMENT GENERATE. TURN RAMP INPUT SIGNALS 
TRNTRU=PULSt <0.0*b0300.0*TRNTIH,T> 

TRNCOl =INTE G< TRN1 ON* T R N 1 R T *TR N TRU « 0 . > 

TRNC02 =INTE GC TR F'2 ON* T R N2R T *TRNTRU*0. » 

TRNCO? =integ<trn;on*trn3rt*trntru»c.) 

TRNCOM = TRNCOl ♦ TRNCU2 ♦ TRNC03 
COMMENT CFLESTIAL SfNSORS 

thetai=rfalpl< taus.alphi*thtai !> 

THFTA2=REALPL<TAUS»ALPH2.THTA2I ) 

THET A?=REALPL(TAUS.ALPH3»THTA31 > 

ERROR 1 =T RNC 0 1 -THE T A 1 
CRR0R2=TRNC02 -THf TA? 
r.RR0R?=TRNC03 -THf T A 3 

COMMENT ****•******••**•••****•****••**•*•***••*•••*• 

COMMENT *••*•• GAS CONTROLLER OYNAMICS SECTION *#«••• 

COMMENT COMPUTE GAMMA CO"HANOS 
F1LX=RE ALPL <TAUF tERROR I «FILXI » 

FILY=RC ALPL<TAUF*ERROR?*FlLYH 
FIL2=REALPL (TAUF,EKROR3*FIL2I) 

SSLEAO LLX*ERR0RULLXlfTAUlLL.TAU2LL 
SSLE AO LLV*ERR0R2*LLYI*TAU1LL*TAU2LL 
SSLEAO LLZ. ERROR?. L L 2 I ♦ T AU ILL . T AU2LL 

rRPPX=FPROR 1-KKRPX*R ATEI \ 
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\ 

COMSAl , 


COMCL • 

NXQOYN • 

• • • 




C QHSA2 • 


CONCN « 

NXOOE R • 

• • • 




TGASGR • 


ANGM • 

NXQTER • 





XGASGR* 


YGASGR • 

ZGASGRt 

• • • 




TPONX « 


TPONY • 

TPONZ « 

• • • 




TTONX • 


TTONY • 

TTONZ 




derivative ciogas 







ci.ntcrval ci 

= O.Ob 







XERROR ... 








GMin =l. r -7 

• GM20 

= 1 

• E -7 «G.M 1 

=l.E-7 

• GM2 

=1 .E-7 

t • • • 

GM30 =l.r-7 

« GM*0 

= 1 

.E-7 «GM3 

=1 .E-7 

• GM4 

= 1 .E-7 

• • • • 

T TON X =l.C-3 

• T T ON Y 

= 1 

.E-3 «TT0NZ 

sl.E-S 

• . . • 



RATfl = l.t-f 

• RATE 2 

= 1 

.E-fc «R AT £ 3 

=l.E-6 




'•fRRO* ... 








r- mi r> =i. f -7 

* GM20 

= 1 

.E-7 *G* 1 

=1 .E-7 

• GM? 

=1 .E-7 

• • • • 

G-30 M.E-7 

• G M 4 Q 

= 1 

• E “7 «GM 3 

=1 .L-7 

• GH4 

= 1 .E-7 

t • • • 

TTONX = l,.F-3 

•TTONT 

= 1 

.E-3 *T TQNZ 

=l.f-3 

«... 




RATEl sl.C-6 .RATE? =l.E-6 *R ATE3 = l.E-6 
tens OR T 


NXCD£R=NXQDf R* 1. 


GMO(l.» -GMIO $ 

GM<1 > 

zSMl 


GMn(2)?Gt»20 l 

GM<2) 

= GM2 

. : 

G M n<3) = G«*3r s 

CM ( 3 > 

= GM3 ' 

( 

C"0<4)=GM4D S 

GM« «j ) 

= GMA 


= E T P 1 1 

t 

E T < 1 * 1 ) = E T A 1 1 

. 

ETO(i,2» =£TD12 

1 

nil,.?) =E T A 1 2 

• 

ETOO.5) =nPlJ 

S 

ET11*3» =ETA13 


CTDU.41 =ET01«» 

\ 

E T < 1 • A ) =ETA1A 

_ — • 

f TOM «5> = ETOl«i 

« 

ETM,!i> =E T A 1 5 

• 

ET0<1*6> sfTOlfa 

1 

E T ( 1 * t ) rETAlb 

: — , 

ETOU.7) = C TP1 7 

t 

FT11.71 =r TA17 

' . 

etdu.R) = F T D 1 9 

% 

CTI1.8) sfTAia 


ETDT1.9) sETPlM 

Y 

ET C 1 *9) = E 7. A 1 9 


ETD(1.10)=r»011J 

5 

r TtU10l = ETA110 

, 

F 7011*11 >s£. TC111 

& 

ETtl *11 »=ETA111 

< 

Eton « 12 1=1 TCI 12 

$ 

E T 1 1 • 1 2) =C T A 1 12 

i 

f TO 11 . 13>=ETD1 13 

Y 

ET1 1*13I=ETA1 13 


non .ia >=r tdi ia 

1 

ETC l«14)s£TA114 


rT0<l«l‘i|=FT0115 

* 

E T 1 1 • 1 M=E T A 1 15 


ET012.1) =rTD21 

$ 

E T (2 • 1 ) = *TA21 


ETOt2*2> =E T 022 

t 

E T C 2 • 2 » =ETA22 


E70 < 2 • 3 1 T 02 3 

V 

ET12*3> =FTA23 


fTO 1 2 *4 > =» T 0? A 

* 

E T 1 2*4 > =ETA2A 


E»D<2.^» =£702* 

* 

E T 1 2 •* 1 =E T A2* 


£T012*M = FTD2f. 

1 

E T 1 2**»> =ETA2b 


E TO 1 2* 71 =fT027 

t 

E T12* 71 =ETA27 


ETn«2*«l =ETD2R 

\ 

E 712*81 SFTA28 


ET012*9> =rT029 

$ 

E T * 2 * 9 1 =ETA29 

■ - r" ~ 

ETD12*in»=ETP210 

• 

ET12*101=ETA210 


CTP<2«11 >= r TD2ll 

« 

FT12*111=ETA211 


rT0«2*12» = f T 02 1 2 

i 

ET12*121=ETA212 


ETD«?*l'»st T0213 

% 

ET<2*1 3) = E TA213 


CT0«2* J4»sETU21.4 

* 

ET<2*1 A)=CTA21A 


f T012*lM=FT021b 

$ 

ETC2,lb)=ETA21b 


gnil)=RATEl » 

UO (21 

= R A TF2 1 «0131=RATE3 * ANG« = MM 



CG*MENT :5 

MACRO MACRO SSL t AO Q 



\ 
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TH* 1 >=TH1 
THC2>=TH2 
TH<3 1-TH3 
TM*A )=THA 

COMMENT TOROUTS ON S/C 
TSC1. = TSCX ♦ TDIST1 
TSC2 = TSCY ♦ TOTST2 
TSC3 = TSCZ ♦ T0IST3 
COMMENT 


COMMENT SOLVE FOR SYSTEM ACCELERATIONS 

CALL MRATE<NC*TH*TH*TS*FB*FS*TF*FF*GM*GMQ*GMD0*ET*ETD* ... 
WO*WnOT*ETOO*HM*CK> 

COMMENT 

COMEENT SYSTEM RATES ANO POSITIONS 


COMMENT 

RATE1 = I N T LG I W 00 T ( I) « R A T E 1 1 ) 

RATE.? =INTEG(WDOTf.?),RATF?I » 

RATE3 =INTEGCM00TC3».RATE3I ) 

ALPH1 =INTEG<RATE1*ALPH1I) 

ALPM2 =INTEG<RATE2*ALPH2I ) 

ALPH3 =INTEC-tRATE3.ALPH3I » 

PZO*P10*P2P»P30=HCKfHrK*PZ*Pl«P2*P3*RATEl#RATE2*RATE3) 
PZ =INTEG<PZ0*1 .030 > 

PI =INTFG*P10,C.?CO> 

P2 -INTEG * ° 2 D * C .C 00 ) 

*»3 =INTEG*P30.0 .000 1 

GM10 =INTrG<UOOT«*l,?.» 

GM2D sZNTFGIWOOT ( ^ ) * 0 • ) 

GM3D =INTEGfyOOTIh>*0.» 

G*AO =INTEGiW0OT*7»,0.» 

GM 1 sINTEGCGMlD.O. » 

GM2 =INTrG«GM20.t'.> 

GM5 =INTEG«GM3P* 0.1 

GMA =INTCGCGMA0.0.) 


ETD11 =INTCG«ET0DC1.1> 
ET012 ^INT^GtETDDCl*?). 
ET013 =INTEg«ETOOC 1 *31 
ET01A sINTEGCETDOCl.A) 
ETD1* s I NT E0IETDD<1*5) 
ETDlft =INTFG iETDO* Ub) 
CT017 = INTEG<ETOO*1*7> 
FT018 sINTEG«ET,00Cl*6> 
ET019 =INTEGIETD0*1*9> 
FTO110sINTEGCETOO(l*10l' 
CTDl 1 1 s I NTEG IE TDD* 1 *11 > 
ETD112 = INTEG <ETDO* 1*12> 
F.T0I1? = INTEG<ETD0<1*13) 
ET011A=1NTEG*ET00«1 «IAI 
tT0115=lNTFG*ETDDIl*15> 
ET021 =INTEG<ET00I2*1> 
FT022 sINTEG CETOO* ?*2> 
ET923 sINTF. G<ETPD<2*3> 
CTD2A =7NTEG*£TDD*?»A> 
ET025 =INTEG«ETPO<?.5> 
ET02fc sINTEG <ETOO«2**»> 
ET027 =INTEG «ETDD*2*7> 
FTD28 sINTEG(ET00(?*8> 


0.) 

» 

ET All =1N7EGIETD11 

• 0. ) 

0.) 


ETA12 =INTEG*ET312 

• 0.) 

0.1 

s 

ETA13 = I NTEG <E TO 13 

*0. > 

0.) 

s 

ETA1 A slNTEGIETDlA 

*0.1 

0.) 

s 

ETA15 slNTFGiETOlE 

• 0. > 

0.) 

s 

ETAlfe sINTEGtETOlb 

• 0. > 

0.) 

t 

ETA17 = I NT EG <£ TO 1 7 

*0.1 

o.> 

% 

ETA18 =1 NT EG IE TO 18 

• 0. ) 

0.) 

% 

ETA19 rINTEGIET019 

• 0. > 

0.) 

% 

ETA110=INTEG*CT0110 

• 0. 1 

0. » 

% 

ETA111=INTEG*ET0111 

• 0. ) 

0. > 

A 

ETA112=INTEGIETD112 

• 0. ) 

0.) 

* 

ETA113=INTtG*ET0ll3 

*0.) 

0. ) 

t 

ETAHAsINTEGIETOllA 

*0.1 

0. > 

% 

ETA115=INTFG*ETD115 

• 0. > 

0. ) 

s 

ETA21 =INTEG*ET021 

*0. » 

0*1 

s 

ETA22 =INTEG«CTD22 

*0.1 

0.) 

s 

ETA23 =INTEG (ET023 

*0. > 

o.> 

s 

ET A2A SINTEGIET02A 

• o.> 

c.» 


ETA2S =INTE6*ET02b 

• 0. ) 

0. » 

% 

CTA2G sINTEG*ET02b 

*0. ) 

o.> 


ETA27 SINTEGIET027 

• 0. > 

0.) 

t 

ETA28 =1 NT EG IETO 2ri 

• 0. > 
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ERPPY=ERR0R2-KKRPY»RATE2 
\ ERPPZ=ERR0R3-KKRPZ*RATE3 

C0ERR1 =LLX*LLAGON ♦ ERPPX*RPPON 
CCERR2=LLY*LLACON ♦ ERPPY*RPPON 
COERR3=LLZ*LLAGON ♦ ERPPZ«RPPON 
GAMXC=lLX«KX*LLAGON ♦ ERPPX «KK X «R PPON 
- G.A1YC = LLY«KV«LLAGOM ♦ ERPPY«KKY*RPPON 

GAM7C=LLZ*K7*LLAG0N ♦ ERPPZ*KKZ*RPPON 

COMMENT COMPUTE FORCES AND TORQUES 
EP1X=(*TP0NY-TP0N?) «FONEJT 
FH1 Y=0. 

FBlZ=*TPONX*FONEJT 

FR?X=<-*TPONY-*TPONZ>«FONEJT 

FP2Y=b. 

FB?Z=-TPONX*FONE JT 
Ff-V=FB 1 X*FB?X 
FHY=FP1 Y*Fr2Y 
FPZ=FP1Z*FP27 

TJTC1X=*FBIZ*R JTC1 <2 >-FBl Y*RJTC 1 « 3) 
TJTC1Y=-FB1Z»R JTC1<1>*FB1X*RJTC1 <31 
TJTC1Z =♦ FR 1 Y «R JT Cl < 1 1-FP1X*RJTC 1 C 2 ) 
TJTC2X=*FB2Z*P JTC2C2l-FH2Y*RjTC2C3) 

TJTC2V = -F n ?Z*RJTC2 <1 » ♦ C B2X*RJTC2 «3 ) 

TJTC2Z = *Fe2Y*R JTC2 < 1 1 -FB2X *RJTC2 <2 ) 
TPX=TJTC1X*T JTC2X 
TPY=TJTC1 Y*TJTC2Y 
TPZ=TJTC1Z*T JTC2Z 

COMMENT COMPUTE TOROUr AROUT S/C CENTER OF MASS 
TJTC1X=*FB1Z*SJTC1 <2 > -F8 1Y *SJTC 1 <3> 

TJTC1Y = -FP1Z*SJTC1<1MFR1X*$JTC1 <3* 
TJTC1Z=^FRIY«SJTC1 (1 »-F91X«SJTCl<2> 
TJTC2X=*FR?7*SJTC2t2>-FB2Y*SJTC2<3> 
TJTC2Y=-FP2Z*SJTC2C1 > ♦FBPX *S JT C 2 < 3 > 
TJTC?Z=*FB2Y*SJTC2(1)-FB2X*SJTC2<2> 

TSCX=T JTC1 X+TJTC2X 
TSCY = TJTC1V-*TJTC2Y 
TSCZ=T JTC1Z*TJTC2Z 

COMMENT GAS EXPENDITURES 

TTONX=INTEG< ABSITPCNXI ,TTONXl) 
TTONY=INTEGCABS<TPONY> *TTONYI I 
TT0NZ=INTEG«ABSCTP0N7>»TT0NZI > 
XGASGR=GASFAC*FONEJT«TTONX 
YGASGR=GASFAC*FONEJT*TTONY 
ZGASGR=G ASF AC*FONE JT*TTONZ 
TGASGR = XSASGR*YGASCR>ZGASGR 


COMMENT GAS CONTROLLFR SECTION ENO 

V COMMENT *•*••••*•*••••**•••***•«.•**•••*** 
COMMENT 

COMMENT •*♦••*••«**•♦••«*••**•****•* 
COMMENT FLEX DYNAMICS 


COMMENT FORCTS ANO TORQUES ACTING ON BUS AND KlWGcSCfc, 

TPl=TB» ♦TDISTl S TB<1)=TB1 " '* 

TH2=TRYMQIST2 S TB<2>=TB2 

TB3=TPZ»TPIST3 S T«M3>=TR3 

F«C1»sFBX 
FR<2I=FPY 

F“C3)=Fr-Z v 

\ 
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APPENDIX C 


C-l 


[CR4BFL/RWHGYR0SC | 

PROGRAM THREE AXES RUH CIO FLEX 2 SOLAR PANELS UNOEPLOYtO MAG BOOM 
TWO DEGREE OF FREEOOM SCAN PLATFORM (BODY B3> 

STOWED HIGH GAIN ANTENNA 
CONDITIONS AT TEMPLE 2 RENDEZVOUS 
COMMENT 

REVISED 10 OCT 79 

METRIC UNITS THROUGHOUT 

BUS DATA FROM H. PRICE'S "OH PROGRAM 

PANEL DATA FROM K. GUPTA'S FINITE EL MODEL ISPAR OCT 79) 
REACTION WHEEL CONTROL LOOP (SEP 79) 

FOR SPECIFIC MANEUVER PARAMETERS SEE LINE 127 
COMMENT 


INTEGER NC,fjF,HC3, 

?),F C?,3),PI(5),J,K, PLTNBR,L,M,N 

,NN, MANNBR 


ARRAY CMBLVC3) 

, CMSLVC 3,3 ) 

HI NGL V ( 3,3) , 

TEMPSTCB) 

t • •• 

Mf<7) 

, MS ( 3, 7 ) 

PB ( 3 , 3 ) , 

PS (3,3, 3) 

t • • • 

G < 4 ,3 ) 

, T H (4 ) 

TB ( 3 ) , 

TS(3,3> 

* • •• 

FBC3) 

, FS (3,3 ) 

GM ( 4 ) , 

6M0 (4 ) 

• • • • 

GMD0C4) 

• E I G ( 2,6, 15 > 

RECC2 r 6,15), 

RF(2,1 ,3) 

* • •• 

WF 12,15) 

, ZF(2,I5) 

TF ( 2 , 1 ,3) , 

FE(2»1,3) 

«••• 

CT<2,15) 

, ET0C2,.' 5) 

W0C.3) , 

WFHZ (2,15) 



CKC?^) 

DOUBLE PRECISION W DOT ( 7 ) *E1 00 ( 2 • 15 ) 

LOGICAL NOPLOT 
■COMMENT 
DR OP 

DATA (H(l,J),Jsl,2)/0,l/ 

DATA (H(2tJ)*J=l«2)/0«l/ 

DATA (H(3,J),J=l,?)/0,2/ 

DATA (F(1,J),J=1,3)/1*1*15/ 

OATA (F (2«J)«J=1,3)/2,1,15/ 

OATA PI/O, 3*1*1, 1/ 

COMMENT 
OATA MB • 

/1824.3C4, 2635.642, 1769.269, -50.133, -114.634, -46.249, 1756.474/ 
OATA <MS<l*J>,Jsl ,7)... 

/22C84.00, 234.61, 21876.00, ♦158.02, 0.00, 0.0C, 176.818/ 

OATA (MS(2,J),J:l,7).«. 

/22P84.00, 234.61, 21876.00, -158.02, Q.CiO, 0.00, 176.818/ 

DATA (MSI3«U>,J:1,7)... 

/I 3.3 , 7.5 , 16.5 , ♦2.1 , -0.2 , -1.0 , 90.19 / 

COMMENT 

COMMENT PCS VEC SPECIF LOCATION OF CM'S i HINGE S-IN S/C COORDINATES 
OATA C"PLV /.0703, .0190, .6266 / 

DATA ( C MSL V (1,J),J-1,3)/ 0.00, ♦19.37645 ,♦0.2/ 

OATA <CMSLV(?,J),J=1,3)/ 0.00, -19.37645 ,♦0.2/ 

OATA (CMSLV(3,J),J=1,3)/-1.3691, -0.370, -0.7893/ 

DATA CHI NGL V (1«J)«J=1*5)/ 0.00 ,0.00 ,♦0.2/ 

OATA CHlNGLV<2,J),jsl,5)/ 0.00 ,0.00 • ♦G.2/ 

OATA (H1NGLV(3,J),U=1,3>/-1.3691, -0.370, -0.7893/ 

COMMENT 


COMMENT HINGE ORIENTATIONS 
OATA (G(1«J)*U-1,3)/ 0. 000000, ^1. 000000, 0.00C000/ 
OATA <GC?,J),J:l,3)/ 0 . 000 0 00 , ♦ 1 • 0 000 0 0 , 0.000000/ 
DATA ( G ( 3, J ) ,J=1 *3 ) / 0.000000, C • 0 0000 0 , ♦! . 0 00 000 / 
DATA IG«4.J).J=1, 3)/ O.OOOuOO, ♦1.000000, 0.000000/ 
COMMENT 



ft 


£ T 029 --INTFG <£ TDDt2*9> ,0*> 
ET0210=INTFG (f T09<??10)»0.) 
CT0211 = 1 NT EG <E TDO< 2,11 >,C.> 
£ T 02 12=lNTFG<E7DOt?.12»»0.> 
ET021?s!NTEG<ETOOC2.13»tO.» 
ET021* = INTCGtE TDD 12*1 
t TO? 15= l NT t GlETOtMP.lM.O. ) 
COM, "FAIT •••#** FLEX DYNAMICS 
COMMENT «••#.*•****•*»**#*«*. 
END 
F VP 


$ ETA29 =1NTFG IETD29 *0.)‘ 
$ ETA210=INT£Gt£TD210*0. ) 

1 ETA211 = lNTfG(ETD211»0. ) 

% ET A212 = INTf G CFT0212.0. > 
l ETA213=INTEG <ETD213«0. > 

S £T <,21* = INTEG CETD21*»0. ) 

\ ETA215 = INTEG(ETD21S*0. ) 
END ••**** 


r v d 

TERMINAL 

FIN.. CONTINUE 
Pf HUG 

NX(3T£k = NXQTf R* 1 • 

l.Tf I £>,F?)T .NXQ I NT .NXOtl YN « NXG DE R « GO 

rhc 



C-3 


00 LI K = 1.3 

PBClf K > = HINGLV(I*K|-CM8LV (K) 
PBI2*K >=HINGLVC2 »K>-CMBLVIK> 

PB< 3*K >=HINGLV<3»K)-CHBLV CK ) 

PS(1 «1 f K>=HINGLV<l#K>- CMSLVI1.K > 
PS C2 *2 *K>=HINCLVC2*K>-CMSLVC 2*K> 
PS t3*3«K»=HlNGLV<3*KI-CMSLV<3*K I 
LI.. CONTINUE 
PICKUP 
COMMENT 


COMMENT 

INITIAL 

CONSTANT TINITL-O.COOOQO 
CONSTANT PIf =3 .14159265 
COMMENT 


OEFINf THE BUS FIXEO AXES AS FOLLOWS: 

AX IS 3 = Z < ION ENGINE PLUME EXHAUST AXIS >. 
AXIS2 s Y < SOLAR ARRAY AXIS ) 

AXIS1 = X 

»••******• ********* ••••**•*••• •*•••••*•**• ** ■ 

SFLECT THE MANEUVER BY ITS NUMBER 

1 - TURN ENTIRE SPACECRAFT ABOUT AXIS1 

2 - TURN ENTIRE SPACECRAFT ABOUT AXIS2 

3 - TURN ENTIRE SPACECRAFT ABOUT AXIS3 

4 - TURN RUS ONLY ABOUT AXIS2 

WHILE KEEPING ARRAYS ON THE SUN 

5 -ACQUISITION OF CEL REFERENCES 

6 - SCAN PLATFORM SLEW OR BOXSCAN 


MANEUV = 2.0 


CONSTANT 

COMMENT 


C-2 


COMMENT RIi-IO ELASTIC COUPLING COEFFICIENTS FOR THF. 2 PANCLS ... 

SPAR MOOEL fit K. GUPTA OF OCT 79 ... 


\ CONDITIONS AT TCMPLE 2 RENDEZVOUS 

DATA ( <REC(1 . J*K> «J = I , 6 ) ,K = 1,12>/... 


-11.18700, 

.06995, 

— .04629, 

-.55239, 

.03969, 

68.52600 ,»•« 

-.51909, 

-.09919, 

.08782, 

1 ,024b0 , 

-.04128, 

-88.99000, .. . 

-3.13480, 

.08922 « 

-. 19750, 

-2.21870, 

.03269, 

— 4.53950, •• . 

-.13282, 

-.91158, 

9.57620, 

194 .93000 , 

-.65980* 

.76180 ,,... 

.80206, 

.06203, 

.1 1 783, 

1.25900, 

.04004, 

25.28500 ,. . . 

1 .46380 , 

-.041b4« 

-.025C3, 

-.24815, 

-1.37720, 

,10.39600 «.. . 

.15877, 

-.00148, 

• 0 C 885 , 

.05069, 

12.17030, 

1.12110,.. . 

. 066; 3 , 

.03457, 

.01928, 

.11259, 

.13013, 

•54573,... 

.37672, 

-00979, 

.02564, 

. 1.302 7 , 

4.03850, 

1 . 1 520C , . . . 

.54915, 

.02260, 

.00615, 

.06512 , 

-.47812, 

11.24900*... 

-.07 393., 

.00837, 

.02850, 

.10434, 

.26855, 

-1.51860, .. . 

-.23045. 

.01450, 

. C 3 C 38, 

.08424 , 

2.25460, 

-2.39780/ 

DATA URECU ,J«K) «J=1 

,6) «K=13, 

15)/... 



.81469, 

-.00625, 

. 0 C 9 57 , 

.00618, 

.67656, 

7.36440,... 

-.1191 7, 

-.01331 , 

03352, 

-.05150, 

-.33515, 

— 1 .13710, *. » 

•29823. 

.00182, 

•00037, 

.00137, 

.23342, 

5.35990/ 

data <<pec< 

7 « J, K ) , J= 1 

t6i 9 K— 



-11.18700, 

- .C 6 q 95 , 

.04629, 

-.55239, 

-.0397.9, 

*68.52600 •«« . 

-.51909, 

.39919, 

-.06782, 

1 .02460 , 

.04128, 

88.99000,... 

-3.13400, 

- • 0 ° 9 2 ? , 

.19753, 

-2.21570, 

-.03269, 

4.53950*... 

-.13282, 

• C 1156, 

-9.57L20, 

104.5300b , 

.65980, 

— .76180,... 

.80208, 

-.0620.3, 

-.11783, 

1.25900, 

-.04004, 

-25.28500,... 

1.46380. 

.04164 . 

• C25 03, 

-.24815, 

1.37720, 

— 10.39600 t ... . 

.15877, 

.3 C l 48, 

-.00985. 

.05069., 

-12.17000, 

-1.12110,... 

• 066 C 3- « 

-.00457, 

-.01928, 

.11299, 

-.13019, 

-.54573. ... 

.07672, 

-.00979. 

-.02564, 

• 1 3 C 2 7 , 

-4.03850, 

-1.15200, ... 

.54915, 

— . 0 2 26 C . 

-.00615, 

• 0661 2 , 

.47812, 

-H.24900,... 

-.37393, 

-.03837, 

-.02850, 

.10434, 

-.76855, 

1.51880*... 

-.23045, 

-.Cl 450, 

-.03038, 

.08424 , 

-2.25450* 

2.39780/ 

OATA C (REC«2,J,K) ,J=1 

,b) ,K=13, 

1^)/ • • • 

• 


.81469, 

.0 0625 , 

-.00957, 

• 0 0 £> 1 8 % 

-.87656, 

-7 • 3 6 *+ 4 G t •• • 

-.1191 7, 

.91331 , 

. .03352, 

• • uSl bO • 

,33515, 

1*1371C«**» 

.29823, 

-.09182, 

-.90037, 

.00137* 

-.23342, 

-6#35<*9Q/ 


COMMENT 

COMMENT PANELS FREQUENCIES <IN HERTZ) AND DAMPING 
DATA <"FHZ< 1 , J ) , jr 1 ,1S >/. . . 

C. 0b*?26 « 3.12614* , 0.192507 , 0.22bbl2 «... 
0.261074 • 0.32529b , 0.333713 « 0.351529 «... 

0.378027 • 0.3937*8 * 0.410395 • C. 443883 «... 

3.453118 • 3.478315 « 0.501820 / 

DATA <«FMZC2t J),J=1 «15 )✓... 

C. 065226 « 0.126145 . 0.192507 « 0.226612 «... 
9.261074 • 0.329296 « 0.33.371 3 • 0.351 529 «... 

S ".379C27 • 0.393748 , 0.4 J0395 « 0.443683 «... 

0.453118 « 0.473315 * 0.501820 / 

OATA <7M1«J)«J=1 ,15)/15«0.004/ 

DATA < ZM 2 » J ) « J- 1 « 1 5 > / 1 5* 0 • C 0 4 / 

DO L 2 K=l,l e 

v WFU .KlrwFMZtl ,«)*2.*PIF 

Wf «2,*0=WFM7«?,*« >*2.«Plt 
L? .. CONTINUE 
COMMENT 

COMMENT COMPUTATION OF P8 AND PS VLCT08S 



. 

GO TO MFNO 



M2.. 

T RN20N = 1 . 0 




' GO TO MFNO 


. 

M3.. 

TRN30N=1 .0 
GO TO MENO 



M4 • . 

TRN20N= 1.0 

t 

S APSUN=1 .0 


GO TO MEND 



M5 • • 

GO TO MENO 


■ 

M . . 

ME NO . . 
COMMENT 

SCANON=t *C 
CONTINUE 

* 

SLE W0N = 0. 0 S BOX ON = 1.0 


COMMENT 

C0«NIM 

COM*rNT **••** R WH CONTROLLER PARAMETERS *•« 
COMMENT RATf ♦POSITION CONTROLLER PARAMETERS 


CONSTANT NX =107.48287 
CONSTANT KRPX = 49.9M418 
CONSENT RU. SPEED CONTROL 
CONSTANT .IKW = 0. 08681 17 
K-=0.C0031 77 
KHU=32.793G6G 
VLTM=900. 

WR W X I = 0 • 
WRUYl=0. 

WKW71 =0. 

»•♦♦ RyH 


• KT = 2.1016216 

• KRP Y s 49.961418 
LOOP PARAMETERS AND 


« K 1 Ry = C .0 3 


CONTROLLER PARAMETER END 


SENSOR PARAMETERS 


CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
CONSTANT 
COMMENT •• 

CO-PENT ♦•♦.♦♦♦•**.• 

COMMENT CELESTIAL 
CONSTANT T A US = 1.0 
CO-MENT SA AND SP CHARACTERISTICS 
COMMENT SOLAR ARRAY SELECTION 
IF«SALOCK.l.T.0.9»SASLf W = 0.3 
IMSAPSUN.LT.P.9IG0 TO SA1 
SALOCK=C. 

SASLEW=0.0 
SAl ..CON? INUf 
CONSTANT A1SLBG= 

CONSTANT A2SLBC.= 

COMMENT 

COmTNT SP CONTROL 


0.00 

0.00 


A 1 SLE N= 
A2SLIN= 


SO .00 
50 .00 


LOOP PARAMETERS 


CONSTANT 

constant 

CONSTANT 

CONSTANT 


KCL 
HCL 
HS A 1 
USAl 


=1500. 

= 148.50 
=15.00 • 

= 83.05 


AND 1C*S 
KCN =1500. 

BCN = 220.25 
KSA2 =15.00 
BSA2 = 63.05 


COMMENT SCAN SELECTION 
IF ISLE WON. GT .0 . 9>BOXON=0 . 0 
IF ISCANON.3T • 0 • 9 ) G D TO SCI 
SirwON=?.P 
PCXON =0.0 
Sd..CONT INUI 
CONSTANT CLRATTsO.. 1 7453 
ClSLPG= C.00 
C N S L h G.= ICC. CO 
CLPXrilr O.CP 
CLBXB2=15.G0 
CNB»BG=l 1 .00 
BCXEND=fc0.C0 


CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 


CNR ATE =0.01 7453 
CLSLENs 30.00 


CNSLEN= 
CLP XF 0 = 


130.00 

30.00 


• CNBXFQs 15. OC 


« KZ =107.48287 
• KRP2=49.961418 

IC*S 


• A 1 SLRT =0.00438 

• ATSLRT =0.00438 


ORIGINAL PAGE tt 
OF POOR QUALITY 


• CLBXWD=9.0 
« CNBXWD=2.00 



1 


C-4 


1 

• 4 • • 
4 

• 

4 

SI T SI PUL AT ION T IP! 

• 

• 

4 

CONSTANT 

4 

f 1 NT 1 P = 1N9.9N- 

* 

COPPf NT 

4 

4 4 4 4 


4 


4 

4 ' 

SIT 1 M 1 I AL 

CONOl T IONS 


• 

4 

AXIS! 

AXIS./ 

AX I S3 

• 

OISTUKOAMCI TOKQUIS 


constant 

T D ! S T t = J . 0 • 

TIHST? = .rC‘M -<*• 

TllfSTJsO.O 

CC*"f M • 

4 

INITIAL K A 1 1 S 

• 


CONSTANT 

K ATI 1 1 ■ , 

HAH Cl =0.0 . • 

RATt M=0.Q 

nO^MINT • 

4 

INITIAL AM-HLAK 

l-USl T ION 1 KROKS 


CONSTANT 

ALP Mil" 0.0 , 

ALIM? 1 =0.0 « 

ALPHMSU.O 

COPPf NT • 

TORN COPPANOS 



CONSTANT 

T R N T I P = ;* S 0 . V 0 




CONST AMT TRNT KT =»*. H’bbM -‘'t IRN.'R l = * . JbJ .VI - ' » T R N.'R T =»« . 7 Shb* f -S 


SOIAK ARRAY 


SOIAK ARRAY CO-PANOS 

SATSUN-l.O . SAlOfR=I.U » SASLIW-1.0 


SCAN I'l A T T OKP 


SCAN PLATT OK- COPPANTS 

S C A N 0 N : I •0 • SL»VON-J.i » HOXON =1.0 


CPPPINT 

1RN10N= -.0 % TKN;'OV = iI.O % TKNJON = C.O 
fArsi«N=C.O S sc anon = 0.0 
— ANN !*R = -AM IIV ♦ . . I- 

r o TO <P| *Pr«H.A,P^,N*-,-i,)« PANNt'R 
-1.. 1KN10N=l.P 




CONSTANT 
CO--* NT 


CONSTANT 
COPPI NT 


• • • 



C-7 


R ATE1 

Tei 

TPCSRX 
T RVX 
VMX 
URUEX 
IWRWEX 

urvcx 

W RWX 

GUI 

GM2 

GM3 

GMA 

TH1 

TH2 

TH3 

THA 

angm 
C0MSA1 
COMCL 
NXOTt R 

PRCPAK 10 

COFRR1, 
ALPMi , 
RATC1 • 
TP1 . 
TDFSRX, 
TRUX • 
VMX « 
WRUCX. « 
WRWX , 
GM1 « 
GM? • 
0*3 • 

GMA , 
TRNCOM . 
ANGM , 


RANGE ... 

TRNCOl 
AIPH1 
ERROR 1 
COCRR1 
RATE 1 
TP 1 

TPCSRX 

TRVX 

VMX 

W RUEX 

TyRyrx 

jrvcx 

URyx 

SMI 

G*2 

GM3 

GMA 


R ATF 2 , 

RATE3 , 

TB2 ♦ 

TBS « 

TOESRY « 

TDCSRZ, 

TRy Y • 

TRyZ , 

V MY , 

VMZ « 

yRUEY • 

y r y e z « 

i yR y e y • 

iyRycz* 

URyCY « 

yRycz « 

yRUY , 

VRVZ • 

C T A 1 1 , 

ETA21 , 

ETA12 • 

ETA22 • 

f T A 1 3 « 

ETA23 • 

etaia • 

E T A2A t 

ETAlb • 

ETA25 • 

FTA16 ,, 

ETA26 • 

ETA17 « 

ETA27 • 

ETA1B , 

FTA28 • 

OERPOY* 

NXQ1NT • 

C0MSA2, 

NXQOVN, 

COMCN • 

NXQOE R « 


COE K R2 t 

COERR 3, 

ALPH2 , 

ALPH3 , 

RATE 2 • 

RATE3 « 

TB2 • 

TB3 , 

TOESRY, 

TOESRZ, 

TRy Y , 

TRyZ , 

w my • 

VMZ • 

WRUCY • 

yRycz , 

y Ry v , 

yRy? , 

TH1 , 

COMSAl, 

TM2 , 

C0MSA2, 

TH3 , 

COMCL , 

THA , 

COMCN , 

FTA11 , 

ETA12 , 

ETA13 , 

ETAIA , 

FTA15 , 

FTAlb , 

FTA17 ., 

ETA18 

TRNC02, 

TRNC03, 

ALPH2 , 

ALPH3 , 

FRROR2, 

ERR0R3, 

COERR?, 

C0LRR3, 

RATE? , 

RATE 3 « 

TB2 , 

TBS • 

TOESRY, 

TOESRZ, 

TRyx , 

TRyz , 

VNY , 

VMZ , 

y RyE y , 

VR VF Z , 

I yRy e y « 

1VRVEZ, 

VRVCY , 

yRycz , 

y RU Y , 

VRVZ , 

E T A 1 1 , 

CTA21 , 

ETA1? , 

ETA22 , 

ETA13 , 

ETA23 , 

etaia , 

E.TA2A , 






*°Otf 
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COMMENT 

COMMENT SET UP PLOTTING INFORMATION 
NOPLOT = .FALSE. 

11 CONSTANT PLTNBRsO «... 

INTRV1=10.0 .... 

NEW INT =100.C .... 

INTRV2=30.0 

SUITCH=NFWINT 

COMMENT 
NC=3 % NF=? 

NX!JOYN = 0. 

NXGDf R -0 • 

NX0TER=0. 

NXQ!NT=NXaiNT«l. 

COMMENT 

CALL NMOYFL <NC«H. MP.MS.PB.PS. G.PI.NF.F.E IGfREC.RF.WF.2F) 

or.nuG 

LNO 

COMMENT 

DYNAMIC 

VARIABLE T = T IN1TL 

CALL 72-0N<*AFTRWG».T.1.0,r,0) 

OFRPOY=NXQDFR-NXOERL 

NXOERLsNXOnKR 

NXQOYN=NXQOYN*I. 

IF <NOPLOT) GO TO NOSAVE 
IF IY.GT. SWITCH) GO TO 1NT2 
YFSPLT=PLTNOR« INTRV1 
GO TO SI 1 
! NT 2 • • CONTINUE 
CHrCKP=PLTNPR. INTRV1 

IF CCHECKP.Gr .NFWINT) NEW I NT =NF. WINT.INTRV2 
YESPLT = Nf W INT 
SL 1 • • CONTINUE 
TlNEP = T..0r,01 

IF <TIMFP.LT. YESPLT) GO TO NOSAVl 
WRTTr <b.F?) PLTNRK.T 
WRITE Ib.F?) P2tPI.P2.P3 
WRITE <f.F2) GH < 1) «GH<2) .GM(3).GM(4) 

DO SL 2 Jsl,? 

LSO 

00 SL2 Mil. 3 

WRITE <ft.F2> <ET<J.L»M>.Mrl,5> 

-*■ F 2 •• FORMAT <1X,5G1A.R) 

L=L.S 


S L 2. • CONTINUE 

WRITE <G.F3) ,<<CK<N.NN).NN=1.3).NsI . 2) 
F3. • FORMAT <1*.3G1A.P> 
PLTNPRsPLTNHR^I 

IF <T.GT. SWITCH) NEWINT=NEWINT. JNTRV2 


NOSAVC.. CONTINUE 
IF<T.GT.riNTIM)GO TO FIN 



OUTPUT 200 

TRNCOI. TRNC02* 

TRNC03 • 

• • 

ALPM1 . ALPH? , 

ALPH3 . . 

• • 

ERROR I « ERR0R2. 

ERR0R3. 

• • 

COERR 1 * C0ERR2 * 

C0ERR3. 

• • 


CLSL0N=PULSE I C LSLBG*50 000 • • «CLSL£N~CLSLBG> *T >*SLEWON 
CNSLON=PULSC (CNSLB6* 50 COO • * CCNSLEN-CNSLB6) *T J ‘SLEWON 
IFtT.GT.BOXENO>BOXON=0.0 

CLBXON* IPULSE (CL8XHl*CL&XFQ«CLBXV0«T)... 

-PULSE <CLBXB2»CLBXFQ*CLBXWD*T>>*B0X0N 
CNBXON= PULSE<CNBXPG*CN8XFQ«CNBXW0«T»*B0X0N 
CCMCL=INT£G< tCLSLON^CLBXON) *CLR ATE • 0 . » 

COMCNslNTEG( CCNSLON^CNBXON > «CNR ATE *0. ) 

COMMENT SOLAR ARRAYS CONTROLLER 
A1 SLON=PULSE <A1SLBG*5COOO.* CA1SLEN-A1SLBG) «T»*SASLEW 
Cf-SAl = t-ALPH2«SAPSUNMNTEGUlSL0N*AlSLRT*C.»>*<l.-SAL0CK> 
A?SLONsPULSE »A ?$LP 0*5 0000. • IA2SLEN-A2SL8G* «T)*SASLCW 
CC"SA?s<- 4LPH2*SAPSUN4INTEG<A2SLON«A2$LRT*O.J>*<1 .-SAL CCK I 
COMMENT SOLAR ARRAY HTNGC TORQUES 
THls-KSA 1* <GM1-C0MSA1 )-eS-l*GM10 
TH2 = -KSA2* (GM2-C0MSA2,)-PSA2*GM2D 
COMMENT SCAN PLATFORM HINGE TORQUES 
TH3 = -KCL««GM3-CO**CL >-BCL*GM30 
TH4S-KCN* <GM4-C0MCN> BCN-GM40 
COMMENT 

COMMENT GENERATE TURN RAMP INPUT SIGNALS 
TRNTRU = PULSE <0.0*5C000.0*TRNTIM*T> 
TRNCOl?lNTEGnRNiON*TRNlRT«TRNTRU«0.) 

TRNC02=INtrC (TRN20N*T*;N2RT*TRNTRU«0.) 
TRNC03=INTtG<TRN30N»TRN3RT*TRNTRU*0.> 

TRNC0M=TKNC01 ♦ TRNC0?^TRNC03 
COMMENT CELESTIAL SENSORS 

THfTAl=Kf ALPLETAUS.ALPH1 «0. I 
THE TA2=R r ALPL<TAUS*ALPH?*0. > 

THETA3=PEALPL«TAUS ,ALPH3*0.» 

ERROR 1 rTRNCOl-THETAl 
ERR0 R?sTRNCO2-THETA2 
ERROR 3= T RNC03-TMET AS 

COMMENT •***•******••***•***•*••*•*•*••***•**•****♦***♦ 
COMMENT **♦*«« RUH CONTROLLER OYNAMICS SECTION •****♦ 
COMMENT RATE*POSITION CONTROLLER 
CrrRRl=ERRORl-KRPX*RATCl 
COERR. ‘* = EPR OR ?-KRPY*RATE? 

C0FRR?=ERR0R3-KRPZ*RATF7 

Tr»CSRX=KX«C0ERRl4jRW*i-RATr3«TACHWY>RATE2*TACHWZ> 
T0ESRY=MY*C0ERR2*URW*<*RATE3*TACHWX-RATEI*TACHUZ> 
T0ESRZ=K?*C0ERR3*JRW*<-R ATE?«TACHUX*RATE1«TACHWY» 
WRWCXsLIMINT E-TOFSRX/ JR W,0. *-2 30 • 3834 5 •♦230. 38345> 

WPWCY = LTMINT<-TOESRY/JRW,0.«-230.3e345 >♦230. 383451 
WPWCZ=LI M INT<-TOESRZ/JRW«0. *-230. 3P345* ♦230. 383451 
COMMENT RW SPEED CONTROL LOOP MODEL 
WRWrx=WKUCY-WRWX 
WRWEY=WRWCY-WPWY 
WRWEZSWRUC7-WRWZ 

lWRUtX=Ll**INT «KIRW*WRWEX*0. *-27. 447 *♦27.447) 
lWRWEY=LIMINT<KlRW«WRWEY*e.,-27.447*^27.447) 

IWRyEZsLIM INTI KIRW«WRWEZ*0. *-27.447* ♦27.4471 
VCX=MPU*<WRwrx ♦IWRUFX* 

VCY=KKW*<WRWTY ♦IWRWEY) 

VCZ=M»W*<WRWEZ^1WRWEZ) 

VLX sPOUNflE - VLIM**VLIM*VCXI 
VLY=90L'NDI-VLIM*^VLIM*VCY> 
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TH1 * 

F.TA15 « 

ETA25 • 



TH2 • 

ETAlb » 

ETA26 • 



TH3 * 

ETA17 • 

ETA27 • 



TH4 • 

CTA18 « 

ETA28 • 

. • 


ANGM • 

OERPOY* 

NXOINTf 



COMSA1 • 

COMSA2« 

NXQDYN, 



COHtL « 

CC"CN « 

NXODER * 




NXQTER 

DERIVATIVE CIO 

CINTERVAL CI=0.05 
XFRROR ... 


GMlO sl.E-R 

« GN2P 

=1 .F-8 

• GM1 

= 1 

.C-8 

* GM 2 

=1.F-B 

• • • • 

GM3D = 1.£ -8 

• GM4D 

-l.E-8 

• G*3 

= 1 

• F -8 

*GN4 

= 1 .E-8 

• • • • 

PATE 1 =l.F-6 

•RATE? 

=l.E-b 

* R A T F 3 

= 1 

.f -b 




MERROR ... 









GMIO =l.f-8 

• GM20 

=1 .E-B 

«6Nl 

= 1 

• L - 8 

• GM2 

-l.E-8 

• • • • 

GM30 =l.E-8 

• GH40 

=1 .F-B 

«GM3 

= 1 

• E-8 

* GM 4 

= 1 .F-8 

t • • • 

RATE1 = l.F-b 

♦RATE? 

=1 .t-b 

«R A T E 3 

= 1 

.E-b 





NOSORT 

NXQOCR=NXO0ER*1. 


GM0C11=GM1D S 

GM< 1 1 

= C,M 1 


C-M0C2 

1 = G**2P » 

GNC21 

= GM2 

* 

GNDC31=GM30 * 

GM < 3 1 

— GM 3 


GM0(»)=GM4D \ 

GM(41 

= (,M4 


FTOC 1 

1 ) sf TOl 5 ' 


r t < i ♦ 1 1 

= E T A1 1 

ETOU 

2> =r T012 


FTC 1*21 

= F T. A 1 2 

FT041 

31 = FT013 


ETC 1*31 

= CT A13 

ETOC1 

4 1 = F. T D 1 4 

< 

FT C 1 « 4 1 

= E T A 1 4 

L TO Cl 

51 =ET015 

X 

FTC1.51 

= C T A 1 5 

FTO C 1 

b> =F TOl b 

X 

F T C 1 • b 1 

= E TAlb 

FTO C 1 

71 =ET017 

X 

L T C 1 • 7 1 

= E T A I 7 

ETDC1 

ei =ETD18 

« 

ET C 1*81 

= E T A 1 8 

noci 

91 =TTP19 

1 

FTC1.91 

= F T A l 9 

ETOU 

10!=FTD11C 

t 

FTC1.101 

= C TAM J 

ETOU 

11 1=F TD1 11 

X 

ETCltlll 

= E T A 1 1 1 

ETOC1 

1? l=t TOl 12 

X 

ETC 1*121 

= E T A 1 12 

ETOU 

13 l=rT01 13 

X 

FTC 1*13) 

=F T A 1 1 3 

f. T D C 1 

14 1=FT01 14 

s 

ETC1 *141 

= E T A1 1 4 

ETOU 

1 51 = F T 01 15 

$ 

FT C 1 • 15 1 

= ETA1 15 

ETDC2 

11 =F T 02 1 

% 

E T C 2 * 1 1 

= E TA21 

FTDC2 

21 =CT022 

t 

F T C 2*2 1 

= F T A 22 

ET0C2 

31 =FT023 

X 

E TC 2*3 1 

~E T A 2 3 

FTO <2 

41 =FT024 

X 

E T C 2 *4 1 

= 1 T A24 

FT0C2 

51 = r T025 


ET C 2*51 

=FTA25 

FTDC2 

bl =ET02b 

t 

FT C 2 «b 1 

=ETA2b 

FTO 12 

71 =ETD27 

* 

CTC2.71 

= E T A27 

ET0C2 

81 =ETD28 

X 

F T C 2 • ft 1 

= T T A28 

CT0C2 

91 =f T 029 

«, 

F T C 2* 9 1 

=F T A 29 

FT0C2 

101=F TO210 


FTC2.1G1 

=ETA21u 

FT0C2 

11 1=1 T0211 


ETC2*11! 

= F T A2 1 1 

ET0C2 

121 = r TD212 


F T C2 *1 21 

=ETA212 

ET0C2 

1 3 1=E T 02 13 


FTC2.131 

= f T A 2 1 3 

FTO C 2 

14 )=FT0214 

X 

CTC2«141 

=FTA214 

FT 0 <2 

151 = FTD215 

X 

FTC2*151 

=f TA215 


V3C11=RATE1 * V0C21=RATE2 X V0C31=RATt3 X ANGM=HM 
COMMENT 

CONSENT SCAN PLATFORM CONTROLLER 


/ 


/ 

/ 
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ETDll = 1 NTE G (E TOO < 1 

*n 

• 0.) 

s 

ETA11 = 1 NTEG <ETD 1 1 

• 0. > 

ETD12 -INTFG CETDD.l 1 

*2) 

• 0.) 

* 

FT A 12 =INTEGtETD12 

• 0. ) 

PT013 =INTEG IFTDOE 1 

• 3) 

• 0.) 

s 

ETA13 = I NT EG <ETO 1 3 

• 0. ) 

ET014 - I NTEG <E TOD< 1 

• 4) 

• 0.) 

s 

ETA14 = I NTEG <E TO 1 4 

• 0. ) 

CTO 15 =INTCG <E TOO< 1 

• 5 > 

• 0.) 

* 

ETA15 sINTEG 1ETD15 

• 0. ) 

CTOlft s I \’TF6 <E TOO< l 

«6) 

• o.» 

t 

ETA1G = INTEG <ET01 6 

• 0. ) 

CT017 = I NTEG <C TOO < 1 

*7) 

• 0.) 

$ 

ETA17 = I NTEG <E TO l 7 

• 0. ) 

CTD18 = I NTFG <E TOO < 1 

• 8) 

• 0.) 

t 

ETA18 =1 NT EG <E TO 1 8 

• 0. ) 

FTD19 = I NTEG <£ TDD< 1 

• 9) 

• 0.) 

s 

ETA19 =INTEG<ETD19 

• 0. > 

CTO110=INTFGCCTO0«l 

*10) 

• 0.) 

$ 

etaiio=intfg<etoiio 

• 0. > 

ETD111=1NTEGIET0DC1 

• 11 > 

• 0. ) 

s 

ETA111 = 1NTF.G(ETD111 

• 0. ) 

rT0112 = JVTCG fETOOn 

*12) 

• 0.) 

5 

tTA112=INTLG<LTD112 

• 0. ) 

) Tnm = INTFG CETODt 1 

• 13) 

• 0.) 

f 

ETAU3 = 1NTEGCETD113 

• 0. ) 

r TDll42lNTEGCETDD<l 

• 1* > 

• 0.) 

t 

ETAU4 = lNT»G(ET0li4 

• 0. ) 

ETD115 = I\TIG tETDDU 

• 15) 

• 0.) 

% 

ETAllb = lNTEG CET0115 

• 0. ) 

F T 0 ? 1 slNT* G lETDCU? 

• 1 > 

• 0. ) 

% 

ETA21 =1NTEG (ETD21 , 

• 0. ) 

FT022 = I NT r G (F TOP < 2 

• 2) 

• 0. > 

$ 

ETA22 sINTEG (ETD22 

• 0. ) 

F T 0 2 3 = I NTEG (E TOO < 2 

• 3 1 

• 0 •> 

% 

ETA23 =1NT C G(LT023 

• 0. ) 

CT024 s I NTE Git TOO i 2 

*4 ) 

• 0.) 

% 

LTA24 =T NT EG 1ETD24 

*0. ) 

ETD25 =INTCGCETOD(? 

• 5) 

• 0.) 

i 

ETA25 = 1 NT EG (E TO 25 

• 0. ) 

F T02 t > 2 1 NTE G <E TOO ( 2 

• 5) 

• 0. > 

* 

ETA2u '= 1 NTEG (FT026 

• 9. ) 

FT027 2 I NTT G iE TDO < 2 

• 7) 

• 0. ) 

s 

ETA27 =1NTFG IET027 

• 9. > 

ET028 r I NTEG IE TOO < 2 

• 8) 

• 0. ) 

•l 

ETA28 =1NTEG (ET028 

• 0. ) 

r T02 Q 2 I NTFG (E TDO ( 2 

.9) 

• G.) 

$ 

ETA29 2 I NTEG (ET029 

• a. ) 

f T021C=I\'T r G<ET0O<? 

• io ) 

*0 . ) 

& 

ETA210=INTEGCET0210 

• 0. ) 

t T0211s|NTCG<CT00<2 

• u ) 

*C. ) 

$ 

ETA211=INTEG<ET0211 

• 0. > 

CT02!?=INTrG«rT00(? 

• 12) 

• 0.) 

t 

CTA212=1NTEG<ET0212 

• 0. ) 

IT 02 15 = 1 NTFG *.ET0tM ? 

• 13) 

• 0.) 

% 

ETA213=lNTEGtET0213 

• 0. ) 

FT0214=INTf GCETDOI'* 

*14 ) 

• 0.) 

s 

FTA214=INTEGIET0214 

• 0. ) 

ET0215=INTEGCET00C2 

.15) 

• 0.) 

s 

ETA215=INTEGCET0215 

• 0. > 


COMMENT ****** FLEX DYNAMICS END ****** 
CO?1fr\'T ••*•.****•**.•**•**.*•****•****** 

FND 
f MO 

r m d 

TERMINAL 

FIN.. CONTINUE 
Ot MUG 

NXOTF R=NXOTFR*l . 

WRITE <F«F?)T «NXQIMT«NXOOYN.NXQOER*NXQTER 

FNt; 

END 


r '■ 

t 


N G , 


l * C«M 


* 
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*aZ = eouNOI-VLlH,*VLlH*VCZ>, 

VHXsSORT <A8S(VLX)) • SI GNU . *VCX> 

VHY^SQRT <A0S <VLY>>*SIGN<1.*VCY> 

VNZsSORf <ABS(VLZ)>*SIGNU**VCZ> 

TRWXSKH* VHX* ABSIVHXJ 
4RUtSKM*VNY«ABSCVNT) 

TRm2sKH*VH2*ABS<VHZ) 

WRWX = INTCG(TR«X/JRW*WRWXI > 

WRWY=INT£GITR«Y/JRU*WRWYI) 

WRWZ = INTrG«TRliZ/JRW»WRMZI I 
T ACHWX=URWX-RATE1 
TAf HWY SWRWY-RATE2 
TACHWZ=WRWZ-RATE3 

TPX=*TRWX-JP«* <-RATE3*TACHWY*RATC2*TACHyZ» * FBX=Q.O 
TPY=-TKWY-JRW* <*RATf 3*TACHWX-RATE1«TACHWZ> * FBY=0.Q 
TPZs-TRWZ-JRW* <-RATE2*TACHWX*RkTE 1*TACHWY> S FPZsQ.Q 
COMMF NT ****** RWH CONTROLLER OYNAHICS ENO ****** 

COHHcNT 

COXXFNT .##*..********♦***********•* 

C OHMF NT ••**«• FLEX OYNAHICS ****** 

COHHCNT FORCfS ANO TORQUES ACTING ON BUS ANO HINGES , . 
TB1=TPX*TD1ST1 4 TBU» = TB1 

TB?=Tpy*TnrSTZ S TP(2I=TP2 

TP3=TbZ*TDIST3 * TBC3I=TB3 

Ff>U> = FP< 

FP<2»=F»>V 
F H(3>spH7 

THUJrihl _ 

TH<2>=T»’2 

TH<? >sTH3 

TMA» = THA 

COMMENT 

COHHCNT SOLVE FOR SYSTEH ACCELERATIONS 

CALL HR ATE <NC» TH,TR*TS»FB*FS*TF*FF . GH, GHO« GHOO *ET • ET D t ... 
VOtUOOT »f TDD*HH*CK ) 

COHHENT 

COHHCNT SYSTEM RATES A NO POSITIONS 
COHHCNT . 

RATEl =INTEG<WOOT«l»*RATClI> 

RATE? rINTCG(WDOT<2>»RATr2n 
RAtrj =INTCG(UOOT < 3) «R AT F3 I > 

ALPHI =INTEG<RATE1*ALPH1I) 

ALPH? =INTF.G*RATE?*ALPH2I) 

ALPH3 slNTEC <R ATE ? • ALPH3I > 

P?0«Pin*P?O f P3DsHCK <HCK f PZ*Pl »P 2 *P3«R ATE I *R A TE 2 * R A TF 3 > 

9 PZ =lNTfGCPZO*l *00C I 

PI slVTro <P1D*0.000 * 

P? =1NT1GIP?0*0.000) 

P? slMTFG <P3O*r..0QC> 

* f.HlP =INTEG(MOOTCA)*0.» 

GH20 =INTf GTUOOT *0 « ) 

GH 30 =lNTEG<ynOT < M ,0. I 

G"AQ = I NT r G IMOOT < 7 ) « C • > 

C.H1 - I NTEG (GRID* 9« ) 

0H2 =INTf GCG**20» J«» 

GH? =lNTEjtGH30*0«) 

GHA =INTFG«GHAO*u*> 




I 
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I 

COHHENT CALCULATE THE S/C CENTER OF HASS 
DO CHI K = 1.3 

CMSCLV<K)=MP<7I*CNBLV<K)4NS<1«7)*CHSLV(1*K> ♦MSI2«7>*CMSLV... 
' # | C2,K|-»HS<3t7>«CHSLVl3*K> 

* CMSCLV<K>=CMSCLV<K >/ CHB < 7 » ♦HS < 1 • 7> *HS I 2* 7 » ♦ MSI 3 • 7> > 

CHI,. CONTINUE 
WRITE ( 6 « F 3 ) CHSCLV 

COHHENT COMPUTATION OF PB AND PS VECTORS AND CHB TO ENG VECTORS 
00 LI Ksl *3 

RJTC1<K)=JTC1LV<K>-CHBLV<K> S R JT C2( K > =JTC2L V < K > -CNBLV < K > 
CJTCI(K)-JTC1LV<K)-CMSCLVIK) $ S JT C2 ( K ) = JTC2L V < K I-CHSCL V I K ) 
PP<l.K )=HINGLV«1.K>-C,PBLV<K» 

PH<?.K>=HlNGLV<2.K)-CHBLVtX> 

PP<3*K»=HINGLV<3*K)-CMBLV<X) 

PSH*l*K)sHINGLVCl*K>- CMSLVC1*K> . 

PS <2,2 .K » SH1NGLV <2.K >-C«SLV< 2.K > 

PS<3«3.K)sHINGLV < 5.K l-CMSLV < 3*K » 

LI.. CONTINUE 
PICKUP 
INITIAL 

COMMENT ♦«**•*.*.**.*.***«.****••.••*.*******•********•********** 


4 


C -s 


OEFINE THE BUS FIXED AXES AS FOLLOWS: 

A x I S 3 = Z * ION ENGINE PLUME EXHAUST AXIS > 
AXIS? -Vi SOLAR ARRAY AXIS > 

A X I S 1 = X 


SELECT THE MANEUVER BY ITS NUMBER 


1 - TURN 

ENTIRE 

SPACECRAFT 

ABOUT 

AXIS1 

2 - TURN 

ENTIRE 

SPACECRAF T 

ABOUT 

AXIS? 

3 - TURN 

ENTIRE 

SPACECRAFT 

ABOUT 

AX IS 3 


4 - TURN BUS ONLY ABOUT AXIS2 

WHILE KEEPING ARRAYS ON THE SUN 

5 - ACQUISITION OF CEL REFERENCES 
ft - SCAN PLATFORM SLEW OR ErOXSCAN 


HANEUV = 1.0 


CONST 6NT 
COHMfNT 




0 


FILMED 



